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自从 六 十 年 代 线 性 多 变量 控制 理论 问世 以 来 ， 如 同一 切 
新 的 理论 出 现 一 样 ， 它 与 原 有 的 经 典 控 制 理 论 之 间 产 生 了 许 
多 的 不 协调 。 这 里 既 有 这 两 方面 理论 本 身 的 原因 , 也 有 一 些 
难以 避免 的 社会 因素 。 这 种 不 协调 在 一 段 时 间 里 甚至 发 展 到 
如 本 书 作 者 所 说 的 “ 裂 症 ”的 程度 ,这 种 状况 的 形成 是 不 奇怪 
的 甚至 可 以 说 是 自然 的 . 

人 们 常 说 的 经 典 控制 理论 ， 一 般 是 指 本 世纪 五 十 年 代 前 
创立 的 单 回 路 调节 理论 .这 一 理论 的 出 现 首先 不 是 由 于 严格 
的 具有 完整 体系 的 一 套数 学 而 是 由 于 工程 与 物理 的 推动 ， 数 
学 在 这 里 是 作为 论证 与 推演 控制 思想 的 工具 ， 而 控制 本 身 的 
结论 及 其 表述 也 远 未 数学 化 .这 一 理论 的 两 个 主要 方法 ,频率 
法 与 根 轨迹 法 也 都 具有 相当 强 的 工程 与 物理 意义 ,而且 这 一 
理论 还 能 针对 实际 的 工程 要 求 ( 例 如 过 渡 过 程 时 间 , 超 调 量 ， 
振荡 次 数 等 ) 提 出 一 些 实际 的 解决 办 法 。 正 由 于 此 ,迄今 在 相 
对 简单 的 控制 系统 ,如 单 回路 调节 系统 与 随 动 系统 设计 中 ,经 
典 控制 理论 所 提供 的 方法 仍然 是 主要 的 也 是 行 之 有 效 的 . 随 
着 工程 系统 的 复杂 化 或 控制 要 求 的 提高 已 不 能 对 系统 采用 过 
于 简化 的 模型 ,而 且 将 控制 理论 向 其 它 应 用 领域 扩展 ,例如 经 
济 系统 ,环境 系统 等 ,都 使 经 典 控 制 理论 所 提供 的 方法 磁 到 了 
难以 克服 的 困难 ， 这 是 由 于 这 些 系统 从 本 质 上 说 来 是 多 变量 
多 回路 的 系统 ， 而 要 将 它们 强行 纳入 单 回路 系统 的 框框 是 很 
难 奏效 的 。 这 表明 经 典 控制 理论 不 作 适 当 的 改造 与 发 展 使 其 
适应 多 变量 多 回路 的 特点 ， 就 不 可 能 在 今天 更 为 复杂 的 控制 
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系统 中 发 挥 作用 。 而 且 为 了 适应 新 的 复杂 系统 的 要 求 , 没有 
深刻 的 理论 思维 (包括 深刻 的 数学 理论 思维 ) 也 是 很 难 取得 成 
功 的 . 

六 十 年 代 兴起 的 线性 多 变量 控制 理论 ， 在 描述 方法 上 一 
开始 就 采用 了 状态 空间 模式 ， 这 种 模式 把 系统 归结 为 三 个 或 
四 个 具有 适当 阶 次 的 矩阵 .由 于 这 种 模式 的 一 般 性 ,不 少 控制 
问题 就 归结 为 这 几 个 甜 阵 或 它们 所 代表 的 映射 应 具有 的 要 求 
和 满足 的 关系 .这 样 控制 系统 的 一 些 问题 经 过 转化 就 成 为 比 
较 纯化 的 数学 问题 特别 是 线性 代数 问题 ,例如 可 控 性 、 任 意 配 
置 极点 、 内 模 原理 等 。 这 些 问题 的 解决 已 经 在 实际 系统 分 析 
与 设计 中 起 到 了 作用 ,特别 由 于 计算 机 与 计算 方法 的 进步 ,又 
使 这 一 理论 得 到 了 新 的 动力 一 一 利用 先进 的 计算 手段 来 开拓 
理论 的 应 用 领域 ， 这 样 不 仅 在 理论 上 而 且 在 各 种 计算 机 辅助 
设计 上 都 得 到 了 发 展 。 尽管 如 此 ,由 于 这 种 模式 的 一 般 化 , 它 
与 实际 工程 系统 的 结合 总 有 一 段 距 离 。 对 于 这 种 模式 我 们 也 
难以 将 具体 的 控制 工程 要 求 提 成 数学 便于 接受 的 形式 ， 而 且 
对 于 这 种 模式 由 于 其 物理 工程 背景 的 不 直接 ， 因 而 物理 与 工 
程 的 洞察 力 在 这 里 也 难以 起 到 作用 。 于 是 对 熟悉 工程 与 物理 
的 控制 工程 师 说 来 ， 在 相当 长 的 一 段 时 间 里 总 是 不 乐于 采用 
线性 多 变量 系统 的 理论 来 解决 控制 工程 问题 ， 即 使 在 磁 到 复 
杂 的 系统 ,他 们 也 仍然 想 用 单 回 路 的 经 典 理论 去 套 , 致 使 磁 到 
一 些 难 于 克服 的 困难 . 

这 就 是 二 十 多 年 来 控制 理论 发 展 过程 中 的 现实 ， 面 对 这 
个 现实 一 些 人 不 是 固守 在 经 典 控 制 理论 或 线性 多 变量 系统 理 
论 的 一 个 方面 ,去 谋求 在 已 有 的 框框 内 进行 发 展 ,更 不 是 站 在 
一 方 去 论证 另 一 方 的 局 限 性 ， 而 是 采用 兼 收 二 者 之 长 以 补 二 
者 之 短 将 这 两 方面 结合 起 来 的 做 法 .这 里 译 出 的 an，Postleth- 
waite £3 Alistair, G. J. MacFarlane 合 著 的 《线性 多 变量 反馈 
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系统 分 析 的 复 变 方法 》 一 书 ,就 是 这 样 的 一 本 著作 。 相信 ， 本 
书 的 内 容 无 论 对 于 从 事 控制 工程 还 是 从 事 控制 理论 的 人 来 
讲 , 都 是 会 引起 浓厚 的 兴趣 的 . 

为 了 将 经 典 控制 理论 中 的 两 个 主要 方法 , 频率 法 或 
Nyquist 方法 与 根 轨 线 法 或 Evans 方法 推广 至 多 变量 情形 ,而 
现在 闭环 频率 与 开 环 增益 之 间 的 关系 已 经 是 多 值 化 的 代数 函 
数 关系 ,出 于 问题 本 质 上 的 复杂 性 ,作者 运用 代数 函数 理论 及 
与 之 相关 的 Riemann 曲面 结构 来 阐述 问题 , 得 到 了 比较 清楚 
的 结果 . 虽然 代数 函数 理论 与 Riemann 曲面 的 结构 对 控制 界 
比较 生疏 ， 但 由 于 作者 对 这 些 理论 作 了 比较 形象 的 叙述 并 且 
提供 了 具体 的 计算 方法 ， 就 使 得 这 些 理论 读 起 来 并 不 枯燥 而 
且 可 以 学 会 如 何 实际 应 用 。 实际 上 即使 对 经 典 控 制 理论 中 的 
根 轨迹 方法 ,为 了 使 理论 清楚 严谨 ,引进 代数 函数 与 Riemann 
曲面 也 是 可 取 的 . 

由 于 本 书 是 这 方面 的 第 一 本 著作 ， 原 书写 作 中 有 其 自身 
的 特点 ,加 之 本 人 业务 水 平 所 限 , 失误 醇 漏 之 处 难免 ， 敬 请 读 
者 指正 . 

祝愿 本 书 中 文 版 的 出 版 能 在 我 国 控制 理论 及 其 应 用 的 事 
业 中 ,在 促进 我 国 控制 工程 与 控制 理论 两 方面 的 合作 中 起 到 
应 有 的 作用 . 

北京 大 学 力学 系 黄 琳 
1984. 8. 25 
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在 五 十 年 代 后 期 与 六 十 年 代 初 的 一 段 时 间 里 ， 最 优 控制 
与 最 优 滤波 技术 在 空间 研究 方面 取得 了 巨大 成 功 ， 这 就 自然 
地 导致 将 这 些 方法 应 用 于 世间 各 种 多 变量 工业 过 程 的 努力 . 
这 种 努力 在 很 多 场合 并 不 能 直接 奏效 ， 特 别 比 起 空间 研究 来 
说 ,这 里 所 采用 的 装置 的 模式 是 不 精确 的 ,或 者 用 来 刻 划 受 控 
装置 品 性 的 品质 指标 也 不 明确 。 MA, 基于 直接 应 用 最 优 控 
制 与 最 优 滤波 方法 的 控制 器 ,其 综合 方法 一 般 讲 是 复杂 的 ; 事 
实 上 如 果 要 装 上 一 个 完全 的 Kalman-Bucy 滤波 器 ,由 于 该 滤 
波 器 应 由 受 控 装置 的 模式 与 其 周 图 的 反馈 所 组 成 ， 这 样 其 动 
态 复杂 程度 就 将 与 受 控 对 象 相当 .与 此 相反 ,对 不 少 多 变量 过 
程控 制 问题 来 说 ,真正 需要 的 乃 是 相对 简单 的 控制 器 , 它 一 方 
面 应 能 在 工作 点 附近 实现 镇 定 且 仅 要 求 受 控 装置 有 一 可 用 的 
近似 模式 ， 同 时 它 还 应 能 缓和 由 于 引 人 积 分 作用 后 低频 干扰 
的 影响 。 对 于 惯 于 用 频率 响应 思维 的 工业 工程 师 说 来 , 应 用 
深奥 的 最 优 控制 方法 是 困难 的 ; 这 些 工 程 师 基本 上 依靠 的 是 
物理 的 洞察 能 力 和 例如 用 微 积分 这 些 直接 方法 去 解决 问题 . 
明显 的 是 在 由 Nyquist [1], Bode [2] 与 Evans [3] 所 竟 基 
且 至 今 还 在 不 少 工业 应 用 中 采用 的 经 典 单 回 路 频率 响应 方法 
和 在 空间 应 用 中 发 展 起 来 的 精美 而 又 有 效 的 多 变量 时 间 响 应 
方法 之 间 已 存在 着 巨大 的 裂痕 . 

由 于 这 些 原因 ， 在 六 十 年 代 中 期 慢 慢 开 始 恢复 了 对 频率 
响应 方法 的 关注 ， 在 经 典 频 率 响 应 方法 与 最 优 控制 方法 之 间 
弥合 裂痕 的 重要 的 第 一 步 是 由 Kalman [4] 给 出 的 ， 他 研究 
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了 最 优 性 的 频率 域 特征 。 对 多 变量 系统 发 展 频率 响应 分 析 与 
设计 理论 的 系统 性 的 研究 ,是 从 Rosenbrock [5] 的 开拓 性 论文 
开始 的 ， 于 是 在 已 复苏 了 的 频率 响应 方法 上 开始 了 增长 兴趣 
的 十 年 。 在 这 一 新 起 点 之 前 , 在 多 变量 控制 问题 上 已 经 形成 
了 一 些 直接 方法 ，Boksembom 与 Hood [6] 提出 了 不 相互 作 
用 控制 器 的 思想 ， 他 们 的 方法 是 选择 一 个 级 联 补偿 器 以 使 受 
补偿 系统 总 传递 函数 矩阵 为 对 角形 ， 如 果 求 得 了 这 样 的 补偿 
器 ,控制 器 设计 就 可 以 运用 标准 的 单 回路 设计 方法 来 完成 . 
用 这 种 方法 求 得 的 补偿 矩阵 必然 是 复杂 的 ， 并 且 该 方法 明显 
的 缺 限 乃 是 仅 为 减少 相互 作用 就 要 进行 如 此 宛 烦 的 工作 .这 
一 解决 多 变量 控制 途径 自然 的 延伸 是 考虑 运用 有 理 矩 阵 进行 
标准 矩阵 计算 方法 能 得 到 什么 。 沿 着 这 个 思路 讨论 问题 的 文 
章 有 Golomb 5 Usdin [7], Raymond [8], Kavanagh [9, 10, 
11] 与 Freeman [12,13]. S&iffi, Rosenbrock [14, 15] 却 用 
一 种 语 熟 的 技术 寻求 将 多 变量 问题 化 至 经 典 方法 可 接受 的 形 
R, 从 而 开辟 了 一 个 全 新 的 发 展 路 线 。 在 他 的 逆 Nyquist A 
阵 方法 [14, 15] 中 , 其 要 领 就 在 于 减少 相互 作用 到 可 以 采用 
单 回路 方法 的 程度 而 不 是 去 完全 消除 相互 作用 . Rosenbrock 
的 方法 就 是 基于 部 分 相互 作用 的 特定 判 据 一 一 对 角 优 势 概念 
的 精心 运用 之 上 的 .这 一 方法 的 成 就 在 于 启示 了 其 他 研究 者 
去 发 展 将 多 变量 控制 问题 化 归 到 单 回 路 问题 ,例如 Mayne[16] 
所 作 的 序列 返 差 方 法 . 

在 将 不 相互 作用 或 部 分 相互 作用 方法 用 于 多 变量 控制 
时 ， 其 活力 在 于 在 设计 研究 的 最 后 阶段 能 分 开 使 用 经 典 的 单 
回路 频率 响应 方法 。 另 一 途径 则 是 把 传递 函数 矩阵 作为 一 个 
单一 对 象 来 研究 并 由 其 本 身 来 回答 : 如 何 把 经 典 单 回路 频率 
响应 方法 的 基本 概念 作 适 当 的 扩展 ? 在 多 变量 情况 下 什么 才 
是 极点 、 零 点 、Nyquist 图 以 及 根 轨迹 图 这 些 概念 的 合适 的 推 


，2 ， 


T? 对 于 这 类 问题 , 这 儿 的 工作 是 建议 并 且 曾 明 复 变 的 思想 
在 多 变量 反馈 系统 的 研究 中 可 以 担当 一 个 重要 的 角色 .将 
Nyquist 图 的 思想 扩展 到 多 变量 问题 的 早期 努力 是 由 Bohn[17， 
18] 作出 的 。MacFarlane [19] 给 出 Nyquist 稳定 性 判 据 的 一 
个 推广 ， 紧 接 这 一 启发 性 结果 ，Barman 455 Katzenelson [20] 
和 MacFarlane 与 Postlethwaite[21] 给 出 了 复 变 函数 的 基本 证 
FA. 随 着 多 变量 情形 Nyquist 稳定 性 判 据 的 推广 , 不 久 就 出 
现 了 根 轨迹 方法 的 补充 推广 [21, 22, 23, 24]. 

本 书 的 目的 是 将 Nyquist, Bode 与 Evans 方法 中 的 概念 
扩展 到 多 变量 系统 。 在 线性 反馈 设计 的 两 个 经 典 方 法 中 ， 
Nyguist-Bode 方法 是 研究 增益 为 频率 的 函数 ,而 Evans 方法 则 
研究 频率 为 增益 的 函数 。 第 三 章 剖 述 了 研究 复 增益 为 复 频 率 
的 函数 和 复 频 率 为 复 增益 的 函数 的 思想 是 怎样 扩展 至 多 变量 
情况 的 ， 这 里 用 到 与 传递 函数 矩阵 《其 行列 数 相同 ) 的 一 对 
解析 函数 : 特征 增益 函数 与 特征 频率 函数 ， 这 些 是 代数 函数 
[25] 且 每 个 函数 均 定义 在 一 合适 的 Riemann 曲面 上 [26]。 第 
二 章 给 出 了 一 些 基 本 知识 ,例如 被 考虑 的 多 变量 反馈 系统 的 
典型 表述 ;稳定 性 的 基本 定义 与 有 关 的 定理 ;基于 返 差 算 子 之 
上 的 开 环 与 闭环 特性 之 间 的 基本 关系 .第 三 章 也 包含 对 在 第 
四 章 中 出 现 的 多 变量 反馈 系统 的 广义 Nyquist 稳定 性 判 据 的 
背景 所 作 的 广泛 的 讨论 。 这 个 判 据 的 证 明 是 基于 运用 在 合适 
的 Riemann 曲面 上 定义 的 代数 函数 的 辐 角 原理 。 第 五 章 给 出 
一 个 多 变量 情况 的 逆 Nyquist 判 据 的 推广 ， 它 是 前 一 章 给 出 
的 广义 Nyquist 判 据 的 补充 .运用 在 第 三 章 中 发 展 了 的 材料 
在 第 六 章 把 根 轨迹 方法 扩展 到 了 多 变量 情形 ， 这 里 很 好 地 运 
用 了 代数 函数 理论 的 结果 .这 同样 阐明 了 代数 函数 的 基本 方 
法 可 以 怎样 用 来 寻求 多 变量 时 不 变 最 优 线性 调节 器 ， 当 二 次 
型 性 能 指标 输 和 人 项 的 权 趋 于 零 时 对 应 闭环 极点 的 渐 近 性 质 . 
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第 为 一 个 推广 ， 显 然 用 到 的 增益 变量 可 以 看 作 是 系统 的 


一 个 参数 ， 这 祥 已 发 展 了 的 方法 就 不 仅 可 运用 于 增益 与 频率 
而 且 也 可 运用 于 任何 参数 与 频率 ,在 第 七 章 用 引入 参数 根 轨 
迹 图 与 参数 Nyquist 轨迹 的 概念 ， 考 虑 了 多 变量 反馈 系统 中 
参数 变化 的 影响 .这 一 章 还 包含 了 一 些 为 进一步 研究 而 作 的 
试验 性 建议 与 设想 . 


次 要 的 信息 不 放 在 正文 而 留 在 附录 上 。 每 一 章 引 用 的 参 


考 文献 列 在 该 章 的 末尾 。 在 本 书 的 最 后 还 提供 了 一 个 总 文献 
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第 二 章 预备 知识 


本 书 对 由 几 个 多 输入 多 输出 的 子 系统 串联 而 后 形成 的 线 
性 时 不 变 反 馈 动 态 系统 ,来 考虑 经 典 的 Nyquist 与 Evans 方 
法 的 推广 。 在 这 一 考虑 下 , 本 章 给 出 了 多 变量 反馈 系统 的 描 
述 ,并 且 还 给 出 了 稳定 性 的 基本 定义 ,一 些 与 之 相关 的 定理 和 
基于 返 差 算 子 之 上 在 开 环 与 闭环 性 质 之 间 的 基本 关系 ， 


21 系统 描述 


线性 时 不 变动 态 系统 的 基本 描述 取 如 下 状态 空间 模式 : 
ål) = Ax(t) + Bule) 
(2.1.1) 
y(t) = Cx(#) + Dule) 
其 中 x(t) 是 状态 向 量 ，y(z) 为 输出 向 量 ，x(2) 为 输入 向 量 ; 
AC) 表示 Ce) 对 时 间 + 的 微 商 ; 4, B, C 与 忆 是 实 常数 和 矩 
BE, 为 方便 起 见 上 述 模式 在 意义 明确 的 前 提 下 简 记 为 SA, 
B,C,D) RS, HAA RA. 


图 1 状态 空间 模式 


一 般 S(4，B，C ,D) 可 认为 是 由 几 个 状态 空间 模式 的 子 
系统 
Si(A;, Bis Cis D;), 4C) = Aixi) + Biu) 


: (2.1.2) 
i= 1,2,-++,h, y(t) = Cat) + Diu) 
串联 而 成 的 ,如 图 2 所 示 。 
有 1 
usu | Yu 加 =y 
1 
HE eee oe RA i 
图 子 系统 的 串联 


例如 , 若 S 由 两 个 子 系统 S 与 9 组 成 , 则 S 用 方程 (2.1.1) 状 
态 空间 模式 表示 ,对 应 有 


0 一 | |， 两 子 系统 的 联合 状态 
a(t) 


u(t) = mle), Si 的 输入 (2.1.3) 
IU) = Ct), S: 的 输出 


A 0 B, 
aS EA ah R= lsn] 
C = (DC Cı), D = DD, 
而 由 几 个 子 系统 串联 的 状态 空间 模式 则 可 以 通过 反复 运用 上 
述 公 式 得 到 . 

由 于 状态 空间 模式 保留 了 系统 的 内 部 动态 结构 ， 因 而 可 
作为 一 种 内 描述 。 RADE (2.1.1) 两 边 作 单 边 的 Laplace 
变换 ,就 得 到 

s8(s) — x(0) = As(s) + BA(s) 
tc; = C&(s) + D&(s) 1A) 


e7 ， 


其 中 ACs) HAR x(t) 的 Laplace 变换 ,这 样 就 可 得 到 系统 的 外 
描述 或 输入 -输出 描述 .如 果 在 时 间 * = 0 时 的 初 条 件 全 为 
零 即 x(0) = 0, 则 输入 与 输出 变换 向 量 由 关系 式 
IC) = G(s)4(s) (2.1.5) 
相 联系 ,其 中 ` 
G(s) = C(sl,n — A)'B+D (2.1.6) 
In ym Bir ME C J| RREN. G(s) EHER 
s 的 有 理 函 数 矩 阵 ， 称 为 对 输入 -输出 变换 集 的 传递 函数 E 
PE, 或 开 环 增益 矩阵 . 传递 函数 矩 咱 可 用 来 描述 系统 在 具有 
指数 * 的 指数 函数 输入 下 的 响应 [2] ,因而 如 同 在 单 输入 单 输 
出 情况 下 一 样 , 复 变 量 s 可 认为 是 复数 的 频率 变量 ， 
Æ SCA, B, C, D) 如 图 2 所 示 ， 由 4 个 子 系统 《5S;(4;， 
Bj, Cis Di):i 二 1, 2,"*，h) 组 成 ,每 个 子 系统 的 传递 函数 
JERE 
Gis) = CiGsIn, — A;)'B; + D; (2.1.7) 
于 是 $ 的 输入 -输出 变换 向 量 由 关系 式 
Is) = G(s) Gyals)-+ GCs) a(s) (2.1.8) 
相 联系 ,因此 开 环 增益 矩阵 具有 明显 的 关系 式 
G(s) = G(s)Gs-s(s) GCs) (2.1.9) 
为 了 可 以 将 输出 由 反馈 连接 到 输入 ， 以 下 设 G(s) 是 一 个 也 
Bak. 


22 反馈 结构 


考虑 如 图 3 所 示 的 一 般 反 馈 结构 ， 这 里 反馈 系统 的 输出 
是 第 4 个子 系统 的 输出 ， 但 在 实际 上 它 也 可 以 是 较 前 一 个 子 
系统 的 输出 , 此 时 其 后 的 子 系统 可 视 为 反馈 补偿 器 。 参数 不 
是 整个 回路 总 的 实 增益 控制 变量 系统 的 输 人 与 输出 用 方程 


e) = rle) — 9@) (2.2.1) 
u(z) = ke) 
和 参考 输入 rO 联系 在 一 起 。 将 它们 和 方程 (2.1.1) 结合 一 
起 考虑 ,就 可 得 下 述 闭 环 系统 的 状态 空间 方程 : 
àlt) = Axle) + Ber(e) (2.2.2) 
y(i = Caxt) + Dre) 


其 中 
4 一 4 一 BCOCUn 十 D)C 
B, = kB — kB(k™I + DD 
C. = Un + kD) 'C 
D, = (klm + D)'D 


find fs, Fo ic 


r(t) et 


图 3 反馈 结构 


23 稳定 性 


稳定 性 是 反馈 系统 最 重要 的 一 个 要 求 。 对 于 一 般 时 变 非 
线性 系统 ,稳定 性 问题 是 很 复杂 的 ,但 对 于 线性 时 不 变 系统 来 
说 要 比 一 般 情形 下 简单 ,这 是 由 于 

O 所 有 稳定 性 性 质 关 于 时 间 是 恒定 的 , 且 

Git) 由 于 对 系统 的 状态 来 说 任何 解 均 正 比 于 在 时 间 零 时 
的 状态 (见方 程 (2.3.2), 于 是 所 有 稳定 性 性 质 就 都 是 全 局 的 . 

在 文献 中 有 不 少 稳定 性 的 定义 ， 大 致 说 来 可 分 为 两 类 . 
第 一 类 定义 是 关于 无 输入 的 自由 系统 的 稳定 性 ; 第 二 类 定义 
则 是 关于 具有 输入 的 受 迫 系统 的 。 两 类 稳定 性 均 在 下 面 讨 


. Qe 


沦 , 其 定义 以 及 与 之 相关 的 完美 的 定理 可 在 Willems 所 车 [31 
中 找到 . 


2.3-1 自由 系统 


对 由 方程 (2.2:2) 描述 的 具有 rl) 一 0 与 C. 一 7 (如 图 
3 所 示 ) 的 闭环 动态 系统 ,其 稳定 性 问题 可 归结 为 考虑 自由 系 
统 


2(t) = A.x(2) (2.3.1) 
方程 231) 的 平衡 态 显 然 是 原点 〈 设 4。 非 奇异 )， 因 而 式 
《2.3.1) 的 任 一 解 若 在 某 个 时 间 经 过 原点 则 在 以 后 的 一 切 时 间 
均 留 在 原点 , 这 样 的 一 个 解 称 为 零 解 。 原点 平衡 位 置 的 稳定 
性 用 下 述 定义 刻 划 ， 
定义 1 自由 系统 (2.3.1) 称 为 稳定 系 指 当 系统 从 原点 被 
扰动 后 ,其 一 切 后 继 运动 总 保持 在 原点 的 对 应 小 邻 域 内 。 
定义 2 ”自由 系统 (2.3.1) 称 为 浙 近 稳定 系 指 当 系统 从 原 
点 受到 小 的 扰动 后 ,其 一 切 后 继 运 动 将 返 向 原点 . 
定义 3 自由 系统 (2.3.1) 称 为 大 范围 或 全 局 渐 近 稳定 ， 
系 指 它 是 稳定 的 且 . 当 + 一 co 时 每 个 运动 均 收敛 至 原点 ， 
方程 (2.3.1) 的 通 解 是 
Ct x(to), t) = expl A(t — to) I] x(4)) (2.3.2) 
它 清楚 地 表明 若 自 由 系统 渐 近 稳定 则 它 也 大 范围 渐 近 稳定 。 
Æ JA 4 的 Jordan 正则 型 且 使 
A, = TJT” (2.3.3) 


J 
h, 
Ji 


其 中 每 个 Jordan 块 具有 形式 


。10 。 


具有 


1,32 A. 的 一 个 特征 值 , 于 是 能 够 证 明 有 [3] 
expl4.(¢ — 10)] 
一 Texp[J( — IT (2.3.4) 
具有 
exp[J(t — 10)] 
exp [fils — 10)] 
exp[Ja(t — to)] 


exp [Jalt — 4)] 


而 
exp [JC — to)] 
i p f/2y ses #1)! 
='° ne se A) exp[4,(2 — 10)] 
Lo 0 ` 0 : 1 


其 中 4 二 (+ — to), r 是 Jordan 块 J; 的 阶 数 。 这 样 自由 系统 
的 通 解 就 可 以 表示 为 
x(t;3 x(t)» to) 
= Texp[J(t — 1)]T x(t) (2.3.5) 
由 此 可 有 下 述 定理 ,其 证 明 见 [3].。 
定理 1 系统 (2.3.1) 的 零 解 是 渐 近 稳定 的 , 当 且 仅 当 和 矩 
阵 A. 的 全 部 特征 值 均 具 有 负 的 实 部 。 
定理 2 系统 (2.3.1) 的 零 解 是 稳定 的 ， 当 且 仅 当 Ai 
有 正 实 部 的 特征 值 ， 若 有 零 实 部 的 特征 值 则 它 只 对 应 1 阶 的 
Jordan $k, 


。11 


2.3-2 ZARA 


考虑 闭环 动态 系统 (2.2.2), 它 的 通 解 是 *[3] 
x(13 x(10), to) = exp[Alz — to)]x(10) 
+ ep [A — 1) Ber(r)dr (2.3.6) 
为 研究 该 系统 的 稳定 性 ,我 们 需要 引入 输入 -输出 稳定 性 的 概 
B: 
定义 4 一 动态 系统 称 为 输入 -输出 稳定 , 系 指 对 任何 有 
界 输入 ,不论 初 始 状态 如 何 均 产生 有 界 输出 。 
按 定理 1, 由 未 受 迫 系统 (2.3.1) 渐 近 稳 定 推 知 A 的 全 部 
特征 值 均 具 负 实 部 ,于 是 存在 正 实数 与 a, 使 
llexp (A <P exp(—ar), Vi>0 (2.3.7) 
其 中 || + || 玫 示 矩阵 或 向 量 的 Euclidean 范 数 [3]。 由 方程 
(2.2.2),(2.3.6) 与 (2.3.7) 则 有 
lly) < lex + Dr 
< IDO ll + Ceexp[ Ae 一 4)x(zo)l 


C. f exp[4.(¢ — r)]Ber(r)dr | 
cbMP 


dM 十 cpllz(ro)l + © Wh 


RH b= |B, c= Cl, d = iD, Ir@l<M, Y> 
o**, 


*) 原 书 误 为 
x(t; £to) to) = exp(Act)*(to) + flexptace 一 TD)]Bcr(r)dr. 
一 一 详 者 注 
**) BRS OISID) + Cepat) 
+e f lexpt AU 一 1) 1) Br lar <dM + cP x(t.) +E 


— RE 


< 


这 一 结果 可 归结 为 
定理 3 若 未 受 迫 系统 (2.3.1) 的 零 解 渐 近 稳定 , 则 受 迫 
系统 (2.3.2) 输 入 -输出 稳定 。 
注意 要 由 输入 -输出 稳定 性 推 知 原点 平衡 位 置 的 渐 近 稳 
定 还 必须 有 系统 (2.2.2) 状态 可 控 且 状态 可 观测 ;或 其 一 切 不 
可 观测 与 不 可 控 模 均 具 负 实 部 。 在 本 书 以 后 的 讨论 中 , 系统 
稳定 性 就 理解 为 输入 -输出 稳定 和 在 原点 平衡 位 置 的 渐 近 稳 
E. 
定理 3 重要 之 处 在 于 线性 时 不 变 系统 的 稳定 性 能 仅 由 系 
统 的 “4” 和 矩阵 的 知识 判定 。 稳定 性 所 关系 到 的 闭环 动态 系 
统 (2.2.2) 的 特征 值 是 指 满足 方程 
det [A], — 4.] 一 0 (2.3.8) 
的 4 的 值 ， 方程 (2.3.8) 的 左 端 称 为 闭环 特征 多 项 式 , 简 记 为 
CLCP(2), Rp 
CLCP (2) Set (41n — A.) (2.3.9) 
相应 地 对 开 环 系统 5S(4，B,C，D), 开 环 特 征 多 项 式 OLCP(4) 
定义 为 
OLCP(4)Adet (17, — A) 
一 det(L1。 一 Aidet(Aln, 一 4 
X det (115, — 4a) (2.3.10) 
FASE FES HB ES DK FES FR 
系 起 来 15]。 


2.4 一 般 反馈 结构 的 开 环 与 
闭环 特征 多 项 式 间 的 关系 


设 一 般 闭 环 结构 的 总 反馈 被 切断 且 子 系统 均 用 传递 函数 
HERR WEA 4 所 示 。 对 这 一 切断 点 说 来 对 应 的 返 差 矩 阵 


Boe 


是 
F()A1, + LCs) (2.4.1) 
其 中 
LCs) = &Gs(s) ++ G(s) = kG(s) (2.4.2) 
称 为 系统 的 返 比 矩 阵 [5], 返 差 算 子 是 由 进入 信号 的 变换 所 
产生 的 进入 与 返回 信号 变换 的 差 。 由 于 建立 反馈 线 的 本 质 在 
于 造成 两 信号 集 恒 等 从 而 使 在 其 中 的 差 恒 为 零 ， 这 样 在 反馈 
理论 中 返 差 算 子 就 起 着 重要 的 作用 。 Fs) 与 LG) 均 为 复 变 
Ks 的 有 理 函 数 和 矩阵 ， 本 书 的 关键 概念 都 是 围绕 这 些 矩 阵 的 
性 质 建 立 起 来 的 ， 返 差 矩阵 在 建立 开 环 与 闭环 特征 多 项 式 之 
间 的 关系 上 的 重要 性 受到 了 重视 ， 现 在 用 一 般 反 馈 结构 来 推 
SE. 


fis) &(s) G(s) ĝis) 


(a) 


(b) 
图 4 KRAH. (G) AR (b) FR 


若 取 方程 (2.4.1) 的 行列 式 且 将 G(s) 用 状态 空间 模式 表 
出 , 则 有 
det F(s) = det [Up + KCCsls 一 4)-1B + kD] (2.4.3) 
机 应 用 分 块 行列 式 计算 的 Schur 公式 则 可 将 上 式 改 写 为 
sla — A B 
一 AC Ig t sl 
+ det(sI, — 4) (2.4.4) 


det F(s) = det [ 


它 等 价 于 
Ini —BUn + kD)" 
det F(s) = det [Ea -| 


x eae B 
—kC lat kD. 
sls— AFBI t DOC| 0 | 


| + at Gta = A) 


-a| 


+ det(sI, — A) 
a det sla — A + BRM + D) Cdan + kD) 
det(sIn — A) 


(2.4.5) 
而 由 方程 (2.2.2) 有 
Ae = A — B(k`™im + D)’C 
因而 在 假定 det F(co) < 0 后 由 式 (2.4.5) 可 有 下 述 关 系 式 : 
det F(s) _ det [sln — Ae] 


det F(co) det [s1, — A] 


= det [s], — Ac] 
det[51 一 Aa] + edet [sln,— Ad 
CLCP(s) 
We 
™ OLCP(s) (24%) 


开 环 与 闭环 特征 多 项 式 OLCP(s) 与 CLCP(s) 的 零点 就 是 
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开 环 和 闭环 的 极点 或 其 对 应 的 特征 频率 . 


关系 式 (2.4.6) 指 出 了 和 矩阵 值 有 理 传递 函数 是 怎样 与 动态 


反馈 系统 的 稳定 性 紧密 相关 的 。 对 这 些 和 矩阵 及 其 特征 值 的 研 
究 开 膀 了 一 个 将 经 典 Nyquist [7] 与 Evans [8,9] 技术 进行 
适当 扩展 的 途径 .这 一 结果 将 在 第 三 章 中 给 出 . 
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第 三 章 ”特征 增益 函数 与 
特征 频率 函数 


在 线性 单 回路 反馈 系统 的 分 析 与 设计 中 ， 两 个 经 典 方法 
均 用 复 变 函数 去 研究 作为 频率 函数 的 开 环 增益 (Nyquist-Bode 
TE) 和 作为 增益 函数 的 闵 环 频率 (Evans 方法 ). 本 章 的 目的 
在 于 表明 这 些 方 法 是 怎样 用 与 复 变量 适当 的 有 理 隙 数 和 矩阵 相 
联系 的 特征 增益 函数 与 特征 频率 函数 扩展 至 多 变量 的 情形 . 


3.1 开 环 增益 与 闭环 频率 闻 的 对 侦 性 


对 图 4 所 示 的 一 般 反馈 结构 ， 由 第 2.4 节 可 以 有 如 下 的 
基本 关系 : 


det F(s) det (sIn — Ac) 
det F(0o) det (slp — A) (3.1.1) 
其 中 返 差 矩阵 FO) 给 定 为 
FCs) = In + G(s) (3.1.2) 


若 将 F(s) 代入 式 (3.1.1), 则 有 
det (sIn ~ Ac) _ det Um + &G(s)) 
det (sl, — A) det (Im + kD) 
— det (RI + G(s)) 
det (Alm + D) he 


若 再 应 用 表达 式 ， 


来 蔡 换 增益 ,其 中 8 可 以 是 复数 , 即 8e C (SPH). WE 
det (sIn — A) _ det (gl m — G(s)) (3.1.5) 
det (sla — A) det (g/, — D) 
而 由 式 (2.2.2) 给 出 的 闭环 系统 和 矩阵 A. 是 
A, = A — BCR" g + DC (3.1.6) 
用 式 (3.1.4) 代 换 k WA 
Ai=A+B(gln — D)-'C 
Ss(g) (3.1.7) 
这 样 表 达 式 (3.1.5) 可 改写 成 


det [+10 — S(g)] det[g1n — G(s)] (3.1.8) 
det [s], — 4] det [g/, — D] 


或 
det [s]a — S(g)] _ det g/m — G(s)] 
det [s]a 一 S(00)] det[gl1n— G(œ)] 

这 个 关系 式 表明 在 复 频 变量 :与 复 增益 变量 8 之 间 通 过 
其 对 应 的 “ 亲 ” 和 矩阵 SCg) 与 G(s) 表现 出 一 种 严格 的 对 偶 性 . 
这 个 在 频率 和 增益 的 法 则 之 间 的 对 偶 性 构成 了 将 经 典 复 变 方 
法 推广 至 多 变量 情形 的 基础 。S(g) 称 为 闭环 频率 矩阵 ; 它 的 
特征 值 是 闭环 特征 频率 且 显 然 依 赖 于 增益 变量 g。 开 环 增益 
矩阵 G(s) 的 特征 值 称 为 开 环 特征 增益 且 显 然 依 赖 于 频率 变 
Rs, EZEBE G(s) 与 S(g) 的 状态 空间 结构 

G(s) = CsI, 一 4)-B +D (3.1.9) 
S(g) = B(glm —D)C + A (3.1.10) 
以 后 ,它们 之 间 的 类 似 性 就 更 明显 了 . 

从 图 5 上 可 以 看 出 , 图 4(a) 在 具有 零 参 考 输 入 下 的 反馈 
结构 ,对 G(s) 的 状态 空间 描述 和 以 式 (3.1.4) 代 换 被 重 画 
出 来 的 闭环 特征 频率 变量 ;与 开 环 特 征 增益 变量 8 之 间 的 对 
BHE. 
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关系 式 (3.1.8) 的 重要 性 在 于 说 明 对 sc(4) 和 gg&afD) 
(等 价 于 假定 de FCco ) 关 0) 的 变量 说 来 ,有 
det [s1, — S(g)] 
= 0< 一 >det [g1,,— G(s)] 
=0 (3.1.11) 
其 中 o( A) 表示 4 的 谱 . 
这 一 点 告诉 我 们 关于 开 环 
特征 增益 作为 频率 的 函数 
的 知识 等 价 于 作为 增益 的 
函数 的 闭环 特征 频率 的 知 
识 , 由 此 可 以 断言 从 GCs) 
的 特征 增益 谱 线 来 断定 反 
馈 系统 的 稳定 性 应 该 是 可 
能 的 。 注 意 到 这 一 点 由 式 (3.1.8) 可 以 有 
_ det[g1», — G(s) JOLCP(s 
CLOP(s) TARET (o) 
这 一 表达 式 将 使 Nyquist 稳定 性 定理 推广 至 作为 频率 函数 
G(s) 的 特征 增益 轨 线 上 的 这 件 事变 得 很 明显 . 


图 5 说 明 : 与 8 之 间 对 偶 性 的 反馈 结构 


(3.1.12) 


32 代数 函数 : 特征 增益 函数 与 特征 频率 函数 


. G(s) 与 S(g) 的 特征 方程 
f Alg, s)Adet (glm — G(s)) = 0 (3.2.1) 
5 
Vs, 8g) 人 det (sla — S(g)) = 6 (3.2.2) 
是 复 变量 * 与 & 的 代数 方程 ， 每 个 方程 均 可 作为 g (或 *) 以 
s (R g) 的 有 理 函 数 为 系数 的 多 项 式 来 考虑 。 若 在 有 理 函 数 
域 上 为 不 可 约 , 则 每 个 方程 定义 一 个 代数 函数 对 [1; 附录 1]: 
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O 特征 增益 函数 g(s)， 它 给 出 开 环 特征 增益 为 频率 的 

Gi) 特征 频率 函数 s(g)， 它 给 出 闭环 特征 频率 为 增益 

一 般 说 来 方程 (3.2.1) 与 (3.2.2) 并 不 是 不 可 约 的 ， 且 每 个 
方程 将 定义 一 特征 增益 与 特征 频率 函数 的 集合 ， 为 了 简单 起 
见 同时 从 实际 提出 的 G(s) 与 SCg) 的 通常 情形 出 发 ， 可 以 认 
为 方程 (3.2.1) 与 (3.2.2) 在 有 理 函 数 域 上 不 可 约 是 合适 的 。 

虽然 方程 (3.2.1) 5 (3.2.2) 定义 的 是 同一 函数 g(*) 与 
SCG), 但 方程 (3.2.2) 一 般 讲 将 含有 较 多 系统 的 信息 ， 在 某 种 
情况 下 将 可 能 出 现 VG, g) 含有 与 & 无 关 的 因子 而 该 因子 在 
Alg, s) 中 并 不 出 现 。 这 类 因子 常 发生 于 下 述 情形 : 

G) 考虑 其 输入 为 第 1 子 系 统 的 输入 ,输出 为 第 4 子 系 
统 的 输出 所 对 应 的 开 环 系统 S(4，, B, C, D), 其 4 和 矩阵 具有 ` 
对 应 系统 不 可 观测 或 不 可 控 模 的 特征 值 ， 注 意 若 对 较 前 子 系 
统 的 输出 测量 是 有 效 的 , 则 实际 上 SCA, B,C ,DD) 的 某 些 不 
可 观测 模 可 以 是 可 观测 的 . 

Gi) 当 开 环 增益 矩阵 GO) 的 极点 与 零点 不 同 于 特征 增 
益 函 数 gC) 的 极点 与 零点 时 ;参见 第 3.3-3 节 . 

使 式 (3.2.1) 与 (3.2.2) 不 同 的 这 两 种 情况 ,清楚 地 表明 了 
把 Nyquist 类 判 据 发 展 到 GCs) 的 特征 增益 轨 线 上 所 出 现 的 
问题 ， 但 是 利用 g(s) 的 极点 和 零点 与 CCO) 的 极点 和 零点 之 
间 的 联系 ,再 仔细 地 考虑 不 可 观测 与 不 可 控 模 ,这 些 问题 就 都 
能 克服 广义 的 Nyquist 判 据 将 在 第 四 章 给 出 。 

在 下 一 节 , 将 给 出 特征 增益 函数 的 详尽 讨论 ,并 以 此 作出 
根 轨 迹 图 的 一 个 推广 。 而 第 3.4 节 将 对 特征 频率 函数 的 相应 
研究 给 出 一 个 广义 Nyquist 图 。 
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33 ”特征 增益 函数 


定义 特征 增益 函数 的 一 个 自然 途径 是 运用 特征 方程 
A(g, s)Adet [glm — G(s)] = 0 (3.3.1) 
一 般 讲 A(g, *) 可 化 成 下 述 形式 : 

Alg, s) 一 Al(g，5)A(8g， 5) Arle, 5) (3.3.2) 
HAT IA(g, Dsi=1,---,/} ReWHSMAHAE s 
的 有 理 函 数 域 上 不 可 约 。 设 不 可 约 因子 Ales) 具有 如 下 形 
A: 

Alg, s) = gii + a Cs)gi™ + ee + ails) = 0 (3.3.3) 
Hp n ER i 个 不 可 约 因 子 的 次 数 , 而 系数 {a;(s)37=1,--, 
1,7=1,2, +++, u) Æ s DARRA. Æ b) 是 系数 
dals) j= 1,2, ttes u) 的 最 小 公分 母 ， 则 方程 (3.3.3) 可 
写成 


bi (s)gii H bi (s)gi + oe 
+ bs) = 0; i= 1, 1 (3.3.4) 

其 中 系数 {6,053 i= 1,2, 067,05 7 1,2, 0015 8} 均 为 
s 的 多 项 式 。 由 式 (3.3.4) 定 义 的 复 变 函 数 gC) 称 为 代数 函 
数 [1; 附录 1]。 于 是 与 开 环 增益 甜 阵 G(s) 相 联 系 的 是 代数 
函数 集合 (el); i 二 1, 2。…，!}， 它 直接 与 GC) 的 特征 
值 相 联系 ，G(*) 的 特征 增益 函数 就 用 代数 函数 {g:(s); i = 
1，2，……，/!]} 的 集合 来 定义 。 

寻求 {Ai(g, 人 ); i 二 1, 2，*…*,!} 的 不 可 约 因子 并 由 此 
而 定义 特征 增益 函数 的 问题 与 Cs) 的 某 种 合适 的 正则 型 紧 
密 相关 ， 若 A(g, s) 可 以 化 简 成 8 的 线性 因子 , 则 GCs) 可 以 
化 成 Jordan 型 [3], 一 般 讲 这 并 不 常见 ,合适 的 正则 型 如 下 : 

令 
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0 0 +++ 0 一 ons) 
1 0 s+ 0 Sanal) 


CCA) 全 |0 1 0 ays) 1,4 >1 (3.3.5) 
0 0 ve 1 —a;(s) 
与 
C(A)B—a;,(s), 4=l (3.3.6) 
则 存在 变换 矩阵 ECs) 使 
G(s) = E(s)O(s)E(s) (3.3.7) 


其 中 OCs) 是 唯一 的 对 角 块 矩阵 ， 它 称 为 G(s) 的 不 可 约 有 理 
正则 型 并 由 下 式 

Q(s)Adiag[C(A) CCA), tt CCA)] (3.3.8) 
给 出 。 很 清楚 在 给 出 了 8(s) 以 后 不 可 约 因子 Alg, 5) RE 
易 得 到 ,而 对 任意 给 出 的 G(s), 如 何 寻 求 OCs) 的 办 法 将 在 附 
录 2 中 给 出 。 


3.3-1 特征 增益 函数 的 零 , 极 点 


研究 定义 特征 增益 函数 o(s) 的 方程 
Cg, 5) Bb(s)gi + blg H ++ 
+ 4,(s) =0 (3.3.9) 
以 下 将 设 
b(s)=0 与 bs) 0 
这 是 由 于 若 这 两 个 多 项 式 系数 之 一 成 为 零 ， 我 们 就 可 得 到 一 
个 首 项 与 零 次 寡 系 数 均 不 为 零 的 简化 形式 ， 于 是 这 一 降 了 阶 
的 方程 又 可 定义 一 个 代数 函数 ， 而 对 这 个 方程 说 来 前 述 对 系 
数 的 假定 则 是 满足 的 . 
可 能 发 生 包 与 刀具 有 公 因 子 的 情形 ,这 样 在 其 个 * 的 集 
合 上 这 二 者 将 同时 消失 ， 在 考查 这 个 影响 之 前 , 先 设 如 54, 
ae 


不 具有 公 因 子 ， 当 


bls) =0 (3.3.10) 
时 代数 函数 将 显然 为 零 , 而 当 
bs) > 0 (3.3.11) 
时 代数 函数 将 变 至 无 穷 。 由 此 就 定义 了 满足 式 (3.3.10) 的 s 
值 为 代数 函数 eCe) 的 零点 ,而 将 满足 
b(s) =0 (3.3.12) 


的 * 的 值 定义 为 代数 函数 gC) 的 极点 。 除 特别 声明 ， 术 语 
“零点 、 极 点 ”一 般 均 指 为 有 限 的 ， 而 对 ， 一 co 将 另 给 特别 注 
意 , 这 将 在 本 小 节 结 束 时 再 谈 。 

为 了 能 使 方程 (3.3.10) 与 (3.3.12) 在 bo(*) 与 bx(*) 具 有 公 
因子 时 同样 可 以 用 来 定义 gC) 的 零点 与 极点 ,我 们 就 必须 说 
明 其 合理 性 。 首 先 来 处 理 一 个 显然 的 情况 , 即 当 所 有 G) 
§=0,1, +++, } 含有 一 个 公 因子 时 ， 这 一 公 因子 可 以 从 方 
程 中 简单 地 约 去 而 留 下 一 个 为 一 适当 代数 函数 所 满足 的 方 
程 。 然 而 若 OCs) 与 6.(s) 有 这 种 公 因 子 而 系数 的 一 个 非 空 集 
{bw《s)， bails), oot, DC) 没有 这 样 的 公 因子 ， 那 么 对 式 
(3.3.9) 左 端 以 bo(s) 除 之 , 则 有 

g + 2 git tee + PhS) pe 十 eee 
bols) 


bs) 
PAS) meg weg OD 一 
+ BY) get +3) 0 (3.3.13) 


若是 bs) 5 和 (5) 的 公 因子 的 零点 , 则 当 s 一 了 时, 系数 集 


Bis) ,bs) lo = 
(22, GO) arte mn te BIR =E HAE eC) 


的 一 个 极点 。 
又 设 名 与 汪 具 有 公 因 子 而 系数 的 某 个 非 空 集 {2(5)，…， 
bm(s)} 不 具有 这 种 因子 , 则 当 * 一 了 时 代数 方程 (3.3.9) 就 可 
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用 下 式 代 替 : 
53) gS) Hee + bl) gi O 一 0 (3.3.14) 
其 中 是 上 述 公 因子 的 零点 , Hbl) = 0, ++, bn) 关 0. 
因此 就 必 有 
2() 一 0 
这 表明 :确实 是 代数 耳 数 的 一 个 零点 。 
由 此 我 们 可 以 断言 方程 (3.3.10) 与 (3.3.12) 是 可 以 用 来 定 
义 代 数 函数 cs) 的 有 限 零 点 与 有 限 极 点 的 . 并 且 这 一 定义 
也 能 使 我 们 处 理 零 点 与 极点 相 重 合 的 情形 。g(s) 的 极点 与 堆 
点 多 项 式 分 别 以 PC) 与 l) 表示 ,并 定义 为 
Pl) Abs) 
za(s) 人 bi(s) 
其 中 如 与 如 是 由 久 与 5 用 各 自首 项 系数 除 后 所 得 的 首 一 多 
FR. 
KABE s 一 co 处 o(s) 的 目的 出 发 , 设 
s= (3.3.16) 


(3.3.15) 


则 有 
Olg, s) = O(g, 2) =z Wg, z) (3.3.17) 
其 中 9 是 g(s) APRA. 在 > = 0 的 任 一 邻 域内 (点 = 
= 0 不 包括 在 内 ) 方 程 @(g,*) = 0 等 价 于 方程 (8,z) 一 0， 
因而 若 考 虑 方程 
yg, xz) 会 co(z)8 + cls)g t+ ++ 
+c(2)=0 (3.3.18) 
则 可 推 知 
Gi) s = co 是 特征 增益 函数 g(*) 的 极点 , 当 且 仅 当 ceo(0) 
=05. 
Gi) s= 00 是 特征 增益 函数 g(s) HER, SAMS 
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DL 


(0) = 0, 

Wik BLE EA ENTE — BA LAST A a E 
阵 说 来 , 按 定义 (参见 第 2.1 节 ) os) 将 不 可 能 有 无 穷 远 的 极 
点 ,事实 上 容易 看 出 * 一 co 的 特征 增益 函数 g(s) 的 值 就 是 刀 
的 特征 值 . 


3.3-2 ”传递 函数 矩阵 极点 和 零点 的 代数 定义 
设 TO) 为 一 复 变量 :的 m X 1 有 理 函数 矩阵 ， 则 有 一 
TCs) 的 正则 型 ( 即 Smith-McMillan 型 ) M (s), 使 
TCG) = H(s)M(s)J(s) (3.3.19) 
其 中 m x mie HC) 与 1 x 1 和 矩阵 J(s) 均 为 单 模 态 矩阵 
《其 行列 式 为 与 无关 的 非 零 常数 )。 若 TORRA r ( 即 对 
几乎 所 有 的 s,，7T(s) 的 秩 为 +), WMG) 具有 如 下 形式 : 


men | Or we | 


1 (1 REEE tor (3:320) 
其 中 
#(s) = diag LECO, E25) ... 8s) 3. 
MPa) = dag |B» Be oe Hey CSD 
并 且 


Gi) EA Els) 可 整除 一 切 ECs); 
Gi) 每 个 OAs) 可 整除 一 切 $i-;(s)。 
用 适当 的 H(s), MGs), JG) 的 分 块 写法 ,可 有 
TO = HH Ao [VO ES 
= H,(s)M*(s)Ji(s) (3.3.22) 
其 中 M*(s) 满足 方程 3.3.21)。 由 此 Ts) 可 表示 为 
TCG) = HAs) [ cise (#9)] Gs) 


= > h(s) HOL Fs) (3.3.23) 


其 中 


G) hl); i= 1, 2, +++, r} ERE AAs) 的 列 ; 
GD (FC) i= 1, 2, +++, r} 是 矩阵 LG) 的 行 。 
我 们 知道 
r <min (J, m) 
AAG), 7(s) 是 对 一 切 :， 均 具 满 秩 m 与 ! 的 单 模 态 矩阵 ， 如 
RTO) 是 对 具 输 入 变换 向 量 as) 与 输出 变换 向 量 I) 的 
系统 所 对 应 的 传递 函数 矩阵 , 则 任何 输入 向 量 Cs) 均 可 


Heys Ones OAG) (3.3.24) 
转换 成 一 输出 向 量 9(:)。 对 单 输入 单 输出 情形 有 
9(s) = 48D) ae) 


Ws) 
其 中 不 为 常数 ,传递 函数 为 
s) = ke) 
g(s) K 


其 零点 是 使 e(s) 为 零 的 * 值 而 其 极点 是 使 J《s) 为 零 的 s 
fi. 由 此 对 于 非 零 的 ACs) 说 来 当 * 是 g(s) 的 零点 时 对 应 
IO 的 模 为 零 ， 而 当 * 是 g(s) 的 极点 时 对 应 90s) 的 模 为 任 
RK. 因此 刻 划 Ts) 的 零点 与 极点 的 一 个 自然 途径 是 采用 
s 的 那些 值 ,对 这 些 值 , TESA AO AT AR MAY IC 为 
零 ， 或 者 对 有 限 的 ACOA HD 为 无 穷 大 , 其 中 I 
表示 标准 向 量 范 数 ， 这 种 标量 情形 下 的 自然 扩展 可 直接 引出 
TC) 的 极点 与 零点 的 定义 , 即 用 Smith-McMillan 型 中 项 
els) 

oils) , 


所 作 的 定义 ,这 也 是 下 述 两 结果 成 立 的 原因 . 
零点 引 理 对 lOl +08 s APR, lOl OS, 4H 
仅 当 某 个 els) = 0, 

极点 引 理 对 lA < 0, [19s] >, 4AM 4H 
个 os) > 0, 

由 这 些 考虑 自然 地 可 以 引出 下 述 定义 [3]. 

T(s) 的 极点 Ts) 的 极点 系 多 项 式 集合 {4(s); i 二 1， 
2,7, r} 的 全 部 零点 的 集合 。 以 下 记 TO) WRAY IA, 
Pas teea Paha & 

Pr(s) = (s — Ps — Pa)t + -Cs — Pa) (3.3.25) 
并 称 它 为 T(*) 的 极点 多 项 式 , 且 有 


PCs) = I gi(s) (3.3.26) 


Tis) 的 零点 ”多项式 fe, (si = 1,2, +++, 7} 集 的 全 
部 零点 均 称 为 T(s) 的 零点 ， 以 下 用 {21t zo} 表示 TCs) 
的 零点 , 令 


ar(s) = (sO—2)°"*(s — za) (3.3.27) 
并 称 其 为 T(*) 的 零点 多 项 式 , 且 有 
zr1(s) 一 II Cs) (3.3.28) 
重要 的 是 应 指出 zz(s) 与 pr(s) 未 必 一 定 互 质 ,因而 简单 地 用 
detT(s) 的 分 母 与 分 子 的 零点 来 定义 TC) 的 极点 与 零点 是 不 
恰当 的 . 
计算 极点 多 项 式 与 零点 多 项 式 的 法 则 
应 用 Smith-McMillan 型 的 办 法 来 确定 TO) 的 极点 与 零 
点 并 非 总 是 方便 的 ， 特 别 在 手 算 时 . 下 述 法 则 能 给 出 与 用 
Smith-McMillan 型 定义 同样 的 结果 。 
极点 多 项 式 法 则 PC) 是 TO) 的 所 有 阶 的 所 有 非 零 
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子 式 的 最 小 首 一 公分 母 . 

零点 多 项 式 法 则 2) ETO RRA ri O È TC) 
WK) 子 式 的 分 子 的 最 大 首 一 公 因 式 。 这 些 子 式 全 都 正好 有 
prts) 作为 公分 母 ， 


3.3-3 ” 开 环 增益 矩阵 G(s) 代数 定义 的 极点 -零点 与 对 应 特 
征 增益 函数 集合 的 极点 -零点 间 的 关系 


在 广义 Nyquist 稳定 性 判 据 的 确立 过 程 中 ， 作 为 至 关 重 
要 的 一 步 是 建立 代数 意义 上 的 极点 ,零点 与 复 变 理论 的 联系 ， 
即 建立 与 对 应 特征 增益 函数 集 的 极点 零点 的 联系 。 
在 展 式 
det [g] m — G(s)] 
= g” tase” +++ + an(s) (3.3.29) 
中 ,系数 a(s) 由 熟知 的 公式 
det [gln — G(s)] = g” — [tr G(s)]g”" 
十 [如 G(s) 的 2 阶 主子 式 ] "° -e 
十 (一 1)"det G(s) (3.3.30) 
推 知 应 均 为 G(s) 的 子 式 的 适当 的 和 ， 因 此 极点 多 项 式 bo(;s) 
是 首 一 多 项 式 , 它 可 以 由 G(s) 的 所 有 阶 全 部 非 零 主 子 式 的 最 
小 公分 母 得 到 。 
现在 一 方 阵 G(s) 的 极点 多 项 式 es) 是 首 一 多 项 式 , 它 
可 以 由 G(s) 的 所 有 阶 的 所 有 非 零 子 式 的 最 小 公分 母 得 到 . 
因此 , A ecl) 是 由 所 有 非 零 非 主子 式 的 最 小 公分 母 得 到 的 
首 一 多 项 式 且 已 消去 全 部 与 b) 公共 的 因子 , 则 可 有 
Pels) = ec(s)bo(s) (3.3.31) 
进而 由 于 


det G(s) = an(s) = ey (3.3.32) 
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又 由 G(s) 的 Smith-McMillan 型 有 


det G(s) =a aa (3.3.33) 
Hp a 为 与 * 无 关 的 标量 ,于 是 我 们 必然 得 到 
ze(s) = ec(s)bm(s) (3.3.34) 


在 不 少 场合 G(s) 的 非 零 非 主子 式 的 最 小 公分 母 可 除 尽 
bo(s)， 此 时 ecl) 将 为 1, 因此 CO) 的 极点 与 零点 多 项 式 相 
应 地 就 是 bo(s) 与 bw(*)。 一 般 讲 一 个 复 变量 方 阵 函 数 G(s) 
将 会 有 一 个 由 工 个 不 可 约 且 满足 方程 (3.3.3) 的 特征 增益 函数 
形成 的 集合 ， 因 而 其 极点 与 零点 多 项 式 的 一 般 形 式 可 以 写成 


Pals) = ec(s) I H09 (3.3.35) 


aos) = eols) JẸ Pils) (3.3.36) 


其 中 第 j 个 特征 增益 函数 的 极点 与 零点 多 项 式 对 应 为 bals) 
Sbu (s). 


阐明 极点 -零点 关系 的 例子 


令 

Giese =、 
(s + 1)(s + 2)(s — 1) 

td 0 ] 
—(s+1)(s +2) (s — 1)(s + 1). 
ane Se, 0 

E s+1 

o} =1 1 
s—1 s+2 


` +296 


对 G(s) 说 来 极点 多 项 式 显然 是 
Pls) = (s #1) +2)s — 1) 


而 对 应 零点 多 项 式 是 
as) —(s—1) 
GC) 的 特征 方程 是 
det [g7 — G(s)] (e-—4)le aa 

因而 不 可 约 特征 方程 是 

Als) =e M 
与 

As, 一 一 
它们 可 改写 成 


(二 lg 一 1=0 
5 
(十 2)g 一 1 一 0 
因而 对 特征 增益 函数 aG) 与 (5) 说 来 极点 与 零点 多 项 式 就 
是 
Pol) = bus) = (6+ 19 sn) = 8D =1 


Pos) = bas) = G H2), a) = HD 1 


现在 由 G(s) 的 非 零 非 主子 式 的 最 小 公分 母 得 到 的 首 一 多 项 

式 中 消去 其 与 56s)( 二 bi(s)bw(s)) 的 全 部 公 因子 , 则 有 
ec() 一 了 一 1 

它 满足 关系 


Pels) = ec(s) Il LAO) 


+306 


2 
zc(s) = eds) Il biG) 


3.3-4 ”特征 增益 函数 的 Riemann 曲面 


特征 增益 函数 g(s) 是 用 一 个 不 可 约 方程 
ZO 十 hs)e + +b(s)=0 (3.3.37) 

定义 的 ， 一 般 讲 它 有 : 个 不 同 的 有 限 根 ， 例 外 的 情况 仅 当 : 

G) bo(s) 一 0， 由 于 方程 的 次 数 降低 , 当 bC) 一 0 时 有 
一 个 或 几 个 根 将 趋 于 无 穷 大 ， 或 

Gi) 方程 具有 重 根 . 

这 后 一 种 情形 可 在 有 限 的 :时 出 现 ， 当 且 仅 当 判 别 式 的 
表达 式 为 零 ， 判 别 式 是 方程 系数 的 整 有 理 函 数 , 记 为 Dels), 
在 附录 3 上 对 其 进行 了 讨论 ， 


特征 增益 函数 的 正常 点 

特征 增益 函数 的 正常 点 [1;6] 系 指 复 平 面 上 任 一 使 bls) 
= 05 DCs) = 0 的 有 限 点 。 
特征 增益 函数 的 临界 点 


eC 的 临界 点 [1;61] 系 指 复 平面 上 任 一 使 
bs) =0 RÈ DC) =0 
成 立 的 点 以 及 点 * 一 00, 


特征 函数 的 分 支点 


D,(s) = 0 
的 解 称 为 特征 增益 水 数 的 有 限 分 支点 ， 无穷 远 点 成 为 分 支点 
是 指 式 (3.3.18) 的 判别 式 DAD) = 0. 
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在 每 个 正常 点 由 于 判别 式 非 零 ， 确 定 特征 增益 函数 的 方 
程 (3.3.7) 将 具 :个 不 同 的 根 。 代 数 函 数理 论 指 出 在 复 平面 
上 单 连通 区 域内 除了 临界 点 以 外 , 特征 增益 函数 g(*) 的 取 值 
是 一 解析 函数 集合 .其 中 的 每 一 个 函数 均 称 为 特征 增益 函数 
&(s) 的 分 支 ， 基 于 标准 的 解析 开拓 方法 联系 到 代数 方程 的 性 
质 , 可 以 指出 这 些 不 同 的 分 支 刚好 组 成 代数 函数 gC) 的 一 个 
整体 ， 这 可 归结 在 关于 代数 函数 的 这 一 基本 定理 之 中 : 一 个 
不 可 约 代数 方程 (3.3.37) 准确 地 在 带 孔 的 平面 上 确定 一 个 : 
值 函数 [7]. 

用 这 样 的 方法 定义 的 函数 称 为 代数 函数 ， 这 可 以 作为 熟 
知 的 复 变 初等 函数 的 一 种 自然 推广 。 一 个 复 变 的 初等 函数 有 
其 定义 域 与 值 域 两 复数 集合 C ， 一 个 代数 函数 其 值 域 为 一 复 
数 集合 C ,但 其 定义 域 RR 则 是 Riemann 曲面 ， 由 于 代数 函数 
的 Riemann 曲面 在 这 一 工作 中 起 着 关键 的 作用 ,因此 重要 的 
是 要 在 其 定义 与 公式 之 外 对 其 思想 有 个 直观 的 理解 

设 对 代数 函数 的 一 分 支 有 一 宕 级 数 表达 式 ， 这 洋 一 个 玫 
达 式 常 称 为 一 个 函数 单元 我 们 可 以 想象 把 这 个 笑 级 数 的 收 
倒 圆 从 纸 上 割 下 来 ， 纸 盘 上 的 个 体 点 构成 了 该 单元 唯一 函数 
值 的 承受 域 ， 现 将 这 个 初始 的 函数 单元 用 第 二 个 震级 数 作 解 
析 开 拓 ， 于 是 可 以 想象 另 一 个 收敛 圆 被 制 了 下 来 且 部 分 地 与 
第 一 个 粘 在 一 起 ,如 图 6 所 示 . 这 个 粘 成 一 块 的 部 分 又 构成 了 
唯一 函数 值 的 承受 域 且 作 为 一 个 值 单 叶 的 区 域 ， 若 再 进一步 
作 和 解析 开拓 , 则 又 可 将 一 个 新 的 盘子 粘 到 原 有 的 上 面 去 ， 重 
复 作 这 样 的 解析 开拓 ， 就 会 出 现 一 个 盘子 搭 上 了 另 一 个 但 并 
不 是 用 解析 开拓 直接 联系 的 那 一 个 , 如 图 7 所 示 。 RE 
的 盘子 与 另 一 个 交 迭 的 盘子 粘 在 一 起 ， 当 且 仅 当 它 们 是 同一 
函数 值 的 承受 域 ; 但 是 , 若 它 们 是 不 同 函数 值 的 承受 域 , 则 就 
认为 它们 之 闻 没 有 粘 在 一 起 ， 而 认为 它们 是 两 片 承 过 着 不 同 
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图 6 解析 开拓 图 7 重复 解析 开拓 

的 函数 值 ,而 且 在 复 平面 的 这 一 部 分 是 重 失 地 放 着 . 
将 上 面 这 一 过 程 继 续 做 下 去 ,一 个 想象 的 覆盖 : 页 的 
曲面 就 形成 了 ,其 中 + 是 代数 函数 的 次 数 ， 在 Riemann 曲面 
的 形成 中 ,这 些 页 可 用 各 种 方法 互相 连通 ,这 也 包括 这 样 两 页 
的 连接 , 即 它们 彼此 被 很 多 页 分 开 着 ， 虽然 这 样 一 种 结构 并 
不 能 在 三 维 空间 中 进行 下 去 ， 但 给 出 一 个 所 要 求 过 程 的 合适 
的 拓 朴 描述 是 并 不 困难 的 。 这样 一 个 想象 的 曲面 就 称 为 多 值 
代数 函数 的 Riemann 曲面 ,在 Riemann 曲面 上 代数 函数 值 
的 整体 完全 是 按 单 值 意义 铺展 开 的 。 在 包括 复 平面 在 内 的 这 
+ 页 的 每 一 页 上 ,每 一 点 是 也 仅 是 一 个 函数 值 的 承受 点 ， 

在 附录 4 上 给 出 了 一 个 构造 Riemann 曲面 的 方法 , 这 包 
括 在 复 平面 上 运用 割 口 的 方法 ， 它 对 于 说 明 问 题 可 能 有 帮 
助 ， 令 代数 函数 eG) 有 7 个 临界 点 ta， a2，…… ,4,}, 设 用 线 
工 将 它们 互相 连接 然后 再 连 至 无 穷 远 点 。 这 样 我 们 在 割 过 的 
HPE C-L EELT + 个 不 同 的 解析 函数 (als), (5)，…， 
本 (9)}， 这 些 函 数 中 每 一 个 均 可 用 标准 的 解析 开拓 办 法 穿 过 
割 口 。 由 解析 开拓 的 基本 原理 可 知 , 若 一 解析 函数 在 其 定义 
域 的 一 部 分 满足 一 代数 方程 ， 则 它 必然 在 其 解析 开拓 的 每 一 
区 域内 满足 这 一 方程 。 因 此 就 一 定 有 
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O RA: 个 不 同 的 解析 函数 在 割 过 的 平面 C- 世 上 满足 
确定 的 代数 方程 . 

Gi) 这些 解析 函数 (B): i= 1, 2，…, 直 的 任 一 个 
经 解析 开拓 给 出 一 个 解析 函数 , 它 同 样 满足 确定 的 代数 方程 
由 此 可 知 联系 割 口 L 一 边 的 解析 函数 的 集合 必然 是 联系 市 
口 另 一 边 的 解析 函数 集 的 一 个 置换 。 因 此 在 割 口 也 两 边 的 解 
析 函 数 上 《用 其 计算 值 集 合 ) 利用 适当 的 比较 方法 ,就 可 以 产 
生出 一 个 合适 的 定义 域 ， 在 其 上 可 以 确定 一 个 单 值 的 解析 函 
数 , 它 将 这 个 定义 域 连续 单 值 地 映射 到 复 平面 上 。 显然 这 个 
函数 就 是 代数 函数 它 已 被 看 作 是 一 个 单一 整体 ， 而 上 述 构 成 
的 定义 域 就 是 Riemann 曲面 . 

对 于 读者 说 来 ， 知 道 对 任 一 给 定 的 代数 函数 均 能 将 其 
Riemann 曲面 构造 出 来 ， 在 其 上 函数 的 值 形成 一 位 置 的 单 值 
函数 , 这 对 于 了 解 本 书 已 是 充分 的 了 .运用 Riemann 曲面 的 
概念 解析 函数 理论 中 不 少 标准 的 关系 与 性 质 可 以 推广 至 代数 
函数 的 情形 , 特别 是 在 Riemann 曲面 上 辐 角 原 理 也 成 立 ， 附 
SK 5 给 出 了 辐 角 原 理 的 推广 。 

作为 特征 增益 函数 定义 域 的 Riemann 曲面 称 为 频率 曲 
面 或 s 曲面 。 当 开 环 增 益 和 矩阵 G(s) 是 m x m, 且 其 对 应 的 
特征 方程 不 可 约 ( 通 常 实际 中 常见 的 情形 ) 时 o M3 H da m 
页 频率 平面 或 :平面 所 组 成 . 


3.3-5 ”广义 根 轨迹 图 


特征 增益 函数 eC) 是 复 变 函 数 , 其 极点 与 零点 均 位 于 频 
率 曲面 定义 域 上 。 在 频率 曲面 上 用 画 出 g《(:) 的 常 位 相 与 党 
辐 值 的 曲线 来 展示 其 性 质 是 方便 的 。 若 将 附录 4 中 勾 划 的 计 
算 方 法 用 来 构 壮 曲 面 ， 则 画 出 常 位 相 与 常 轻 值 的 曲线 就 是 一 
个 简单 的 过 程 。 频 率 曲 面 能 够 想象 成 是 具有 复 增 益 参 数 & 的 
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所 有 可 能 值 的 所 有 可 能 的 闭环 特征 频率 的 集合 ， 一 旦 曲面 已 
用 常 位 相 与 常 辐 值 曲线 勾 刘 好 ， 我 们 就 有 了 在 闭环 特征 频率 
与 开 环 增益 闻 的 一 个 直接 对 应 而 由 于 曲面 是 由 zr 页 复 频 率 
平面 构成 ,于 是 对 每 个 * 值 将 有 和 个 对 应 的 特征 增益 . 

由 方程 (3.1.4) 有 


g(s) = ar (3.3.38) 


所 以 随 着 实 控制 变量 《的 变化 闭环 极点 (特征 频率 ) 勾 划 出 
一 个 图 形 , 该 图 形 等 价 于 g(s) 的 180° 等 位 相 曲 线 ”。 方程 
(3.3.38) 是 对 确定 单 回路 根 轨迹 方程 的 直接 推广 .g(s) 的 180° 
位 相 图 是 一 个 多 变量 根 图 ， 即 闭环 极点 随 增 益 控 制 变量 变 
化 的 图 。 多 变量 根 图 画 在 Riemann 曲面 上 这 一 事实 与 单 输入 
单 输出 的 情况 比较 表明 其 复杂 化 了 的 特性 ， 而 后 者 的 根 图 是 
画 在 简单 的 复 平面 上 的 ( 即 单 页 的 Riemann 曲面 )， 多 变量 和 根 
图 有 时 也 称 特征 频率 图 或 多 变 景 根 轨 迹 . 

在 第 3.2 节 就 已 指出 一 般 对 G(s) 与 5(g) 而 言 特征 方程 
将 有 所 不 司 , 对 5(g) 的 方程 来 讲 可 能 包含 某 些 与 8 无 关 的 s 
的 因子 ， 因 而 这 些 因子 对 应 闭环 极点 也 独立 于 8 或 等 价 地 独 
LFR 从 根 图 的 观点 这 些 是 孤立 点 对 应 退化 根 图 ， 退 化 根 
图 不 能 通过 用 频率 曲面 上 gs) 的 180° 等 位 相 图 求 得 ， 实际 
上 特征 频率 轨 线 是 作为 一 个 轨 线 集合 在 单 页 复 平 面 上 产生 出 
来 的 , 即 用 & 在 增益 平面 沿 负 实 轴 变 化 时 对 应 5S(g) 的 特征 值 
描 出 , 这 一 方法 将 自动 挑 出 退化 的 轨迹 。 与 经 典 Evans 根 图 
方法 的 共同 点 是 特征 频率 轨 线 常用 增益 控制 变量 《二 一 8 
KEB. 


© 这 儿 指 K 为 正 实 数 .一 一 译 者 注 
。35 。 


33-6 ”频率 曲面 与 特征 频率 轨 线 的 例子 
作为 一 个 例子 ,考虑 图 3 所 示 的 一 般 多 变量 反省 结构 ,其 
开 环 传递 函数 为 (增益 矩阵 
G(s) 一 


1 s—1 s 

(1.25) (s 十 1)G + 2) (ve Res 

EME 2 阶 因而 合适 的 曲面 可 视 为 由 两 页 复 * POR. A 
8,9 上 画 出 了 在 两 页 上 gC) 取 常 位 相 与 常 辐 值 的 曲线 ， 可 办 
识 不 连续 性 的 割 口 在 图 上 用 黑 粗 线 画 出 ,而 特征 频率 轨 线 , 唱 
是 gC) 的 180° 等 位 相 的 曲线 ， 它 是 用 方块 符号 标 出 的 。 特 
征 频 率 轨迹 指出 增益 控制 参数 从 零 开 始 上 升 引起 系统 从 经 验 
稳定 、 不 稳定 到 又 重新 稳定 。 这 一 现象 与 在 右 半 平 面 存在 一 
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图 9 第 2 页 的 频率 曲面 
个 支点 (在 * +) 的 情况 十 分 吻合 . 


注意 ,由 于 我 们 用 开 环 增益 矩阵 来 刻 划 整个 反馈 结构 , 因 
而 就 不 考虑 任何 不 可 观测 与 不 可 控 的 成 分 。 


3.4 ”特征 频率 函数 


确定 特征 频率 函数 (e) 的 自然 途径 是 利用 特征 方程 
Vs, g)Adet [s1 — S(g)] =0 > (3.4.1) 
它 是 一 个 代数 函数 ， 在 前 一 节 详 细 讨 论 特征 增益 函数 的 办 法 
可 以 直接 应 用 于 此 而 只 需 将 * Se 互 换 一 下 . 
作为 特征 频率 函数 的 定义 域 的 Riemann 曲面 称 为 增益 曲 
面 或 8 曲面 , 它 由 2# 页 复 增益 平面 或 8 平面 组 成 ,这 是 由 于 对 
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每 个 给 定 的 8 值 闭环 特征 频率 有 ”个 值 (闭环 极点 )。 增益 
曲面 可 以 看 成 是 与 所 有 闭环 特征 频率 相 联 系 区 开 环 增益 矩阵 


G(s) 的 所 有 可 能 的 开 环 特征 增益 。 用 类 似 于 增益 函数 


g(s) 


的 做 法 合适 的 是 在 增益 平面 上 ,用 在 曲面 上 画 出 s(g) 的 常 位 
相 与 常 辐 值 的 曲线 的 办 法 来 展现 :(g) 在 增益 曲面 上 的 性 质 
Ael) 一 样 8) 也 具有 极点 与 零点 但 其 作用 却 是 很 不 同 的 . 


3.4-1 广义 Nyquist 

增益 曲面 的 每 一 页 上 均 用 s《g) 的 常 位 相 与 常 辐 值 
画 好 后 ， 就 可 分 成 对 应 左 半 平面 与 右 半 平面 闭环 特征 频 
区 域 ， 因而 给 定 这 样 一 个 描画 好 的 曲面 后 , 一 眼 就 可 看 
《或 等 价 地 A) 的 那些 值 对 应 稳定 的 闭环 极点 。 很 清楚 稳 


RAY 
率 的 
出 8 


定 与 
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不 稳定 区 域 间 的 边界 就 是 *(8g) 的 +90 ERR. (e) 的 
土 90? 等 位 相 线 是 单 环 Nyquist 图 的 自然 推广 , 称 为 特征 增益 
Re. 

实际 上 特征 增益 轨 线 是 在 复 增益 平面 上 用 GG) 的 特 


图 11 第 2 页 的 增益 曲面 


Q 右 半 平 面 区域 


OC) 太 半 平面 压 域 


revere 特征 增 答 轨 线 
coc 分 发 点 间 的 害 口 


图 12 图 10 的 略图 ;以 强调 右 半 与 左 半 平 面 区 域 
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533 19-2 


图 13 图 11 的 赂 图 以 强调 右 半 与 左 半 平 面 区 域 
征 值 当 * 在 ， 平面 上 走 过 所 谓 Nyquist D 回 线 时 描 出 的 . 
如 果 我 们 考虑 虚 轴 的 一 部 分 ， 则 我 们 可 以 算出 对 应 特征 什 
AGO), ++, gm(jw) 的 轨 线 的 集合 (在 本 书 中 i = V 一 D ,其 
算法 如 下 : 


ice s 


y 


(a) 


图 14 第 1 页 的 增益 曲面 . 


G) 选择 一 第 频率 的 值 ,例如 ov; 

Gi) 计算 复数 矩阵 CG); 

Gii) 运用 标准 的 算法 计算 G(jw,) 的 特征 值 ， 它 是 一 复 
数 集合 ,表示 成 {zo} 

Gv) 在 复 平面 上 标 出 这 些 数 { 有 (jos)}; 

O) 对 角 频 率 的 其 它 的 值 ©, we，…' 重 复 前 述 过程 ， 然 
后 用 基于 特征 函数 各 分 支 的 连续 性 进行 整理 的 办 法 将 这 些 标 
出 的 点 连 成 曲线 . 

对 于 第 四 章 中 发 展 的 广义 Nyquist 稳定 性 判 据说 来 ， 
Nyquist D 回 线 是 用 标准 的 顺 时 针 方向 扫 过 的 . 


《a)“ 完 全 "页 ;《(b) 一 0.2 附近 的 小 区 域 
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34-2 ”增益 曲面 与 特征 增益 轨 线 的 例子 


作为 一 个 例子 考虑 第 3.3-5 节 中 讨论 过 的 开 环 增益 矩 
阵 , 它 具有 下 述 最 小 状态 空间 实现 : 
一 2 0 1 0.6 
4=[; iP z=; asl 
24 一 1.6 o o 
c= i ae DS h a] 
系统 具有 两 状态 因而 合适 的 增益 曲面 将 由 两 页 复 8 平面 构 
R. EA 10 与 图 11 上 分 别 画 出 了 这 两 页 (e) 的 等 位 相 与 等 
辐 值 曲线 。 特 征 增益 轨 线 是 *(8) 的 土 90° 等 位 相 曲 线 , 在 图 
上 用 带 有 “xX” 的 线 表示 . 
图 12 与 13 是 图 10 与 11 的 略图 ， 以 突出 说 明 闭 环 极点 
的 右 与 左 半 平 面 区 域 。 从 这 一 略图 上 易于 看 出 增益 控制 参数 
处 的 稳定 性 界限 。 由 于 g =—k", 当 丰 由 零 沿 正方 向 增加 ， 
则 对 应 8 的 值 在 每 页 上 就 由 一 oo 沿 实 轴 向 原点 移动 。 在 第 
1 页 当 1.25<k < 2.5 时 对 应 的 8 在 右 半 平面 区 域 ,而 在 第 2 
页 ,对 正 值 的 《对 应 的 & 不 会 进入 右 半 平面 区 域 , Alt, 闭环 
系统 在 Kk < 1.25 与 2.5<k < co 时 都 是 稳定 的 。 
若 锡 由 零 沿 负 向 增加 ， 则 8 的 临界 值 在 每 一 页 上 均 由 co 
沿 实 轴 向 原点 移动 .在 这 两 页 上 当 一 00 k<— 1.875 时 8 的 
值 位 于 右 半 平面 。 因 此 对 应 正 反 馈 的 《的 负 值 当 一 1.875<<% 
< O 时 对 应 闭环 系统 将 是 稳定 的 . 
车 将 有 标 度 的 增益 曲面 投影 至 复 增益 平面 C 上 , 就 得 到 
特征 增益 轨 线 正规 的 一 种 描述 ， 加 上 表示 右 半 与 左 半 平面 两 
这 域 的 曲线 图 ， 稳 定性 就 可 以 通过 考虑 右 半 平 面 区 域 与 单个 


临界 点 (二 十 10 ) 间 的 关系 预测 出 来 。 但 由 于 图 形 重 选 的 
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表示 方法 难于 理解 。 因 此 虽然 在 广义 Nyquist 稳定 性 判 据 ( 第 
四 章 ) 中 计算 绕 数 如 同 Saeks [10] 指出 那样 对 系 绕 稳定 性 并 不 
是 基本 的 , 但 它 确实 还 是 提供 了 一 个 在 增益 平面 C 上 预测 闭 
环 系 统 稳定 性 的 最 简单 的 方法 。 “ 

由 增益 曲面 上 画 好 的 图 可 以 确定 闲 环 极点 ， 因 而 去 确定 
闭环 系统 的 相对 稳定 性 是 可 能 的 。 现 在 通过 对 具有 单位 的 
前 述 例子 利用 和 寻找 优势 闲 环 极点 的 办 法 来 叙述 它 。 为 此 具有 
如 图 14 与 图 15 所 示 的 具 常 实 部 值 与 常 虚 部 值 曲线 的 增益 曲 


面 是 合适 的 。 
令 优势 闭环 极点 是 
实 回 线 标志 : 
-09 -tO -H 2 -13-14 -1 is 
< BARRE: os 

特征 增益 轨 线 woroor 

07 10 Ho es 

706, 

-os 

-04 Os. 

+3. 

-02~ 


0 
os 
ov 
or -10 
o 
7 


8 
O9 io m i2 a isis o- 


图 15 第 2 页 的 增益 曲面 
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sg = atjp 
则 “是 通过 临界 的 , 一 1 与 & 的 值 中 任 一 个 的 最 小 (在 模 上 》 
负 实 部 回 线 , 而 8 是 对 应 的 虚 部 回 线 。 由 图 14 与 15 有 
cx 全 一 0.05 与 8 一 0 


由 此 
sa = —0.05 
用 手 算 的 优势 闭环 极点 是 
sa = —0.0528 


由 增益 曲面 同样 可 以 确定 特征 频率 轨 线 。 除 可 能 的 单 点 
轨 线 外 , 根 轨 线 是 在 特征 增益 轨 线 具有 位 相 180? 处 的 :的 值 . 
因此 由 增益 曲面 上 画 好 的 图 ， 根 轨 线 就 由 常数 值 回 线 在 每 一 
页 上 穿 过 负 实 轴 时 的 值 确 定 。 类 似 地 ， 从 一 画 好 图 的 频率 曲 
面 上 ， 由 常数 值 曲线 穿 过 虚 轴 时 的 值 确定 特征 增益 轨 线 也 是 
可 能 的 。 
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BOR 广义 Nyquist 稳定 性 判 据 


在 线性 反馈 系统 理论 中 Nyquist 稳定 性 判 所 [1] 是 一 个 最 
基本 的 结果 , 它 在 多 变量 情况 下 的 推广 具有 很 大 的 意义 ,这 一 
推广 是 由 MacFarlane[2] 作 为 对 反馈 系统 分 析 与 设计 的 特征 轨 
线 法 的 一 部 分 给 出 的 [3] ， 而 当时 并 无 合适 的 证 明 予 以 支持 . 
广义 Nyquist 稳定 性 定理 的 证 明 是 由 Barman 55 Katzenelson 
给 出 的 [4,5] ,他 们 的 方法 忽视 了 所 含量 的 某 些 关键 的 代数 性 
质 , 且 过 分 地 依赖 于 复 平面 上 的 割 口 ,使 得 这 种 提 法 在 技术 上 
复杂 但 却 又 掩盖 了 其 结果 在 本 质 上 的 简单 性 ， 这 一 章 的 目的 
在 于 给 出 一 个 一 般 反 馈 结构 下 广义 Nyquist 式 稳定 性 判 据 的 
严谨 证 明 , 其 基础 在 于 复 变 理论 的 基本 结果 :， 运 用 于 定义 在 
合适 的 Riemann 曲面 上 的 代数 函数 的 辐 角 原理 [附录 51. 

本 章 提出 并 证 明了 两 个 Nyquist 式 稳定 性 的 判定 ， 每 一 
个 的 适用 性 依赖 于 子 系统 的 表征 ， 判 据 的 两 个 提 法 在 第 4.1 
节 中 给 出 而 证 明 在 4.2 节 中 给 出 . 
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若 图 3 所 示 的 一 般 反馈 结构 的 子 系统 均 用 状态 空间 表 
示 , 则 判 据 的 下 述 提 法 是 可 用 的 。 

提 法 1 一 般 反 馈 结构 是 闭环 稳定 , 当 且 仅 当 : 

O 特征 增益 轴线 集合 绕 临界 点 (一 上 + i0) 反 时 外 方 
向 图 数 的 净 和 等 于 GO) 在 右 半 平面 的 极点 数 ; 


Gi) 特征 增益 锁 线 不 经 过 临界 点 (一 二 + i0); 


GD 通 向 无 穷 远 的 特征 增益 训 线 的 分 支 妆 等 于 CO) 在 
He BLA BS 

(iv) 开 环 系统 SCA, B, C, D) (HERE 4 的 对 应 系统 不 
可 控 与 不 可 观测 模 部 分 的 特征 信 全 在 左 半 平面 ， 此 时 开 环 系 
统 系 指 输入 即 第 一 个 子 系统 的 输入 而 其 输出 即 第 4 个 于 系统 
的 给 出， 

REARS AREER, REMETTR 
统 均 不 存在 不 可 观测 与 不 可 控 模 所 对 应 包括 虚 轴 在 内 的 右 半 
平面 的 特征 值 , 则 可 采用 下 述 提 法 : 

提 法 2 一 般 反馈 结构 是 闭环 稳定 , 当 且 仅 当 : 

O 特征 增益 轴线 集合 绕 临 界 点 (一 起 + 50) 反 时 外 广 
向 图 数 的 净 和 等 于 GCG), GXs)，.…，G4(s) 在 右 半 平面 极 
点 数 之 和 ; 

(ii) 特征 增益 轴线 并 不 经 过 临界 点 (—F 十 10); 

(ii) 通 向 无 穷 远 的 特征 增益 轴线 的 分 支 数 等 于 G(s) 在 
iA AB 

ER AREO 或 (i) 不 成 立 ， 则 闭环 系统 有 一 个 或 


几 个 极点 在 虚 轴 上 ,因而 此 时 不 是 输入 -输出 稳定 但 平衡 位 置 
仍 可 能 是 稳定 的 . 
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在 第 2.4 节 已 经 指出 ， 图 4 所 示 切 断 点 对 应 的 返 差 算 子 
与 开 环 和 闭环 特征 多 项 式 之 间 的 关系 为 


det F(s) _ CLCP(s) 
det F(oo) OLCP(s) G2) 


+466 


这 个 基本 关系 式 是 基础 ， 广 义 Nyquist 稳定 性 判 据 的 证 明 就 
在 此 基础 上 给 出 . 
证 明 的 第 一 阶段 是 考虑 返 差 矩阵 FO) 的 特征 值 方程 ， 
它 是 
det [fm — F(s)] = 0 (4.2.2) 
如 同 G(s) 的 特征 方程 一 样 ,一 般 讲 它 可 表示 成 形式 为 
BoC s + dF!l+ + dis) = 0 
1 一 1，2，……，!/ (4.2.3) 
的 不 可 约 代数 方程 的 乘积 ， 由 这 些 方程 确定 出 代数 函数 
{fi(3):i = 1,2, +} RA. 因此 ， 如 同 在 第 3.3-3 节 那 
儿 关 于 GG) 的 极点 ,零点 与 特征 增益 函数 的 极点 、 零 点 间 的 
关系 所 作 的 讨论 一 样 , F(s) 的 极点 和 零点 多 项 式 与 代数 函数 
{ACi = 1,2, +91} 的 极点 和 零点 多 项 式 之 间 的 关系 可 表 
述 如 下 : 


1 
Pr(s) = ex(s) TI da(s) 
fei 
5 (4.2.4) 
1 
zp(s) = er(s) II dui(s) 
i=1 


由 定义 一 般 反馈 结构 下 开 环 特征 多 项 式 为 
OLCP (s)Adet (sla — A)- 

=det (sIn, — Ax) det (sIn, — Ar), © *det(sIn,— Aa) (4.2.5) 

容易 指出 [6] 有 
det [s15, — Ai] = pe(s)pals) (4.2.6) 

其 中 po,(s) 是 传递 函数 矩阵 GCO) 的 极点 多 项 式 而 首 一 多 项 
式 p(s) 的 零点 是 子 系统 的 那些 解 看 零点， 这 些 解 看 零点 联 
系 到 A, 对 应 不 可 控 与 不 可 观测 模 的 特征 频率 (特征 值 )， 开 
环 特征 多 项 式 可 以 写成 


OLCP(s) = Pe,(s)Po,(5)+ + *Po,(s)Pa,(s)Pa.Cs)* > Pag) 
(4.2.7) 
对 于 开 环 增益 矩阵 CO) 的 极点 多 项 式 pe(*)， 它 与 子 系统 传 
递 函数 的 极点 多 项 式 的 关系 为 
Po(s)Px(s) = Po,(s)Pa,(s)+ + *Po,(s) (4.2.8) 
其 中 户 (*) 的 零点 是 GCs), Gls), +--+, GG) 的 那些 极点 ， 
它们 在 形成 G(s) 时 就 消失 了 ，P:(s) 的 零点 实际 上 是 系统 
SCA, B, C, D) 不 可 控 不 可 观测 模 的 子 集合 ， 这 个 系统 的 输 
人 是 第 一 子 系统 的 输入 ， 而 其 输出 是 第 4 个 子 系统 的 输出 . 
FE SCA,B,C,D) 的 不 可 控 不 可 观测 模 的 完全 集合 是 pa(s) 
的 零点 集 , 它 是 


Pals) = Pr(s)Pals)* Ba,C) (4.2.9) 
于 是 开 环 特征 多 项 式 可 以 改写 或 
OLCP(s) = Po(s)Puls) (4.2.10) 


或 者 
OLCP(s) = pe)pe(s)pals) as) (4.2.11) 
若 将 式 (4.2.10) 与 基本 关系 式 (4. 2.1) 联 系 起 来 , 则 有 


a )=pls det F(s) 3 
; GLCP) Pal. Ors FC) (4.2.12) 
由 F(s) 的 Sn 标准 形 有 


m g 226). 
W BOs 6 pre) (4.2.13) 


其 中 8 是 与 :无 关 的 标量 ; n RO) 的 结构 与 方程 (3.1.2) 可 清 
楚 地 知道 首 一 多 项 式 ESS zi(s) 将 有 相同 次 数 , 从 而 有 


` det F(0o) = 8 (4.2.14) 
这 样 联合 方程 (4.2.12 一 I4) ,闭环 特 征 多 项 式 可 以 给 出 为 
CLCP(s) = pal(s)pe(:) re (4.2.15) 
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若 由 方程 (4.2.4) 与 (3.3.35) 对 上 式 进 行 替换 , 则 这 一 表达 式 可 
重 写成 


I dai(s) 
CLOP(s) = pals)ecls) II bols) =- 一 一 (4.2.16) 
i II ads) 


车 对 方程 (3.1 D 用 特征 值 移动 定理 , 则 可 得 FCs) 的 特征 函 
数 ae ):i = 1,2, +++, 7} 与 特征 增益 函数 {8;(s); f= 1, 
2, °° 0} Saas 

f(s) =1 +k g(s); i= 1,2,+2%,2 (42.17) 
Et I eh AGRA AL E 


Il da(s) = I AO) (4.2.18) 
于 是 由 式 (4.2.16) 则 可 得 到 


CLCP(s) = pi(s)ecls) Il LAO) (4.2.19) 


注意 多 项 式 I ig Cs) 是 依赖 于 增益 控制 变量 的 ， 


由 关系 式 (4.2. 19) 可 以 推 知 (参见 第 2.3 OR 闭环 系统 稳 
定 必 要 而 又 充分 的 条 件 是 
G) es) =0, RAAY FHP 


件 Gi), 而 这 同 经 典 Nyquist AY 
对 联系 到 返 差 算 子 FOIN 

有 一 对 应 的 Riemann 曲面 集合 ,在 其 上 合适 的 特征 代数 函数 

{f(s):i = 1,2, + +55} 就 是 单 值 的 ,由 这 些 曲面 映射 到 思平 
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HE 1 到 1 的 并 且 是 连续 的 。 如 果 考 虑 方程 (4.2.3) 定 义 的 集 
合 的 第 i 个 方程 ,方程 的 次 数 为 1;, 因 而 对 应 的 Rieman 曲面 
R, 是 由 4 页 复数 :平面 组 成 . 现 设 对 C 的 4; 页 在 它们 还 未 
粘 成 Rj 之 前 的 每 一 页 上 将 Nyquist D 回路 画 好 ”如 图 16 所 
示 , 这 样 在 曲面 粘 成 后 ，Nyquist 回路 就 并 合 组 成 绕 过 Ry, 的 
右 半 平面 的 闭 Jordan 回路 的 集合 ,这 样 就 可 以 将 扩展 了 的 加 
角 原 理 运用 到 R, 的 每 个 右 半 平 面 之 上 .因此 , 对 一 特定 的 、 
未 必 单 连通 的 、 但 边界 系 由 Nyquist D 曲线 组 成 的 右 半 平面 
区 域 ,代数 函数 f(s) 在 该 区 域内 的 零点 与 极点 数 之 差 就 等 于 
该 区 域 的 边界 曲线 在 映射 i(s) 下 ,其 象 在 C( 复 太平 面 ) 上 对 
原点 的 顺 时 针 绕 数 。 由 此 ， 若 考虑 R 上 的 全 部 右 半 平 面 区 
域 ,再 将 扩展 了 的 辐 角 原理 运用 于 这 些 区 域 , 则 有 

N(fis 0) = Zy, — Py, (4.2.20) 


其 中 
G) NC, 0) EE C 上 用 象 曲 线 集合 对 原点 的 顺 时 针 绕 
数 的 净 和 ,这 些 象 曲线 是 映射 js) 对 粘 合成 的 Rj, 上 的 曲线 
集 的 象 ,该 曲线 集 是 在 R, 未 粘 好 时 先 画 出 而 后 随 Rj 粘 成 而 
形成 的 。 
Gi) Zy 是 f(s) 在 右 半 平面 的 零点 数 . 
Gii) Py 是 f(s) 在 右 半 平面 的 极点 数 ， 
注意 ,在 D 回 线 上 的 由 缺 ， 仅 从 扩展 的 辐 角 原理 的 发 展 理 
论 的 观点 看 是 必要 的 ,而 在 实际 上 DD 回 线 可 以 直接 取 作 虚 轴 . 
闭环 稳定 性 的 条 件 ii) SNE 
(Gti = 1,2, 6+, 7} (4.2.21) 
在 右 半 平 面 或 在 虚 轴 :上 没有 零点 ,从 而 可 以 用 
Zy = j= 1,2,..,1 
或 (4.2.22) 
NC, 0) = —Py3 j = 1,2, -1 
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5 平面 


a RRA 


图 16 理论 的 Nyquist D 曲线 


加 上 代数 函数 式 (4.2.21) 在 虚 轴 上 没有 零点 的 条 件 。 由 条 件 
式 (4.2.227 可 推出 


i 了 
> N(f;, 0) = =>; P; (4.2.23) 
i=. =I 
同时 式 (4.2.22) 也 可 由 式 (4.2.23) 推出 ,由 此 闭环 系统 稳定 性 
的 必要 充分 条 件 可 以 改写 成 


G) ecl) 一 0， 没 有 右 半 平面 的 根 ; 
Gi’) ec(s) = 0， 在 虚 轴 上 没有 根 ; 


Li 1 
Gi) PNG, 0) = — >; Pys 
G) O i 一 1，2，…,!} 在 虚 轴 上 没有 零点 ; 
O) Pals) 仅 有 左 半 平 面 的 零点 . 


现在 ,正如 方程 (3.3.35) 与 (4.2.18) MERRI, ecl) 连同 
特征 水 数 集 的 极点 多 项 式 一 起 组 成 了 66) 的 极点 多 项 式 , 这 


. 5] 。 


使 得 我 们 去 考虑 组 合 条 件 (1 与 (证 ) 为 一 单个 的 特征 值 条 件 : 


i 
Dy N(f;, 0) = 一 Pc 
is 


(4.2.24) 


其 中 Pc 是 GC) 在 右 半 平 面 的 极点 数 。 对 于 闭环 系统 稳定 性 
ERT (这), Gv) 5 (v) 已 成 立 的 条 件 下 , 式 (4.2.24) 的 充 


分 必要 性 可 证 明 如 下 : 
由 方程 (3.3.35) 与 (4.2.18) 可 以 有 


Po=e+ > Py; 
其 中 。 为 ec(s) 在 右 半 平面 的 夫 点 数 ,同样 可 有 
INGO = > Z= > Py; 
联合 这 两 个 表达 式 , 则 有 
DNs0) = > Z; +e— Pe 


其 中 D Zp e RIP WHE RMR, 
为 确立 式 (4.2.24) 的 必要 性 , 设 
> N(fis 0) < —Pe 
因而 和 由 式 (4.2.27) 就 可 推 得 
3 Z = 0 Ble 0 


这 就 断定 了 条 件 (4.2.24) 对 闭环 稳定 是 必要 的 ， 
对 于 充分 性 , 设 


£ 
DNG, 0) = —Pe 
i=l 
则 由 方程 (4.2.27) 就 必 有 


。52 。 


(4.2.25) 


(4.2.26) 


(4.2.27) 


t 
Saats 


因而 闭环 系统 是 稳定 的 , 即 条 件 (4.2.24) 的 充分 性 成 立 . 

因此 现在 我 们 可 以 指出 对 于 质 环 系统 的 稳定 性 下 述 条 件 
是 充分 必要 的 : 

Li 

G) >) NG 0) 一 一 Pa; 

i= 

Gi”) AOS 1, 2, +++, 1) EMA BAS 

Git”) ec(s) 在 虚 轴 上 没有 零点 ; 

Gv’) pa(s) 仅 有 左 半 平面 零点 . 

由 方程 《4.2.17) BIA, Nyquist D 回 线 在 映射 f(s) 与 
Reis) F, C 上 的 象 曲线 集 间 以 单位 移动 互相 联系 , 于 是 稳定 
性 条 件 GO) 可 代替 为 

> N(kg;, —1) = —Pe (4.2.28) 
i=] 
其 中 N(hgi, 一 1) 是 在 C 上 用 倍 后 的 特征 增益 轨 线 绕 点 
(一 1 十 10) 的 顺 时 针 绕 数 的 净 和 ,如 同 在 经 典 Nyquist 判 据 


时 一 样 , 倍 可 以 用 计算 特征 增益 轨 线 绕 临 界 点 (—F+70) 
的 绕 数 加 以 哈 开 , 即 用 
上 
SNe 一 证) =—P, (4.2.29) 
来 代替 条 件 (4.2.28). 
由 方程 (4.2.17) 同 样 可 知 条 件 G) 等 价 于 
Bio) 送 一 二 一 2， sl (4.2.30) 
实际 上 这 个 对 应 特征 增益 轨 线 不 通过 临界 点 (一 二 + 70), 


条 件 Gil”) 也 能 用 一 个 比较 实际 的 条 件 来 代 蔡 . 
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出 方程 (3.3.35), G(s) 的 极点 多 项 式 给 定 为 
Pas) = ccls) [I bas) (4.2.31) 


其 中 bC) 是 第 i 个 特征 增益 函数 gi(s) 的 极点 多 项 式 . 
此 这 一 点 是 清楚 的 , 即 ecl) 在 虚 轴 上 没有 零点 当 且 仅 当 特 
征 增益 轨 线 在 无 穷 远 的 分 支 数 就 等 于 G(s) 在 虚 轴 上 的 极点 
RM. 注意 , 若 轨 线 的 一 个 分 支 走 向 无 穷 , 而 同时 有 另 一 分 支 由 
无 穷 远 返回 ,此 时 应 仔细 不 能 把 这 个 算 成 两 个 无 穷 分 支 。 

于 是 闭环 稳定 的 必要 充分 条 件 可 改写 成 


(iD Sn(e, 一 a) 一 一 Po 
Gi) 特征 增益 轨 线 不 通过 临界 点 (去 +0); 


(说 》 特 征 增 益 轨 线 通 向 无 穷 的 分 支 数 刚好 是 G(s) 在 虚 
轴 上 的 极点 数 ; 
Gv) es) 仅 有 左 半 平 面 的 零点 . 
这 就 完成 了 广义 Nyquist 稳定 性 判 据 提 法 1 的 证 明 . 
如 果子 系统 完全 由 其 传递 函数 矩阵 表述 ， 则 它们 的 状态 
空间 表示 每 一 个 子 系统 都 是 可 观测 与 可 控 的 ,于 是 
Pals) = Pas) = "= Pa C) IL (4.2.32) 


因而 有 

pals) = PCs) (4.2.33) 
于 是 Nyquist 判 据 提 法 1 的 条 件 (iv) 等 价 于 2.) 仅 有 左 
半 平 面 的 零点 ”如 果子 系统 用 状态 空间 表述 , 但 我 们 已 知 子 
系统 没有 不 可 观测 或 不 可 控 模 位 于 右 半 平 面 或 虚 轴 上 时 ， 上 
述 情况 也 能 发 生 . 在 这 种 情况 下 条 件 (i) 与 (iv) 可 合并 给 出 
一 个 不 需要 解 耦 零点 信息 的 判 据 . 

为 说 明 这 一 点 我 们 首先 证 明 , 当 
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G) pels) = PGs) (4.2.34) 
或 

Gi) Pals) = PCs) P2,(s)P2,(5) ++ PaCS) (4.2.35) 
其 中 {ps3):i = 1,2, ---, A} 的 零点 均 在 左 半 平面 . 
则 条 件 G), Gi’) 5 Ww) 等 价 于 


SNG.) = -5 Po; (4.2.36) 


其 中 Pe; 是 GG) 的 右 半 平面 极点 数 。 对 闭环 系统 稳定 说 来 
在 条 件 Gi’) 5 Gv) 已 满足 时 , 式 (4.2.36) 的 充分 必要 性 可 按 
如 下 证 明 : 

由 方程 (4.2.9),(4.2.18) 与 (4.2.31) 有 


> Pg =e + > Py, + P, (4.2.37) 
Heh P, 是 p(s) AEF MS AR, 将 它 与 式 (4.2.26) & 
并 就 有 
S40) = 5 Zy +e +P, — 5 Po; (4.2.38) 
同样 ,注意 7 = 


， Ps) 一 pa(s)bo(s) + Pas) 
所 以 若 用 Ps 表 pa《s) 在 右 半 平面 的 零点 数 , 则 


Pa = P, (4.2.39) 
且 方程 (4.2.38) 变 成 为 
SD NC,0) = > Zp +e +Pi— 5 Po; (4.2.40) 
为 确立 条 件 (4.2.36) 的 必要 性 , 设 = 
do0) x -$ Pa 
于 是 由 方程 (4.2.40) 推 知 


1 
DZ, 0 Re +0 MR Pr 0 
A 


或 上 述 的 任 一 组 合 ， 由 此 系统 将 是 闭环 不 稳定 的 。 从 而 证 明 
了 条 件 (4.2.36) 对 闭环 稳定 说 来 是 必要 的 . 
对 充分 性 设 
N50) = => Po; 
于 是 由 方程 (4.2.40) 则 必然 有 
了 
> 2, = ¢ = Ps =0 
i=1 
因而 系统 是 闭环 稳定 。 即 式 (4.2.36) 的 充分 性 得 以 确立 . 
因此 我 们 已 证 明了 当 条 件 (4.2.34) 或 条 件 (4.2.35) 成 立 
时 对 疡 环 稳定 性 说 来 下 述 条 件 是 充分 必要 的 : 
了 h 
G) DNG, 0) 一 一 > Pos 
i=l i=l 
GE) {EiS 1,2, A) 在 虚 轴 上 没有 零点 ; 
Citi”) ecl) EMR ERAS A, 
但 我 们 在 证 明 提 法 1 中 已 经 证 明 


ZNG) = P(e 一 证) 
条 件 Gi”) 等 价 于 特征 增益 轨 线 不 经 过 (一 证 十 10] 这 一 点 ; 


以 及 条 件 Gi) 等 价 于 特征 增益 轨 线 的 无 穷 分 支 数 正好 等 于 
G(s) 在 虚 轴 上 的 极点 数 。 由 此 就 完成 了 广义 Nyquist 稳定 性 
判 据 提 法 2 的 证 明 
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为 叙述 稳定 性 判 据 , 考 虑 在 第 3.3-5 节 与 3.4-2 节 已 用 过 
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PGF, BIS JL — A I RR Hes EH A a FE 2 ERI AD 


be s—1 s 
O= TG HDGFD| —6 aol 


由 定义 该 系统 没有 不 可 观测 或 不 可 控 的 模 ， 校 验 一 下 问题 的 
提 法 则 可 以 认为 系统 结构 是 由 一 个 子 系统 Gs) = OC) 所 
构成 ,因此 可 以 直接 应 用 稳定 性 淹 据 的 提 法 ! 或 所 法 2。 

GCs) 的 极点 多 项 式 是 

PCs) = (6 + 1)G +2) 

它 在 右 半 平面 或 虚 轴 上 没有 零点 ， 因 此 闭环 稳定 性 成 立 仅 需 
以 特征 增益 轨 线 对 临界 点 的 反 时 针 绕 数 为 零 、 特 征 增益 轨 线 
不 通过 临界 点 和 不 存在 无 穷 分 支 ， 特征 增益 轨 线 画 在 图 17 
上 ,于 是 可 得 到 下 述 稳 定性 条 件 : 

G) -o< 二 一 0.8， 特 征 增益 轨 线 不 绕 过 也 不 经 
过 临界 点 ,因而 闭环 系统 对 0 <k 二 1.25 是 稳定 的 ; 

Gi) aS 二 一 0.8, 特征 增益 轨 线 通过 临界 点 ， 因 而 闭 
环 系统 对 k= 1.25 是 不 稳定 的 ; 

Gii) —0.8 < — $ < 0.4, 临界 点 有 一 个 顺 时 针 绕 数 ， 
因而 闭环 系统 对 1.25 <k < 2.5 是 不 稳定 的 

Gv) = = 一 0.4, 特征 增益 轨 线 通过 临界 点 ,因而 闭环 
系统 对 = 2.5 是 不 稳定 的 ; 

O) -04< 二 二 0， 特征 增益 轨 线 不 绕 过 也 不 经 过 
IFA Aim AMAR 2.5 << k< co 是 稳定 的 ; 

(i) 一 rs = 0, 特征 增益 轨 线 通过 临界 点 , 因而 闭环 系 
BEN! k = oo 是 不 稳定 的 ; 
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图 17 特征 增益 轨 线 
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(ii) 0< = < 0.533, 特征 增益 轨 线 绕 过 临界 点 顺 时 


针 两 圈 , 因 而 闭环 系统 对 一 00 < k< 一 1.875 是 不 稳定 的 ; 


(viii) =F = 0.533, 特征 增益 轨 线 通过 临界 点 ， 因 而 闭 


MAAN k= 一 1.875 是 不 稳定 的 ; 


(ix) 0.533 < a 过 0, 特征 增益 轨 线 不 绕 过 也 不 通过 


临界 点 ,因而 闲 环 系统 对 一 1.875< <0 EREN. 


[6] 


[7] 


18] 
[9] 


注意 《 取 负 值 对 应 的 为 正 反 馈 . 
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SAR 广义 逆 Nyquist 稳定 性 判 据 


在 这 一 章 , 我 们 将 把 单 输入 单 输出 反馈 系统 的 逆 Nyquist 
判 据 推广 至 一 般 反 馈 结构 ， 并 以 此 作为 前 一 章 Nyquist 稳 
定性 判 据 的 补充 。 这 一 发 展 是 以 开 环 增益 矩阵 G(s) 与 逆 特 
征 增益 函数 集合 之 间 的 联系 为 基础 的 ， 逆 Nyquist 图 的 对 应 
集合 将 称 为 逆 特征 增益 轨 线 。 


5.1 逆 特 征 增益 函数 


Nyquist 稳定 性 判 据 推广 的 一 个 主要 方面 ,在 第 四 章 给 | 
的 是 代数 函数 集合 一 特征 增益 函数 与 传递 函数 矩阵 G(s) 
间 用 特征 方程 
A(g,s) = det[gln — G(s)] = 0 (5.1.1) 

建立 的 联系 。 若 G(s) 秩 为 m 则 其 逆 函 数 CG)? 存在 , 且 有 
一 对 应 的 特征 方程 

AGs) = detlgIn — GC)™] = 0 (5.1.2) 
EAG O HEL RAH s 的 有 理 函 数 系数 的 多 项 式 ， 则 式 
(5.1.2) 确定 了 代数 函数 (G3 i= 1, 2。.…，1 的 集合 ,以 
后 称 其 为 逆 特 征 增益 函数 。 若 4(8, s) 在 s 的 有 理 画 数 域 上 
不 可 约 , 则 逆 特 征 增益 函数 g(?) A eO 一 样 ,有 一 由 页 复 
* 平面 以 适当 方法 粘 接 的 Riemann 曲面 为 其 定义 域 。 事 实 上 
由 于 矩阵 的 特征 值 是 其 逆 的 特征 值 的 个 数 , 有 
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> 1 
O= 5 (5.1.3) 


因而 8(:) 5 g(9) 具 有 相同 的 分 支点 从 而 也 具有 相同 的 Rie 
` mann 曲面 作为 定义 域 . 
逆 特 征 增益 函数 是 一 基础 ,在 其 上 可 以 建立 着 Nyquist 稳 
定性 判 据 , 


52 极点 /零点 关系 


在 这 一 节 将 建立 一 些 在 广义 逆 Nyquist 判 气 证 明 中 将 要 
用 到 的 极点 /零点 关系 
由 于 矩阵 的 特征 值 是 其 地 的 特征 值 的 倒数 ， 于 是 逆 特 征 
增益 函数 集合 的 极点 就 是 特征 增益 函数 集合 的 零点 ， 且 反之 
DA OL Teese ARIMARA G) 10), 
因此 有 
Pes) = zes) (5.2.1) 
5 
RC) = Pal) E= 1,2, el (5.2.2) 
由 第 3.3-2 节 已 清楚 ， 对 开 环 增益 矩阵 G(s) 说 来 ,存在 
Smith-McMillan 正则 型 [2] 使 
Gs) aa H¢(s)M (s)Jc(s) (5.2.3) 
其 中 Hels), Jels) 均 为 mX m 单 模 态 矩阵 ，Me(s) 具有 
Smith-McMillan 型 , 它 是 
[as AG 
Ms) = dag ES gait ty, rae (5.2.4) 
于 是 G(s) 的 极点 与 零点 定义 为 


Pfs) 一 TI bi(s) (5.2.5) 
i=1 


. 6l» 


m 


zs) = [I ss) (5.2.6) 
i=l 
的 零点 。 现 设 G(s) 秩 为 m, 则 对 式 (5.2.3) 求 逆 有 
G(s) = Jos) Mos) *He(s)* (5.2.7) 
38 FA mE E, 其 定义 为 
0 Orr-0 1 
00-.-1 0 
E=|: : it: :|=E- (5.2.8) 
0 
1 0 *s*0 0 


可 将 GO) KSA 
G(s)* = Jels) EEMGo(s) EEHc(s)™! (5.2.9) 


或 
G(s)" = H&(s)M&(s)J&(s) (5.2.10) 
其 中 
H&(s) = Jo(s)1E (5.2.11) 
Jé(s) = EH<s)* (5.2.12) 
和 


M$(s) = EM (s)7E 
= diag [$2 als) ... hs) 
aus [S Set» BES] 24s) 
CEE G(s) 的 Smith-McMillan 型 ,将 传递 函数 矩阵 极点 
与 零点 的 定义 应 用 于 Gs)", W 


Pes) = Te = zc(s) (5.2.14) 


26s) = 1 pils) = Pels) (5.2.15) 
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其 中 CG), 28s) 表示 GOT 对 应 的 极点 与 零点 多 项 式 。 

在 第 3.3-3 节 给 出 了 G(s) 的 极点 、 零 点 多 项 式 与 特征 
增益 函数 集 {g;(s):i = 1, 2，…, 1} 的 极点 、 零 点 之 间 的 关 
RA 


Pols) = eos) [| bios) = eds) JE Pei(s) (5.2.16) 
f= i= 
和 
L 1 
zc(s) = ec(s) Il bi,(s) = els) II zei(s) (5.2.17) 
izi i=1 


如 果 将 整个 特征 增益 函数 集合 的 极点 、 零 点 多 项 式 表 示 为 
Pels) 与 xz(*)， 而 对 应 逆 特 征 增益 函数 集 的 用 PX《;) 与 25(s) 
表示 , 则 关系 式 (5.2.16) 5 (5.2.17) 以 及 关系 式 (5.2.14) 与 
(5.2.15) 联 合 可 给 出 

Pols) = 26(s) = ec(s)Pe(s) = ecls)az(s) (5.2.18) 
与 

zc(s) = Pas) = ccs)zsls) 
= ec(s)pzH(s) (5.2.19) 

这 些 极点 /零点 关系 在 以 后 证 明 逆 Nyquist 稳定 性 判 据 时 将 会 
HFI. 


5.3 ” 逆 特 征 增益 轨 线 一 一 广义 逆 Nyquist 图 


逆 特 征 增益 扫 线 是 在 复 平面 上 用 开 环 增益 矩阵 GCs) 的 
每 个 特征 值 的 倒数 在 s 沿 Nyquist DD 回路 按 标准 ( 顺 时 针 ) 方 
向 走 过 时 描画 出 的 轨 线 。 用 计算 特征 增益 轨 线 的 算法 可 以 计 
算 逆 特征 增益 轨 线 〈 第 3.4-1 节 )， 但 步子 (iv) 应 改 为 将 
aCe.) 的 倒数 标 在 复 平面 上 . 

现在 提出 并 证 明 广义 逆 Nyquist 判 据 。 
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5.4 广义 逆 Nyquist 稳定 性 判 据 


若 一 般 反 馈 结构 中 子 系统 都 用 状态 空间 型 描述 ， 则 下 述 
提 法 的 判 据 是 可 用 的 . 

提 法 1 一 般 反 馈 结构 是 闭环 稳定 的 , 当 且 仅 当 : 

Cla) 逆 特 征 增益 轨 线 集合 沿 反 时 针 方 向 绕 过 临界 点 
(HA + 10) 的 圈 数 净 和 与 逆 特 征 增益 轨 线 反 时 针 绕 过 原点 的 
疾 数 净 和 之 差 等 于 G(s) 在 右 半 平 面 的 极点 数 . 

O) 逆 特 征 增益 轨 线 不 通过 临界 点 (一 十 j0). 

(3a) 逆 特 征 增益 轨 线 通过 原点 的 分 支 数 等 于 G(s) 在 虚 
轴 上 的 极点 数 ， 且 

(4) FARES, B,C, D) 的 矩阵 4 的 特征 值 中 对 
应 系统 不 可 观测 与 不 可 控 模 的 特征 值 均 在 左 半 平 面 ， 该 开 环 
系统 是 指 输入 为 第 一 子 系统 的 输入 ， 输 出 为 第 个子 系 统 的 
输出 所 对 应 的 部 分 。 

条 件 Cla) 与 (3a) 可 以 改换 为 

(ib) 逆 特 征 增益 轨 线 绕 临 界 点 (一 & + 70) 沿 反 时 针 方 
向 因数 的 净 和 等 于 G(s) 在 右 半 平面 的 零点 数 . 

Gb) 逆 特 征 增益 轨 线 过 无 穷 远 的 分 支 数 是 G(s) 在 虚 
轴 上 的 零点 数 。 

如 果子 系统 均 由 它们 的 传递 函数 矩阵 表 出 ， 或 已 知 每 个 
子 系统 均 没有 位 于 右 半 平 面 或 虚 轴 上 的 不 可 观测 或 不 可 控 
模 , 则 判 据 可 采用 下 述 提 法 。 

提 法 2 一般 反 馈 结构 闭环 稳定 , 当 且 仅 当 : 

(Cla) 逆 特 征 增益 轨 线 集合 绕 过 临界 点 (一 & 十 10) 反 时 
针 轿 数 的 净 和 与 逆 特 征 增益 轨 线 绕 过 原点 反 时 针 圈 数 净 和 之 
差 应 等 于 Gs), Gs), +++, Gs) 在 右 半 平 面 极点 之 总 数 ， 
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O) 逆 特 征 增 益 轨 线 不 通过 临界 点 《一 + j0). 

Ga) 逆 特 征 增益 轨 线 经 过 原点 的 分 支 数 等 于 GOO ER 
轴 上 的 极点 数 ， 

Ga) TURRE 1 的 条 件 《3b) 来 代替 .而且 若 子 系统 
均 为 方 的 , 即 有 相同 数目 的 输入 数 与 输出 数 ， 则 条 件 (1a) 可 
以 代 之 为 

《1b) 逆 特 征 增益 轨 线 绕 过 临界 点 (一 人 十 70) 的 反 时 针 
圈 数 的 净 和 等 于 GCs), +++, Gs) 在 右 半 平 面 零 点 的 总 数 ， 

注意 , 若 条 件 (2) 或 条 件 (3a)/ (3b) 不 满足 , 则 闭环 系统 
有 一 个 或 多 个 极点 在 虚 轴 上 ， 因 而 它 不 能 是 输入 -输出 稳定 ， 
但 平衡 位 置 仍 可 能 稳定 . 

对 DD 矩阵 为 零 的 严格 正则 系统 ， 其 逆 特 征 增益 轨 线 当 s 
趋 于 无 穷 时 将 趋 于 无 穷 ， 因 而 * = co 是 逆 函 数 的 一 个 极点 . 
这 样 在 实际 上 为 了 得 到 封闭 曲线 仅仅 沿 虚 轴 描 画 是 不 够 的 而 
必须 沿 整个 乙 回 线 描画 逆 特 征 增益 回 线 ， 用 这 样 的 办 法 采用 
逆 特 征 增益 回 线 绕 过 临界 点 与 原点 的 园 数 的 净 和 才 可 能 得 
到 ， 作 为 一 例 , 在 图 18 上 画 出 了 


G(s) = 


1 
Gr (5.4.1) 
的 逆 特 征 增益 轨 线 . 

在 稳定 性 判 据 的 两 个 提 法 中 ,如 果 应 用 (la) 而 不 是 (tb) 
时 , 则 D 回 线 中 离开 虚 轴 那 一 段 的 描画 是 不 必要 的 ,这 是 因为 
D 回 线 上 大 圆周 那 一 段 对 应 的 逆 特 征 增益 轨 线 绕 过 原点 与 临 
界 点 的 圈 数 正 相同 ， 因 而 彼此 可 以 消去 . 因此 当 应 用 条 件 
Cla) 时 一 个 有 用 的 法 则 是 用 一 个 大 半圆 连接 轨 线 散 开 的 端 
点 而 不 管 可 能 存在 的 外 面 的 绪 图 数 。 究 竟 是 用 右 半 平 面 还 是 


堪 半 平面 进行 半 个 大 圆 连接 ,应 考虑 这 样 的 事实 , 即 , 8(s) 的 
Riemann 曲面 上 形成 的 Nyquist 曲线 集合 〈 粘 成 了 一 合适 的 
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Nyquist D 回 线 集合 ) 在 平面 的 映射 SC) 应 是 保 角 睦 射 ， 这 
意味 着 车 我 们 想象 两 个 人 ，X 与 Y，X 沿 着 定义 域 中 曲线 
行走 ，Y 沿 着 对 应 象 曲线 行走 ， 并 在 每 一 步 X 确 定 对 应 的 象 
AY, RAX 转向 它 自己 的 右边 ，Y 也 将 转向 它 自己 的 右 
边 。 


9 ngEoq] 


图 18 GG) 一 SSS AEH AS AR 


1 
G+1)° 
(a) 原点 附近 小 范围 ; (b) 大 范围 


. 666 


5.5 广义 逆 Nyquist 稳定 性 判 据 的 证 明 


对 图 3 所 示 的 一 般 反 馈 结构 ,其 闭环 传递 函数 矩阵 RCs) 
由 关系 式 


RCs) = &G(s) [Im + RG(s) 17 (5.5.1) 
所 给 出 , 它 的 逆 应 为 
Rs) = a GCs) + In (5.5.2) 


若 用 {7,(s):i = 1,2, «0+, 表示 确定 RCO) 特征 值 的 代数 
函数 ,将 特征 值 移动 定理 [3] 运 用 于 式 (5.5.2), 则 有 


* hed 
r; = + li 1,2,-++,2 (5.5.3) 


若 用 首 一 多 项 式 PG), ee) 表示 对 应 代数 函数 {7,5): ;一 
1, 2，.….， H 的 完全 集 的 极点 与 零点 多 项 式 , 则 有 
PFs) = p¥(s) (5.5.4) 
MER KOO) 后 乘 方 程 (5.5.2), 则 有 
RR(s)*G(s) = Im + &G(s) = FG) (5.5.5) 
这 是 系统 的 返 差 矩阵 , 取 行 列 式 则 有 


ee P(e) = pr SERO] y ORG) 
det FO) = Gel” y OTO 


其 中 > 是 与 * 无 关 的 常数 ， 但 是 


二 92zG) g a4) 
det F(s) ao) 8 70) (5.5.7) 


Beth p(s), 2/(s) 是 代数 函数 {fi(s); i 二 1,2, 06+, Tb 集合 
的 首 一 极点 与 零点 多 项 式 , 因 而 
p=7 (5.5.8) 


且 
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zs) = mat oe ps) (5.5.9) 


由 方程 (4.2.18) 与 (5.1.21) 有 


PICs) = Pls) = 23s) (5.5.10) 
因而 结合 方程 (5.5.4) 与 (5.5.9) 则 有 
2(s) 一 z%(s) (5.5.11) 


第 四 章 的 方程 (4.2.19) 表 明 闭 环 特征 多 项 式 可 以 表示 为 
CLCP(s) = PACs )ec(s) I] LAO) (5.5.12) 


或 用 本 章 采 用 的 记号 等 价 地 写成 
CLOP(s) = Pals ec(s)z,/(s) (5.5.13) 
因而 联系 到 式 (5.5.11) 则 有 
CLCP(s) = Pals ec(se*(s) (5.5.14) 


这 样 就 可 推出 闭环 系统 稳定 的 充 要 条 件 为 

G) els) 一 0 仅 有 左 半 平面 的 根 ; 

Gi) a*#(s) 一 0 仅 有 左 半 平面 的 根 ; 

Cii) ps) = 0 仅 有 左 半 平 面 的 根 。 

现 设 如 图 16 所 示 的 Nyquist D 回 线 画 在 未 粘 成 Riemann 
曲面 R* 前 的 那些 对 应 的 复 平面 上 ， 在 这 个 Riemann 曲面 上 
7 已 定义 ， 考 虑 对 应 nC) 的 第 ;个 曲面 RY, 一 且 曲 面 已 经 
粘 成 , 则 Nyquist’ D 回 线 就 形成 一 闭合 的 Jordan 曲线 集 , 它 
在 R* 上 包围 右 半 平面 区 域 ， 这 样 扩展 的 辐 角 原理 就 可 以 应 
用 到 每 个 这 样 的 右 半 平 面 区 域 之 上 。 对 每 一 个 特定 的 右 半 平 
面 区 域 说 来 ,( 它 未 必要 单 连通 但 需 以 Nyquist D 回 线 组 成 其 
边界 ) 都 可 以 求 得 7G) 在 该 区 域 零点 与 极点 数 之 差 , 将 等 于 
对 该 特定 区 域 的 边界 曲线 在 WO) 之 下 其 象 在 C ( 复 7 平面 ) 
上 绕 过 原点 的 顺 时 针 圈 数 ， 若 考虑 所 有 右 半 平面 区 域 并 对 其 
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每 一 个 运用 扩展 的 辑 角 原理 , 则 有 
NGCris0) = Z% — P$ (5.5.15) 
其 中 


O NG), 0) 是 在 C 上 用 象 曲线 集合 绕 过 原点 的 顺 时 针 
赂 数 的 净 和 ,这些 象 曲 线 是 映射 7(:) 对 粘 合成 的 RY 上 的 曲 
线 集 的 象 ,该 曲线 集 是 在 R* 未 粘 好 时 先 画 出 而 后 随 RE ARR 
而 形成 的 ; 

Gi) 2* 是 1G) 在 右 半 平 面 的 零点 数 ; 


Gi) PEE rs) 在 右 半 平面 的 极点 数 。 
如 果 我 们 用 同样 的 方法 考虑 对 应 的 逆 特 征 增益 函数 


8(s)， 则 可 以 求 得 


(Ei, 0) = Z}, — P}, (5.5.16) 
其 中 
Ci) N(%，0) 是 用 对 应 和 (5) 的 逆 特 征 增益 雪线 在 C 上 
绕 过 原点 的 顺 时 针 圈 数 的 净 和 3 


Gi) Zh 是 BGs) 在 右 半 平面 的 零点 数 ; 

Gii) Ph 是 &(s) 在 右 半 平面 的 极点 数 . 

合并 方程 (5.5.15) 与 (5.5.16) 可 以 给 出 

NGis 0) 一 NG 0) = Z$ — PL — Z% + Pt, (5.5.17) 
运用 方程 (5.5.4) 则 有 

N(ri0) — NG, 0) = Z$— Z$% (5.5.18) 
SEAM (rei 1,2, epl} bai = 1,2, 
Ty 的 整个 集合 , 则 有 
e69: 


1 * 1 
SNG: 0) — XI NG, 0) 
i=l i=1 


L i 
=>) 24-3) 24, (5.5.19) 
i=1 


izi 


对 闭环 稳定 性 的 条 件 (i) 可 以 等 价 地 说 是 O) 在 右 半 平面 
或 虚 轴 上 没有 零点 ,这 样 上 式 就 可 改写 成 


£ * 1 * i 
,NMCr 0) — >) NCE, 0) = — DS} Zt, (5.5.20) 
f=) i=l i=l 


加 上 {7s):i = 1,2, 在 虚 轴 上 没有 零点 ， 因 此 闭环 
系统 稳定 性 的 必要 充分 条 件 就 可 改写 为 

G) els) 二 0 没有 右 半 平 面 的 视 ; 

Gi’) ecls) = 0 没有 虚 轴 上 的 根 ; 


上 x, ‘ * t 
> Gi?) D NC, Ò) — NG, 0) = — >} Zis 
i=l i=] i=1 


Giv) WON =1, 2,1) 在 虚 轴 上 没有 零点 ; 
Cv) eds) 仅 有 左 半 平 面 的 零点 . 
极点 ,零点 闻 的 关系 式 (5.2.18) 将 引导 我 们 去 考虑 把 条 件 
G) 与 〈 证 ) 合并 成 一 单一 的 等 价 条 件 : 
SING 0) PNG, 0)=—Pe (55.21) 
fet 


其 中 Pc 是 G(s) 在 右 半 平面 的 极点 数 。 在 条 件 (i), Gv) 与 
Cv) 成 立 的 前 提 下 , 式 (5.5.21) 可 以 作为 闭环 稳定 性 的 充分 必 
要 条 件 , 其 证 明 如 下 : 

从 关系 式 (5.2.18), 可 以 有 


t 
Po=et+ >) Zi (5.5.22) 
i=l 


se 70. 


aA 


Sh e Bes) 在 右 半 平面 的 极点 数 , 将 它 与 式 (5.5.19) 合 并 
就 有 

> NG: 0) — 5 NCE, 0) = DY Z3 + e — Pe (5.5.23) 

2 

ING 0) — E NG, Op 
则 由 方程 (5 3 23) 就 可 推出 
> Ztx0 或 ce0 
由 此 可 知 条 件 (5.5.21) 对 闭环 稳定 是 必要 的 . 
对 充分 性 , 设 


l * t x 
DJ NCri, 0) — >) NCB, 0) = — Po 
i=l 1:=1 
则 由 方程 (5.5.23) 可 推 知 
5 zx 一 一 0 


因而 系统 是 闭环 稳定 的 , 即 条 件 (5.5.21) 的 充分 性 成 立 . 
这 样 就 可 以 将 闭环 系统 稳定 性 的 充分 必要 条 件 表述 如 
F: 


E) SENG, 0) — SENG, 0) = —Pes 
f=] f=1 

GP) Cryst 一 1,2, oo, 中 在 虚 轴 上 没有 零点 ; 

Gii’) ecls) 一 0 在 虚 轴 上 没有 根 ; 


Gv") es) 仅 有 左 半 平面 零点 . 
由 方程 (5.5.3) 可 清楚 地 知道 


i 4 1 ry 
SONG: 0) = $, NCB, —4) (5.5.24) 
i=l isi 


于 是 条 件 (i”) 将 等 价 于 稳定 性 判 据 提 法 1 中 的 《1a). 


同样 由 方程 (5.5.3) 可 知 C) = 1，2，.…，1 EME 
轴 上 有 零点 ， 当 且 仅 当 逆 特 征 增益 轨 线 经 过 临界 点 (k+ 
10), 因此 条 件 Cit”) 等 价 于 稳定 性 判 据 提 法 1 UREO). 

由 极点 \ 零 点 关系 趟 (5.2.18) 与 (5.2.19), 条 件 Gii”) 显然 
等 价 于 稳定 性 判 据 两 提 法 中 条 件 Ga) 与 (3b) 中 任 一 个 。 

这 样 为 完成 稳定 性 判 据 提 法 1 的 证 明 ， 留 下 的 是 证 明 条 
件 (la) 与 (1b) 等 价 。 为 做 到 这 一 点 我 们 将 证 明 当 条 件 (if)。 
Git”) 与 Gv") 成 立时 ,条 件 


SING, 0) = SNE, A = —Ze (5.5.25) 


是 闭环 稳定 的 充分 必要 条 件 , 其 中 Zs% G(s) 在 右 半 平面 的 
FAR. 
由 关系 式 (5.1.22) 可 有 


Zoe + XP} (5.5.26) 
i=] 
其 中 Pt, Be) 在 右 半 平面 的 极点 数 。 但 由 方程 (5.5.3) 知 
C) 与 8(s) 的 极点 相同 ,因而 
Zo=et 5 P% (5.5.27) 


f=} 
若 联合 方程 (5.5.27) 与 (5.5.15), 其 中 后 者 对 每 个 ;一 1, 2，…， 
1 均 成 立 , 则 有 
SING, 0 =J Zt +e— Ze (5.5.28) 
f=1 i= 
为 建立 条 件 (5.5.25) 的 必要 性 , 设 


SNG 0) = 一 Ze 


。72 。 


则 由 方程 (5.5.28), 上 式 可 以 导致 
> Z* <0 和 /或 < 关 0 


而 由 此 则 可 断言 条 件 (5.5.25) 对 闭环 稳定 性 是 必要 的 .” 
对 充分 性 , 设 


DING 0) = 一 ze 
则 由 方程 (5.5.28) 可 有 


t 


zeae =0 


i=1 


AMERA G”), Gi) 与 (iv”) 成 立 的 条 件 下 可 推 知 系统 
是 闭环 稳定 的 。 从 而 式 (5.5.25) 的 充分 性 得 证 . 
这 就 完成 了 广义 逆 Nyquist 稳定 性 判 据 提 法 1 AEBS, m 
稳定 性 判 据 的 提 法 2 应 用 于 这 样 的 系统 [参见 第 四 章 ]， 
在 其 中 或 
Gi) Pals) = Px(s) (5.5.29) 
或 
ii) Pals) = PCs) Pa,(s)Pa,(s)** *Pa,(s) (5.5.39) 
其 中 {pa(s):i = 1,2, ++, 7} 的 零点 全 在 左 半 平 面 ， 
在 上 述 这 些 情形 ,可 以 联合 提 法 1 的 条 件 (Ca) 与 (4) RAM 
一 的 等 价 条件 
1 * t * a 
DNG =) — >) NG 0) = 一 之 Po, (5.5.31) 


i=] 
在 条 件 (2) 与 (3) 已 成 立 的 前 提 下 ,条 件 (5.5.31) 对 闭环 系 
统 稳 定性 的 充分 性 与 必要 性 可 以 证 明 如 下 . 
由 方程 (4.2.37 ) 与 (4.2.38) ,有 
© 原文 误 为 式 (5.5.28) 对 闭环 稳定 性 是 必要 的 .一 一 译 者 注 
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> Po, =e + 5 P;, + Pa (5.5.32) 
i=} i=l 


但 {f(s):i = 1, 2, +++ D ORAS {e:ls):i = 1,2,*… ,1} 


BBABA, MBBS {2C6): i= 1,2, 4,1} 的 零点 要 
同 。 因 而 方程 (5.5.32) 可 改写 成 


六 可 -六 (5.5.33) 
i=l iel 
现在 式 (5.5.33) 可 以 与 式 (5.5.19) 联 合 给 出 
PN, 0) — ENG, 0) 
ist isi 


t a 
=D ZH+etP— D>) Po (5.5.34) 
i=1 i=l 
应 用 式 (5.5.24) 它 可 变 成 
SG, =) — SG, 0) 
i=] 


fet 


上 A 
= $, Zh +e +Pi— D>) Po; (5.5.35) 
i=) 


m 
为 确立 条 件 (5.5.31) 的 必要 性 , 设 
DME) 一 NG;, 0) = -5 Pa 
则 由 方程 (5.5.35), 可 有 
2 e 20 R P20 


或 它们 的 任何 组 合成 立 , 于 是 系统 将 是 闭环 不 稳定 。 这 就 可 
以 断言 条 件 (5.5.31) 对 闭环 稳定 性 是 必要 的 。 
对 于 充分 性 , 设 
。74 。 


SG, -b- YG, D=- Pa 


i=1 


则 贝 方程 (5.5.35 ) 可 以 有 
[i 
D Zhe =P, = 0 
从 而 系统 在 已 满足 《2), G) 的 条 件 下 是 闭环 稳定 的 .因此 式 
(5.5.31) 的 充分 性 得 以 确立 ， 

于 是 为 完成 稳定 性 判 据 提 法 2 的 证 明 ， 留 下 的 是 要 证 明 
条 件 Cla) 与 (lb) 的 等 价 性， 为 此 ,我们 将 证 明 当 条 件 (2) 
与 (3) 成 立时 ,条 件 (1b) 即 

5 NCr 0) = Dv, 一 人 ) = -5 Zo (5.5.36) 


是 闭环 稳定 的 必要 充分 条 件 ， 其 中 Zo 是 G6; CEH 半 平 面 的 , 
零点 数 
方程 (5.5.15) 对 于 j = 1, 2，……，! HRL BIE 


5 NGis 0) = >i NCE, —k) 


= > zz- DP (5.5.37 

又 由 于 {C95i = 2 n DIRAS EO i= 1, 2, 
eo N) 的 极点 相同 ,方程 (95.37) 可 以 改写 成 

EN -D= 5zi- z (5.5.38) 

EIRA 【G1(5):i 一 1,2,- A) 每 个 都 是 方 的 ， 则 


P:( =P.) 是 在 G(s) 形 成 时 零点 \ 极 点 相 消 中 在 右 半 平 面 的 零 
点 (或 极点 ) 数 。 因 此 , 若 G(s) 在 右 半 平 面 零点 数 给 定 为 


i 
Zo=e+ Ď Z; (5.5.39) 
gi 


1=1 
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出 必 有 
Si =e+ È Zu 十 Ps (5.5.40) 
或 由 于 LO) 的 极点 恒 等 于 C) 的 零点 , 则 
Sza = e+ SP; +P, (5.5.41) 
于 是 方程 (5.5.38) 与 (5.5.41) 可 联合 起 来 给 出 
SING’ = pa Ze +e +P, 一 z Ze (5.5.42) 


为 确立 式 (5.5.36) 的 必要 性 , 设 
ENG, W = — > Zo 
再 由 方程 (5 5.42), 就 可 扒 得 


t 
J Z~ 0 R ex0 R P0 


iel 
或 上 述 的 任何 组 合 , 从 而 系统 将 是 闭环 不 稳定 的 。 这 就 证 明 
了 条 件 (5.5.36)* 对 闭环 稳定 是 必要 的 . 
对 于 充分 性 , 设 
t * h 

ENE, k) 一 一 >) Za 

be #=1 
则 由 方程 (5.5.42 ) ,就 可 推 得 

5 Z*# 一 e 一 Pi 一 0 


因而 在 (2),(3) 条 件 已 成 立 的 前 提 下 ,系统 是 闭环 稳定 的 ， 从 
而 条 件 (3.3.36) 的 充分 性 得 以 确立 . 
.这 就 完成 了 广义 逆 Nyquist 稳定 性 判 据 的 证 明 。 
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为 说 明 稳 定性 判 据 , 考 虑 在 第 三 \ 四 两 章 中 用 过 的 例子 
G(s) = G(s) 一 


1 s—i s 
1.25(s + 1)(s + 2) E y | 
开 环 增益 矩阵 G(s) 的 极点 与 零点 多 项 式 分 别 为 
Pos) =(s+1Xs +2), Zs) =1 
FRA 
Po =0, Zo=0 
其 逆 特 征 增益 轨 线 如 图 19 所 示 


图 19 逆 特 征 增 益 轨 线 


由 于 不 存在 不 可 观 或 不 可 控 模 ， 因 而 提 法 1 与 提 法 2 均 
可 直接 运用 。 实 际 上 ,由 于 有 效 地 仅 有 一 个 子 系统 ,两 个 提 法 
在 这 里 将 无 区 别 . 

SES MB Nyquist 稳定 性 判 据 列举 的 条 件 ， 可 以 发 现 
对 于 
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—1875<k <0 
0<k< 1.25 


5 


2.5 < k < œ 


系统 是 稳定 的 。 在 增益 控制 变量 & 上 所 加 的 限制 同 第 4.3 节 
的 一 致 ,在 那里 用 的 是 广义 Nyquist 判 据 . 


有 意义 的 是 指出 逆 Nyquist 式 的 判 据 现 在 已 被 Mess 与 


Rapp 应 用 于 建立 多 变量 非 线性 反馈 系统 的 稳定 性 判 据 [6]. 


[1] 
[2] 
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第 六 章 多 变量 根 轨迹 


Evans 方法 |1,2,3] 是 在 线性 时 不 变 单 输入 、 单 输出 反馈 
系统 分 析 与 设计 中 的 一 个 很 好 的 图 解 方法 ， 它 可 对 系统 闭环 
极点 作为 增益 控制 变量 的 函数 进行 估计 . 将 它 推广 至 多 变量 
情形 是 有 意义 的 ， 而 且 从 状态 空间 观点 去 考虑 特征 频率 轨 线 
的 渐 近 性 质 [4,5] 和 特征 频率 轨 线 的 出 发 角 与 抵达 角 [6] 这 两 
方面 , 就 可 找到 达到 这 一 目的 的 路 子 。 本 章 发 展 了 一 个 在 代 
数 函数 理论 中 已 有 结果 的 方法 ， 在 这 里 特征 频率 轨 线 的 渐 近 
性 质 与 出 发 角 及 抵达 角 均 可 由 系统 的 Laplace 传递 函数 矩阵 
来 确定 [7,8,9]. 

它 同 样 可 以 阐明 该 方法 如 何 能 用 来 寻求 多 变量 时 不 变 线 
性 调节 器 ， 最 优 闭环 极点 在 品质 判 据 中 输入 的 权 趋 于 零 时 的 
渐 近 性 质 [10] 

对 于 具有 少数 输入 的 系统 ,由 于 计算 简单 可 以 进行 手 算 ， 
因而 该 方法 特别 有 用, 


6.1 理论 回顾 


对 图 3 所 示 的 一 般 反 馈 结构 ， 在 第 3.3-4 节 已 表明 特征 
频率 轨 线 (多 变量 根 轨迹 ) 是 特征 增益 函数 gC) 的 180° 等 位 
相 曲 线 , 其 中 g(s) 是 由 方程 

Alg, s)Adet[gI m — G(s)] 一 0 (6.1.1) 
或 方程 
P(g, s) = bols)Alg, s) = 0 (6.1.2) 
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定义 的 ,后 者 已 用 * 的 有 理 函 数 系数 的 最 小 公分 母 乘 过 .为 叙 
述 简单 且 对 通常 传递 函数 常 满足 的 情形 , 设 A(g,*) 在 :的 有 
理 函 数 域 上 已 为 不 可 约 多 项 式 、 特 征 增益 函数 g(s) 与 增益 控 
制 变量 & 间 的 关系 用 表达 式 


g 一 一 天 (6.1.3) 
联系 ,因而 在 式 (6.1.2) 中 替换 掉 & 就 有 
OC— K, s) = KTR, s) = 0 (6.1.4) 
且 
rk, s) = 0 (6.1.5) 


用 正 实数 不 对 * 求解 就 确定 了 闭环 极点 * 依赖 于 增益 控制 变 
量 & 的 关系 。 因 此 ,除了 可 能 的 单 点 轨迹 外 ,这 样 一 个 依赖 关 
系 的 图 形 组 成 了 给 定 系统 的 特征 频率 轨 线 。 TER, BREE 
单 点 轨 线 则 方程 (6.1.5) 直 接 等 价 于 
CLCP(s)4 det [s1, — A.} 
= detls1s — S(—k")] = 0 (6.1.6) 
且 特 征 频 率 轨 线 具有 个 分 支 
特征 频率 轨 线 中 的 每 一 个 分 支 在 理论 上 均 可 表示 为 下 述 
形式 : 
sik) = wih) + inik); F122, +n (6.1.7) 
其 中 j = MV 一 1, 而 标 码 i 表示 分 支 的 序号 。 轨 线 第 ， 个 分 支 
TE so = sko) 的 切线 ， 定 义 为 通过 点 o 和 曲线 上 另 一 点 
si = s, (kot 8k) 的 直线 在 Ok — 0 时 的 极 跟 位 置 。 复数 sso 
可 以 用 由 s% 至 “的 向 量 表示 (图 20), 而 对 应 Cr — s) + ôk 
的 向 量 , 当 ok > 0 时 与 一 具有 同样 的 方向 。 因 此 可 推 得 
对 应 于 


d ds(k) | E So 
(Ro) dk lr lim 
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= im Ho + 8k) — Ch) ED 
bho dk 
的 向 量 , 它 是 分 支 在 % 的 切线 ,而 这 一 向 量 和 正 实 轴 之 间 的 夹 
角 为 
6 = 辑 角 (5(h)) (6.1.9) 
式 〈6.1.9) 对 确定 特征 频 
率 轨 线 的 渐 近 性 质 以 及 轨 线 
的 出 发 角 与 抵达 角 来 说 是 基 
本 的 .对 =0, 方程 (6.1.3) 推 
知 特征 频率 轨 线 由 开 环 系统 的 
极点 出 发 ， 因 而 轨 线 分 支 从 这 
一 极点 的 出 发 角 将 由 式 (6.1.9) MM 
中 心 一 0 来 确定 . 对 4 = oo， 5 
方程 (6.1.3) 表 明 特 征 频 率 轨 线 
将 终止 在 系统 的 零点 上 。 因 而 。 。 图 20 特征 频率 纵 线 切线 
对 于 轨 线 终止 在 有 限 零点 的 分 PEHAR 
支 说 来 ,公式 (6.1.9) 在 名 = om SHH EAA AE UNAG fa 
度 。 着 一 个 分 支 终止 在 无 穷 零 点 , 则 如 二 oo 的 公式 (6.1.9) 将 
给 出 该 分 支 渐 近 线 与 正 实 轴 的 夹 角 . 
式 (6.1.9) 不 能 认为 是 很 适用 的 ， 这 是 由 于 一 般 不 能 明显 
地 找到 区 别 特征 频率 轨 线 分 支 的 表达 式 ， 但 是 ， 在 代数 函数 
论 中 有 着 一 种 实际 构造 方法 [7,8,9], 它 可 对 一 代数 函数 在 给 
定点 邻 域内 给 出 该 函数 的 分 支 的 级 数 表达 式 .这 一 方法 可 直 
接应 用 Newton 医 的 方法 来 做 以 求 得 级 数 中 下 一 个 最 有 效 的 
项 .因此 用 一 次 Newton 图 ,就 可 求 得 特征 频率 轨 线 在 极点 或 
零点 (有 限 或 无 限 ) 处 的 近似 : 
sR) =a + ok (6.1.10) 
其 中 a 与 5 为 复数 。c HAMM, He HRI k ER 
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Si= S (kot Ok) 


是 被 理解 的 ”. 对 于 k TESA MAO EE OK, oh 
将 总 是 实 的 ,因而 应 用 式 (6.1.9) 至 近似 式 (6.1.10) 则 有 


0 二 辐 角 (b)+r180° (6.1.11) 
其 中 
<- 站 > 4ard 
1， 当 ee 一 0 
这 些 是 由 对 《微分 而 得 来 的 。 


下 一 节 将 说 明 如 何 得 到 极点 或 零点 附近 的 近似 式 
(6.1.10), 从 而 确定 特征 频率 轨 线 的 渐 近 性 质 以 及 轨 线 的 出 发 
角 与 抵达 和 角 . 
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Newton 图 是 一 个 在 原点 附近 的 图 形 构 造 , 它 可 以 用 来 寻 
求 代数 函数 g(v) 特定 分 支 级 数 表示 中 每 个 最 有 效 的 项 ,其 中 
so) 具有 当 。 趋 于 零 时 其 值 趋 于 零 的 性 质 ， 因 此 , 无 论 是 寻 
求 轨 线 的 渐 近 性 质 还 是 出 发 角 与 抵达 角 ， 其 第 一 步 经 常 是 在 
特征 方程 中 用 适当 的 变量 替换 来 简化 问题 ， 使 对 代数 函数 的 
这 种 分 支 能 求 得 一 近似 ， 对 该 分 支 当 其 独立 变量 趋 于 零 时 亦 
BTE. 

为 了 确定 特征 频率 轨 线 的 渐 近 性 质 ， 我 们 需要 得 到 一 个 
当 A 趋 于 co 时 轨 线 在 点 * 一 co 周围 的 分 支 (或 几 个 分 支 ) 的 近 
Wh. 为 此 ,在 方程 (6.1.2) 中 令 

saat (6.2.1) 

于 是 就 有 


O 此 处 若 «= 8B/7Y，B，Y 为 整数 且 互 相 无 公 因子 , 则 k= (kV, jth 
AVY Hak 的 主根 。 一 一 译 者 注 
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O(g,5) = O(g, 2) = 2 Wg, 2) (6.2.2) 
其 中 9 委 > 是 es) 的 极点 数 。 TETH z = 0 的 任何 邻 域 


A (z 一 0 本 身 被 排斥 在 区 域 之 外 ) 方 程 B(g, s) = 0 等 价 于 
方程 


Plgz) = Ddxygtz” (6.2.3) 
对 于 一 个 严格 正则 系统 , 即 其 中 忆 一 0 MAD HAF 
Yoo = 0 (6.2.4) 


否则 gw 将 为 非 零 ,从 Newton 图 将 会 看 到 , 对 应 的 将 不 会 
渐 近 性 质 . ERE DAER 则 es) 具有 与 极点 数 相同 的 
零点 数 , 这 样 就 不 能 指望 有 任何 闭环 极点 趋 于 无 穷 。 
下 一 步 是 对 亚 (g,， z) 构造 Newton 图 , 形 如 
z= cg" (6.2.5) 
的 近似 将 由 此 得 到 ,其 中 < 为 复数 ,4 为 有 理 实数 .构造 合适 
的 Newton 图 的 方法 如 下 [8]: 
E (8,2) 平面 上 我 们 标 出 了 方程 (6.2.3) 中 dey = 0 所 对 
应 的 点 (x*，?)。 作 为 例子 ,在 图 21 LEHT 
(1— z + 22? — 2523 十 29z4)82 
+ (z — 2227 + 1992? — 2102")g 
— (3323 一 5942") = 0 (6.2.6) 
的 对 应 点 。 将 一 直 尺 放 在 与 水 平 轴 重 合 的 位 置 然后 绕 最 小 & 
轴 点 Po (图 21 上 即 为 点 (2,0)) 作 顺 时 针 旋 转 , 直至 直 尺 通过 
网 格 的 另 一 点 P, BRP, 与 P, 作出 一 直线 。 过 已 由 已 水 
平 指向 铅 直 轴 放置 直 尺 , 然后 绕 P 点 作 顺 时 针 旋 转 直 至 通过 
网 格 的 另 一 点 Po ERP, P,, 作出 又 一 条 直线 . 重复 上 述 过 
程 直至 达到 铅 直 的 > 轴 。 方程 (6.2.6) 的 完全 Newton 图 作 在 
22 E. PP, 连 线 与 铅 直 = 轴 所 夹 锐角 的 正切 确定 了 4 的 
第 一 个 可 能 值 , 而 PP, 连 线 与 铅 直 z 轴 所 夹 锐角 的 正切 确定 
了 “的 第 2 个 可 能 值 , 等 等 。 对 图 22 所 示 的 Newton 图 , 有 
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° 1 了 


图 21 方程 (6.2.6) Newton 图 的 点 ”图 22 方程 (6.2.6) 的 完全 Newton 图 


m=1 5 m= 1/2, 
对 一 特定 指数 wm， 可 以 有 几 个 有 逼近 形式 : 
Zi 一 Cug" (6.2.7) 
注意 若 8 的 某 个 根 由 a 定 出 , 即 py; 是 分 数 , 则 对 应 考虑 的 是 
主根 。 为 了 确定 系数 cus 必须 在 方程 (6.2.3) 的 对 应 网 格 中 ， 
CHRP oP 上 的 点 的 项 中 代入 z= <8g% ， 然 后 使 其 为 零 再 
去 解 结果 方程 ， 设 方程 (6.2.3) 中 有 关 项 的 和 具有 形式 
Proyog og + 十 TT (6.2.8) 
其 中 yo De > ys o 系 一 正 整数 且 小 于 系统 无 穷 零 点 的 数 
目 ， 由 于 其 项 均 对 应 在 同一 线 上 则 可 有 
LT mi WO Font TO ws a AT 
Yon You — Nn Ya Mn 
因此 式 (6.2.8) 表 达 式 的 全 部 项 由 于 用 z= co” 的 结果 替换 而 
成 为 相 类 似 的 ,因此 为 确定 系数 c* 我 们 得 到 一 个 形 如 


™ F 


Pigra” Hore + pan” = 0 (6.2.9) 
的 方程 或 以 oR C 关 0), 则 有 
Dagr tH ts + dey, = 0 (6.2.10) 


. 84. 


CH y 一 MBL 注意 ye EP GP, ERZE. 
因此 它 应 等 于 线段 P,_,P, 在 纵 坐 标 轴 上 的 投影 。 由 于 它 能 使 
我 们 看 到 对 应 指数 wm 究竟 有 多 少 系 数 ， 因 而 这 是 有 用 的 事 
实 ， 同 样 ， 最 后 一 个 线段 达到 铅 直 轴 时 的 点 的 纵 坐 标 将 给 出 
对 考虑 所 有 指数 的 所 有 近似 系数 的 总 数 。 因此 在 最 后 线段 上 
这 一 铅 直 轴 点 的 纵 坐 标 就 是 系统 无 穷 零点 的 数目 . 

在 得 到 方程 (6.2.7) 那样 的 近似 后 ， 由 方程 〈6.1.3) 与 

(6.2.1) 直 到 我 们 所 需 的 形式 

Sit = bukt (6.2.11) 
只 是 一 个 替换 的 简单 过 程 ， 而 这 正 是 特征 频率 轨 线 对 于 4 一 
co 时 有 s 一 co 的 近似 描述 。 如 果 应 用 公式 (6.1.9) 我 们 就 可 以 
得 到 渐 近 线 与 正 实 轴 间 的 夹 角 。 

渐 近 线 正 如 在 第 6.2-1 节 阐 明 的 ,将 组 成 一 个 Butterworth 
结构 。 每 个 模型 均 有 公共 的 截 距 , 它 被 Kouvaritakis 和 Shaked 
[4] WA ZERRA” (multivariable pivot)， 他 们 也 给 出 了 
对 给 定 状态 空间 描述 计算 多 变量 枢 点 的 办 法 。 在 附录 6 中 将 
阐明 如 何 由 特征 方程 Ag, s) = 0 推导 出 多 变 倒 枢 点 ， 


6.2-1 Butterworth 模型 


在 这 一 小 节 证 明了 闭环 极点 走向 无 穷 远 是 沿 着 组 成 几 个 
Butterworth 结构 的 渐 近 线 的 。 证明 是 基于 代数 函数 理论 中 已 
确立 的 结果 [8]. 

考虑 由 方程 (6.1.5) 定 义 的 闭环 特征 多 项 式 , 即 

rT(k,s)=0 (6.2.12) 
其 中 人 允许 4 为 复数 . 设 在 这 一 点 由 方程 (6.2.12) 定 义 的 代 
数 函 数 s(《) 具有 P— FIBA seo 设 s(*) 在 的 邻 域 NN 内 的 
个 分 支 为 {5(R):i 王 1,2,，… Pa etant AR {sk = 
Seti = 1,2, +++ P} 
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设 对 分 支 {5(《):i = P+ 1, +, a} 中 任 一 个 在 N 内 围 
绕 进行 解析 开拓 .并 且 我 们 思 绕 着 临界 点 对 i(k) 解析 
开拓 , 则 在 一 个 绕 圈 以 后 我 们 得 到 函数 {5:(%):i 二 2,3,…, P} 
中 的 一 个 ,例如 5(*), 而 在 另 一 个 绕 圈 以 后 得 到 函数 {5:(*): 
k= 3，4，…… sp} 中 的 一 个 ,例如 ss(《)。 在 有 限 数 v 个 绕 图 
以 后 回 到 原来 的 函数 CX)， 在 这 样 的 情况 下 将 称 函数 

Ck), cts Sn Ck) (6.2.13) 
组 成 了 代数 函数 (A) (E k. BRNA) 的 分 支 的 一 循环 系 . 

Hn<p, WHR Sy,41€)» 然后 对 它 在 六 内 围绕 人 BE 

行 解析 开拓 ,重复 上 述 做 法 .得 到 第 二 个 分 支 的 循环 系 , 记 作 
5y41€R) > sey Sn Ck) (6.2.14) 
最 后 ,所 有 函数 aA), +++» so《%) 均 被 收入 循环 系 . 

当 绕 着 多 重点 转动 时 , 每 个 循环 系 的 函数 按 循 环 律 
依次 转 为 另 一 个 函数 ， 每 个 分 支 的 循环 系 在 k IPREN E 
义 一 个 解析 函数 ， 该 解析 函数 在 邻 域 N 的 点 的 值 的 数目 等 于 
考虑 组 成 这 个 循环 系 的 分 支 数 ， 若 一 个 分 支 就 组 成 一 个 铂 环 
系 时 , 则 对 应 的 解析 函数 是 单 值 的 且 恒 等 于 该 分 支 。 

设 邻 域 N 是 用 个 这 样 的 区 域 块 堆 起 来 的 ,函数 {e(4): 
§=1,2, +++ sn} 中 的 每 一 个 都 有 一 块 ,而 且 使 第 i 页 上 的 点 
表示 为 ky s:《《)] 对 。 又 设 每 个 领域 呈 贺 形 且 具有 一 径 向 由 
临界 点 至 圆 边界 的 割 口 ,如 图 23 所 示 。 则 邻 域 可 以 沿 这 个 割 
口 的 边 进行 连接 , 以 便 在 第 i 页 上 围绕 进行 解析 开拓 ， 第 
i 个 分 支 (i < p) 经 过 制 口 进入 第 (i 十 1) 页 ,… EAFA 
显然 又 可 回 到 第 ;j 页 。 邻 域 的 堆砌 显然 形成 包含 一 页 或 多 页 
的 一 些 循环 如 图 24 所 示 ， 且 每 个 循环 是 定义 在 的 邻 域 N 
上 的 解析 承 数 的 定义 域 ， 

为 表述 由 一 循环 确定 的 解析 函数 ,我们 将 采用 罗马 数 标 . 
例如 对 应 循环 系统 (A) 2 os A) 将 用 符号 a(《) RZ. 
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图 23 临界 点 的 邻 城 N 图 24 对 临界 点 MOBI HAE 
当 我 们 考虑 特征 频率 轨 线 的 渐 近 性 质 时 ， 我 们 关心 的 是 
代数 函数 (CA) 的 这 样 的 分 支 : 它 当 趋 于 无 穷 时 趋 于 无 穷 . 
如 果 对 于 《二 oo 我 们 有 个 无 穷 分 支 ， 则 按 前 面 的 讨论 这 了 
个 分 支 将 安排 在 几 个 循环 系 中 。 每 个 循环 系 由 一 解析 函数 确 
定 , 它 的 形式 由 下 述 定理 给 出 [8, p39]. 
定理 4 代数 函数 A) 在 其 临界 点 = co 的 领域 内 构 
成 一 循环 系 的 那些 无 穷 分 支 在 总 体 上 组 成 一 个 由 一 形式 级 数 
51 = (Oy bk HR? Hee) (6.2.15) 
SHAT BK, Heh» 是 该 循环 中 的 分 支 数 ，4 为 一 正 整 
数 ， 
于 是 在 & 一 co 的 邻 域内 对 应 一 特定 循环 系 的 无 穷 分支 ， 
将 可 以 近似 地 用 
sR) = bk” (6.2.16) 
特征 频率 轨 线 是 和 沿 正 实 轴 由 0 变 至 无 穷 时 对 应 代数 函 
Be (O 的 值 。 因 而 频率 轨 线 的 无 穷 分支 将 出 现在 能 近似 用 
式 (6.2.16) 表 示 的 循环 系统 中 。 在 方程 (6.2.16) HEY RAD v 
个 根 , 我 们 可 以 找到 对 应 循环 系 的 渐 近 线 , 它 们 自己 被 安排 在 


一 个 Butterworth 结构 中 。 此 即 在 复 平面 上 有 个 用 (5) 


v 
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角 进 行 等 分 的 渐 近 线 . 

一 般 说 来 正如 已 经 讨论 的 那样 ， 在 不 一 o 的 邻 域 内 将 存 
在 几 个 分 支 的 循环 系 ， 因 此 特征 频率 的 渐 近 线 将 被 安排 在 几 
个 Butterworth 结构 之 中 ， 


63 出 发 角 与 抵达 角 


为 了 确定 特征 频率 轨 线 的 出 发 角 与 抵达 角 ， 我 们 需要 得 
到 轨 线 在 系统 的 极点 与 有 限 零 点 那里 的 分 支 的 近似 式 。 我们 
首先 若 虑 由 系统 开 环 极点 出 发 的 角度 . 

设 对 于 特征 方程 

O(g,5)=0 (6.3.1) 

在 * 一 8 处 有 一 极点 (或 多 重 极点 )。 对 于 方程 (6.3.1) 可 取 
s =s —p 为 新 的 独立 变量 ， 同 时 用 8 = ad 替换 使 方程 
(6.3.1) 变 成 

Plg, s) = O(d", s +8)=d "sa(d,s)=0 (6.3.2) 
我 们 正在 考虑 的 情况 现 已 化 简 为 对 应 4 = s = 0 的 情形 , 而 
E s = p (B 本 身 排斥 在 外 ) 的 邻 域内 方程 p(g,s) = 0 将 等 
价 于 方程 


ald, s) = DE syd*s'” = 0, Eq = 0 (6.3.3) 
若 对 方程 (4,s) = 0 构造 Newton 图 , 则 得 到 一 个 形 如 
s = ed* (6.3.4) 


的 近似 (或 多 重 极点 情况 下 的 近似 ), 其 中 。 为 复数 , "为 有 理 
实数 ， 由 方程 (6.1.3), 改变 变量 并 作 替 换 g = d7', 就 得 到 了 
从 极点 8 出 发 的 特征 频率 轨 线 的 一 个 分 支 或 几 个 分 支 的 近 
似 ,其 形式 为 

s = B+ bakes (6.3.5) 
若 现 在 应 用 公式 (6.1.9), 则 出 发 角 为 


. 8B 。 


6. 一 辐 角 {ba} (6.3.6) 
现在 我 们 将 考虑 到 达 系 统 有 限 零 点 的 角度 ， 设 特征 方程 
《6.3.1) 有 一 零点 (或 多 重 零 点 ) 在 s= jb. W "= 一 7 为 
新 变量 ,于 是 方程 变 为 
O(g,5)=P(g,5 +7) =Xg,5)=0 (6.3.7) 
这 样 考虑 的 情况 简化 至 g 一 * = 0 的 情形 。 而 在 “一 7 的 邻 
Bir 本 身 排斥 在 外 ) 内 方程 8(g, s) = 0 等 价 于 方程 


Xle, s) = Drsyg's’ = 0, xm = 0 (6.3.8) 
如 果 对 方程 (6.3.8) 构 造 Newton 图 ,一 个 形 如 
s & pg? (6.3.9) 


的 近似 就 得 到 了 ， 其 中 ?是 复数 ,7 是 有 理 实数 。 由 方程 
(6.1.3) 与 变量 改变 ， 我 们 就 得 到 了 抵达 零点 :二 7 的 特征 频 
率 轨 线 的 一 个 分 支 或 几 个 分 支 的 一 个 近似 ,其 形式 为 


ser + bk" (6.3.10) 
如 果 现 在 应 用 公式 (6.1.9), 则 抵达 角 9. 为 
6, = 辐 角 {5。} 土 180° (6.3.11) ` 


注意 n 将 总 是 负 的 ， 因 而 在 公式 (6.3.11) 中 出 现 180° 这 
一 项 .同样 要 注意 的 是 这 一 抵达 角 的 定义 与 通常 的 定义 相差 
180° ,后 一 定义 是 站 在 零点 去 看 轨 线 到 达 时 才 可 看 到 的 方向 ， 


6.4 例 1 


考 菩 开 环 增益 矩阵 : 


G(s) = 5 


S +5 — 28— 445 + 40 
PE ATA 8 — 24s + 16 
+ 79s? + 44s — 868 4 一 49: 十 405 一 32 


它 的 特征 方程 为 
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O(g5s) = (st + 58 — 28? — 445 + 40) 2 
+ (8 + 116s — 432)g — 126? — 2s +2) =0 
G) 确定 渐 近 性 质 
在 特征 方程 中 设 s = 2", 则 得 到 
Tlg, 2) = (1 + 5z — 227 — 442° 十 40z4)8 
+ (z + 1162? — 4322")g 
— 12(2? — 2? + 22) =0 
Ulge, z) 的 Newton 图 画 在 图 25 上 , 由 此 可 知 p= 二 1， 于 是 
近似 为 


zXcg 
算出 系数 。 具有 了 两 个 值 一 二 与 二 《用 第 6.2 节 的 办 法 )。 因 
此 对 zx 进行 回 代 并 且 g = 一 -!, 则 对 特征 频率 轨 线 可 求 得 下 


z 述 近 似 : 
A su 4k 5 sœ —3k, 
当 k->00 
” 因此 特征 频率 轨 线 的 两 分 支 以 
2 对 正 实 轴 的 角度 0° 与 180° 走 
IES. 
\ GD) 求 出 发 角 
系统 有 四 个 开 环 极点 
o 1 k- g s=1, s=2, s=—4+2j 
图 25 (gy =) 的 Newton 图 与 5 一 一 4 一 2 
在 极点 * 一 1 处 


在 特征 方程 中 令 * = 一 1 与 8 一 4 WE 
Ed, s’) 一 (人 十 8 二 18 一 289) 

+ (59 4 3s? + I1gs’ — 315)d 

— 126s? + 1)a? = 0 
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(dss) 的 Newton 图 画 在 图 26 上 ,由 此 可 有 o = 1, 因而 
近似 为 
s œ ed 
系数 e 经 算出 为 一 10.86《 用 第 6.2 节 的 方法 )， 最 后 特征 频 
率 轨 线 在 极点 * 一 1 局 围 的 近似 为 
s œ 1 + 10.86% 

因此 在 极点 * = 1 的 出 发 角 为 0%. 
在 极点 * = 2 处 

在 特征 方程 中 令 =: 2 5g = d, WA 

(d,s )=(s" + 138 + 5257 + 405°) 


s s 
4 
4 
3 
3 
i 2 
1 2 
1 
1 
1 
i A 
L ; 7 
o i 2 r 


图 26 对 全 (ds) 的 Newton 图 图 27 对 (df) 的 Newton 图 


图 28 对 sdyr) 的 Newton 图 图 29 对 Xess) 的 Newton 图 
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+ (5 + 6s? + 128s’—192)d—12(s? + 2s + 2)0? = 0 
Ed, 人) 的 Newton 图 画 在 图 27 上 ， 由 此 可 有 o = 1, 因而 
近似 为 

$ Med 
系数 。, 算出 为 4.8， 最 后 特征 频率 轨 线 在 极点 = 2 的 周围 
具有 下 述 近 似 : 
s=2—4.8k 
从 而 由 极点 = 2 的 出 发 角 是 180°, 
在 极点 ;== 一 4 十 21 处 
在 特征 方程 中 令 =s +4 2i 5 g= d>, WE 
(ays) = Ls" + (—11 + j8)” 

+ (10 — j60)s? + (88 + 7104)s’] 

+ [5 + (—12 + j6)s + (152 一 148)7 

十 (一 912 十 1320)]d 一 12[s 十 (一 10 十 14)s 

+ (22 一 120)]2 = 
Bid, s) KJ Newton 图 画 在 图 28 上 ,由 它 可 知 o 一 1, 从 而 
近似 为 


7.0| -60 -50 -40 -30 20 -l 0 50 60 70 


图 30 完全 特征 频率 轨 线 
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Fi e 算出 是 
912 — 4320 
88 + i104 
最 后 特征 频率 轨 线 在 极点 一 4 十 2; 周围 具有 下 述 近 似 : 
(一 912 + 1320) 
(88 + 4104) 
因而 由 * = 一 4 + 27 这 一 极点 的 出 发 角 是 


一 912 + 7320) L 1109。 
ma ( 88 + j104 ) 


ERA s = 一 4 一 2j 处 
应 用 对 称 性 ,由 极点 * 一 一 4 一 27 URE —110.9°, 
Gii 求 抵达 角 
系统 具有 两 个 有 限 零点 
5 一 1 十 1 5 s=1-7 
在 零点 一 1 十 i 处 
在 特征 方程 中 令 * =s—-1-7, WE 
XCgos’) = [s'* +19 + 74)? + (13 十 129)57 
+ (一 58 十 146)5 十 (一 18 — 738) 1g 
+ [59 + (3 + 73)s7 + (116 十 18)7 
十 (一 318 + {118)]g — 121s” + j2s'] 


之 一 4 十 2i 十 


=0 
X(g, s’) 的 Newton 图 如 图 29 所 示 , 由 此 有 ?一 工 且 近 似 为 
s = pg 
系数 ?经 计算 为 
一 318 + f118 
124 


最 后 特征 频率 轨 线 在 零点 * 一 | 十 i 周围 的 近似 为 
expe y+ GIS) 人 
124 
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图 31 频率 曲面 (a) 第 1 页 ; (b) 第 2 页 


因此 对 零点 “一 1 十 7 的 抵达 角 是 


辐 角 (mme a 180° = 69.64° 
1 


在 零点 * 一 1 一 7 处 

由 对 称 性 得 零点 * 一 工 一 1 的 抵达 角 是 一 69.64° 

为 了 核对 结果 ， 把 完全 的 特征 频率 轨 线 投影 在 单一 复 
平面 上 并 画 在 图 30 上 ， 轨 线 的 分 支 在 极点 :二 1 附近 重合 
造成 对 羽 环 极点 运动 理解 上 的 困难 。 es) 的 常 位 相 与 党 
辑 值 曲线 措 画 了 的 频率 曲面 表示 在 图 31(a)，(b) 上 ， 以 给 出 
特征 频率 轨 线 的 清楚 的 表述 ， 极 点 :二 一 4 十 2i 与 零点 
5 二 1 + 1 处 领域 内 的 小 区 域 画 在 图 32 与 33 上 ， 以 验证 计 
算 的 出 发 角 与 抵达 角 。 


图 32 由 极点 :一 一 人 十 27 图 33 到 零点 ?一 1 ty 
的 出 发 角 的 抵达 角 


65 最 优 闭环 极点 的 渐 近 性 质 


在 这 一 节 “Newton 图 "方法 被 用 来 确定 多 变量 时 不 变 线 
性 调节 器 的 最 优 闭环 极点 在 品质 指标 中 控制 项 权 趋 于 零 时 的 
渐 近 性 质 。 其 方法 是 基于 最 优 特征 频率 轨 线 与 合适 的 代数 函 
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数 分 支 之 间 的 联系 的 。 虽然 这 一 方法 等 价 于 Kwakernank!111 
的 办 法 ， 但 在 应 用 代数 函数 理论 以 后 该 方法 的 简单 性 就 得 以 


突出 了 . 
考虑 可 镇 定 又 可 检测 的 时 不 变 线性 系统 : 
a = Axle) + Bule) (6.5.1) 
t 
y) = Cale) (6.5.2) 
与 品质 指标 


Veo) 一 人 LOVE) + our(DORuCD1a (6.5.3) 

其 中 2 与 R 为 正定 对 称 和 矩阵 .符号 工 表示 和 矩阵 或 向 量 的 转 置 . 
于 是 熟知 的 最 优 控 制作 用 是 

u(t) 一 一 1 R™'BTPx(2) (6.5.4) 


Hh p AMARE Riccati 方程 
—PA — ATP + J PBR™B'P =CTOC (6.5.5) 


的 唯一 半 正 定 解 。 
Kwakernaak 【11, 方 程 (16)] 求 得 最 优 调节 器 下 闭环 特征 
多 项 式 (用 DAs) 表示 ) 与 开 环 特征 多 项 式 (用 Dols) ER) 
之 间 的 关系 是 
OAs) SBA—s) 
= Dols) gw 一 9da| In 十 二 R-GT 一 0)9GC9)] 


(6.5.6) 

其 中 In 为 m 阶 单位 矩阵 ,m 为 系统 输入 个 数 , 而 且 
G(s) = C(s1 — A)-B (6.5.7) 
是 系统 的 开 环 传递 函数 或 增益 矩阵 。 由 式 (6.5.6) 明显 看 出 ， 
最 优 闭环 系统 的 极点 依赖 于 输入 权 ， 其 在 左 半 平 面 的 * 值 应 
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WE ; i 
det | In + HG) =0 (6.5.8) 
其 中 
H(s)8R7'G61(—s)0G(s) (6.5.9) 
现 考虑 确定 As) 特征 值 的 特征 方程 , 它 是 
ACn, s) Sdet (nl, — H(s)] = 0 (6.5.10) 
或 用 展开 行列 式 得 
Ans s) = n” + a(i H e 
+ anl) = 0 (6.5.11) 
其 中 系数 (aC): = 1, 2,，…, m) E WARAN. 系数 
是 ?的 函数 的 理由 是 
A(n, s) = detlnln 一 R“G™—s)OG(s)] 
= det[n1n 一 R7#G™(—s)OG(s)R™*] (6.5.12) 
MA RAG —s)OG(s)R-* 是 平行 Hermite 和 矩阵 就 能 推 知 
Aln, =A, 一 s)[11, 附 录 A]， 而 这 必然 发 生 在 系数 是 
ARRA. 
EO 是 系数 (al) = 1,2,…,m} 的 最 小 公分 
母 , 则 日 式 (6.5.11) 有 
bl)ACn, 5) = balsy” + bl + + 
+ bnl) = 0 (6.5.13) 
其 中 系数 (b)i = 1,2, --+, mp 是 的 多 项 式 。 由 式 
(6.5.13) SCAU EAR o(s) 是 一 代数 函数 。 一 般 讲 ， 式 
《6.5.11) 与 (6.5.13) 可 以 分 解 成 在 * 的 有 理 函 数 域 上 的 几 个 不 
可 约 方程 ,因此 定义 了 一 代数 函数 集合 。 但 是 , 为 了 表达 简 
单 ， 将 假设 式 (6.5.11) 与 (6.5.13) 均 为 在 s 的 有 理 函 数 域 上 不 
可 约 . 
如 果 我 们 比较 式 (6.5.8) 与 (6.5.10), 则 最 优 闭环 极点 能 由 
方程 


2972 6 


n(s) = —e (6.5.14) 

的 左 半 平 面 解 确定 ， 即 最 优 特征 频率 轨 线 是 代数 函数 在 左 半 

平面 的 180° 等 位 相 曲 线 . 若 将 ? 考虑 为 复 变量 且 在 式 (6.5.13) 
中 替换 m(s), WA 

b(s?)(—p)” + b,(s*)(—p)™7 tS easy 

+ dy (s?) = 0 (6.5.15) 

它 是 一 个 用 来 确定 代数 函数 Co) 与 p(s) 的 代数 方程 .对 代数 

PAR Co), 4o 为 正 实 数 时 其 在 左 半 平面 的 分 支 即 为 最 优 特 

征 频率 轨 线 [注意 方程 (6.5.15) 也 可 由 方程 (6.5.8) 直 接 推 得 ]. 


为 了 确定 最 优 特征 频率 轨 线 的 渐 近 性 质 ， 需 要 给 出 代数 
PAIK CoE s ==co 的 邻 域内 当 。 BF SMH KAO, 
可 作为 对 (o) 的 分 支 在 点 。 一 co 周围 当 趋 于 零 时 的 级 数 
展开 的 首 项 来 求 得 .为 此 目的 在 方程 (6.5.15) 中 令 * 一 at 并 
HER 6.2 节 中 所 做 的 工作 . 


66 例 2 


为 了 阐明 上 述 过 程 ,对 
E 1 
G= l +1 F2 + 3X5 +4) 
x[s+2 6 
l 3 十 3 1 | 
与 
eR =pl 
计算 最 优 闭环 极点 的 渐 近 性 质 ， 从 这 些 数 据 我 们 有 
H(s)AR™G™(—s)0G(s) 


本 1 
(è — 1)(? — 4)(2 — 9) — 16) 
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ae tB —7s +15 i 


7s +15 37 
因而 
Aln, s) Bdet([n 7, — H(s)] 
ge (一 22 + 50)n 
(a — 1)? — 4)(3 — 9G? — 16) 
(一 2592 + 256) 


ta — 1)(3 — 4)? — 9)(?— 16)? 
一 0 
EREN n= 一 p 且 乘 其 系数 的 最 小 公分 母 则 有 
(st — 1) 5? — 4s? — 9)*(s? — 16)p’ 
+ (è — 1)(8? — 4)(s? — 9)(s? — 16) 
x (一 22 十 50)p 十 (一 252 + 256) = 0 
Hs = xz: 替换, 则 有 
a= zya 一 422)2(1 一 9z2)2(1 — 1627)'p? 
十 zi(1 — 2?)(1 — 427)(1 — 927) 
x (1 — 1627)(—2 + 5027p 
+ 2(—25 + 2562?) = 
这 最 后 一 个 方程 的 Newton 图 表示 在 图 34 上 ,由 此 有 


， 和 一 二 ,因而 近似 为 > Satza ept, 


e. 
aes 
6 


HR a Zils e 一 ae RA, 

即 对 应 Newton 图 上 在 第 一 个 
线段 上 的 点 的 项 等 于 零 

就 可 得 e 一 /05 0.891exp PiE; R=0,1,2, .5 


Pp'—2zp= of 


9) 此 处 分 母 原 书 为 (? 一 IF — 2N 一 3)" 一 4) 系 错 的 ?随后 的 式 子 
西 由 此 导出 亦 错 * 但 这 些 并 不 影响 Newton ARMs. 一 一 译 者 注 
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Ce een aad 
-20 20 4C Rels) 


一 -一 一 


+ 
pee 20 40 Re{s) 


(b) 


图 于 Newton 图 图 35 WE (5) 的 Riemann 
曲面 上 的 最 优 特 征 频率 轨 线 (加 上 
右 半 平面 的 象 ) 的 渐 近 性质 
(2) 曲面 的 第 工 页 ; 
(b) 曲面 的 第 2 页 


若 用 > 一 ett A, 
; 即 对 应 Newton 图 上 在 第 二 个 

—252" -2r = of) : op: 

RELNANMSTS 


就 可 得 


ea = 9 二 008 =0.729exp IE 二 2 


k=0,1,2,.,7 
如 果 再 替换 回 * 且 仅 取 左 半 平 面 解 ， 就 可 得 到 下 述 对 最 优 特 
征 频 率 轨 线 渐 近 性 质 的 表达 式 : 


$1.12 [ee(i 1) p; k=2,3,4 


5 
\] -4 
5 1.37 [exp (ste) p *3k=2, 3,4,5 


注意 应 用 Kwokernaak 的 表示 法 [11], 上述 这 些 就 确定 
了 3 阶 与 4 阶 的 Butterworth 模型 。 最 优 特征 频率 轨 线 和 它 
们 的 右 半 平面 的 象 表示 在 图 35 中 ,注意 它们 已 作为 代数 函 
A aC) HY 180° 等 位 相 线 描画 在 它 的 Riemann 曲面 区 域 上 , 
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第 七 章 ”参数 稳定 性 与 进一步 的 研究 


如 果 一 反馈 系统 的 全 部 闭环 极点 都 位 于 左 半 平 面 ， 则 该 
系统 就 称 为 稳定 。 控 制 系统 的 稳定 性 是 依赖 于 其 所 联系 的 一 
些 参数 的 。 有 时 在 控制 系统 中 由 于 可 能 出 现 的 老化 、 损 坏 与 
变质 等 使 参数 的 值 是 不 确定 的 ,在 另 些 时 候 , 为 了 经 济 的 理由 
又 可 能 希望 去 改变 参数 值 。 在 这 样 两 种 情形 下 预测 关于 给 定 
参数 的 系统 的 相对 稳定 性 是 很 有 用 的 

在 前 几 章 一 个 占 优势 的 论题 是 系统 与 两 个 代数 函数 集 
合 一 一 特征 增益 函数 与 特征 频率 函数 的 联系 在 这 一 章 引 人 
了 特征 参数 函数 并 且 将 它 用 来 发 展 参数 根 轨迹 与 参数 Nyquist 
轨 线 的 思想 。 由 此 系统 关于 单个 参数 的 相对 稳定 性 就 可 得 以 
确定 。 

本 章 的 最 后 一 节 , 是 为 进一步 的 研究 提出 的 意见 和 建议 。 


7.1 特征 频率 与 特征 参数 函数 


考虑 的 反馈 结构 如 图 36 所 示 ， 其 中 Alkas kss ts ka)s 
BCA2 kss tts ka), CCkas kss ot ka) 与 Dk Ray ka) 
是 状态 空间 和 矩阵， 它们 依赖 于 (4 一 1) 个 实 的 时 不 变 参数 而 
是 标量 时 不 变 增 益 参 数 为 整个 回路 所 共有 . 

这 个 结构 的 闭环 极点 是 方程 


det[s1 — S(k)] = 0 (7.1.1) 


的 解 , 其 中 
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图 36 为 进行 参数 分 析 的 反馈 结构 


SCR)EA lkiya kul 一 Blais -+*, kal 

XRT m + DR, `>, Ra 一 CC + +9 ka) 
是 闭环 频率 矩阵 (参见 第 3.1 节 )， 如 果 除 外 的 所 有 参数 的 
数值 均 已 代入 方程 (7.1.1), 而 将 k 考虑 成 一 复 变量 ， 则 最 后 
的 代数 方程 (为 简单 起 见 将 认为 它 是 不 可 约 的 ) 将 确定 一 对 代 
数 函数 (Ai) 与 RC). 代数 函数 Chi) 称 为 关于 的 特征 频 
率 函 数 ,而 代数 函数 kC) 为 对 k 的 特征 参数 函数 .( 注 意 ,在 
第 三 章 引 进 的 特征 频率 函数 (e) 与 特征 增益 函数 g(s)， 它 
们 对 应 地 等 价 于 (K) 和 一 As) 

SCR) 的 分 支 在 取 实 值 时 清楚 地 决定 闭环 极点 相对 ki 
的 变化 ,因此 就 称 为 参数 根 轨 迹 。 换言之 ,参数 根 轨迹 可 看 作 
是 在 ki(s) 的 频率 曲面 区 域 上 《i(s) 的 0° 等 位 相 曲 线 . 

与 参数 根 轨 线 相对 偶 的 是 参数 Nyquist 轨 线 或 特征 S 
RAR, TEAS) 的 分 支 在 * 走 过 虚 轴 描画 的 轨迹 ， 换 言 
之 ,特征 参数 轨 线 可 看 作 是 在 (ki) 的 Riemann 曲面 定义 域 
上 xf) 的 土 90° 等 位 相 回 线 , 该 Riemann 曲面 称 为 对 心 的 
参数 曲面 . 

如 朱 对 一 特定 系统 ,对 于 系统 的 参数 有 一 名 义 值 的 集合 ， 
则 我 们 可 以 用 察看 参数 曲面 的 办 法 来 确定 任何 相对 稳定 性 是 
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RESA GRAE, SUE Tas SEE), 
7.2 ”增益 与 位 相 容 度 


在 心 的 参数 曲面 上 *(f) 的 90° 等 位 相 回 线 示 出 了 稳 
定 的 与 不 稳定 的 闭环 极点 间 的 边界 ， 由 此 我 们 就 能 确定 一 个 
稳定 工作 点 名 对 的 参数 增益 与 位 相 裕 度 , 它 给 出 系统 关于 
的 相对 稳定 性 的 一 种 测量 . 

参数 增益 容 度 ”参数 增益 次 度 是 相对 于 稳定 工作 点 ki, 
定义 为 使 系统 成 为 不 稳定 所 要 求 的 在 参数 增益 上 的 最 小 改 
变 ， 令 di 是 沿 着 实 轴 由 一 稳定 工作 点 名 到 的 第 i 页 参数 
曲面 上 稳定 性 边界 (特征 参数 轨 线 ) 的 最 小 距离 . 于 是 参数 增 
益 裕 度 定义 为 


min{d;?i = 1,2, -++,n} 


参数 位 相 容 度 Zk 的 ”页 参数 曲面 的 每 一 页 上 ,想象 
作 一 弧 线 , 其 圆心 在 原点 , 起 点 在 稳定 工作 点 名 , 另 一 端 达到 
稳定 性 边界 (特征 参数 轨 线 ). 令 py 是 在 原点 上 述 第 i 页 上 弧 
对 应 的 中 心 角 , 则 参数 位 相 裕 度 定义 为 


min {8:7 = 1, 2, -sn} 
7.3 例 


在 这 节 考 虑 一 倒 摆 位 置 系统 ( 见 图 37), 并 分 析 相 对 于 其 
单 参数 车 质量 的 系统 的 稳定 性 . 

这 个 系统 同样 被 Kwakernak 与 Sivan [2]，Canaon 13] 
与 Elgerd[4] 利用 过 .系统 可 采用 下 述 线性 状态 微分 方程 模 
All: 


© 105+ 


0 1 
-i 
ile) = 0 0 
-£ 0 
L 


0 0 


0 0 


2 + | Ia) (7.3.1) 
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其 中 ue) 是 小 马达 在 车 上 产生 的 作用 力 ，M 为 车 子 质量 ，F 


i 
枢 轴 
车 
图 37 DAERA 
可 知 [2] 


为 联系 到 车 运动 的 摩擦 系数 ; 而 L 
给 定 为 


,+mnL’ 

Litn (732) 
Ht mAH, LRA 
至 摆 的 重心 间 的 上 距离 ， 而 了 是 摆 相 
对 重心 的 惯性 和 拢 . 

应 用 形 如 

u(t) =—Kx(t) (7.3.3) 

的 反馈 系统 是 能 稳定 的 。 若 取 数值 


(7.3.4) 


' = 0.842m 


K = [86.81, 12.21, —118.4, —33.44] (7.3.5) 
可 镇 定 线性 化 系统 且 使 闭环 极点 为 一 4.706 土 j1.382 与 


— 1.902 + j3.420, 
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39 第 2 页 参数 (质量 ) 曲 面 


-two 


40 第 3 页 参数 ( 质 量 ) 曲 面 
Tm(M) 


图 41 第 4 页 参数 (质量 ) 曲 面 


现在 来 研究 M 的 参数 曲面 看 看 小 车 质量 在 稳定 工作 点 
lkg 周围 其 质量 的 改变 是 如 何 影响 系统 稳定 性 的 。 四 页 质量 
曲面 (如 图 38 一 41) 上 画 好 了 (M) 的 常 位 相 与 常 辐 值 曲线 ， 
由 此 可 以 得 到 下 述 稳定 性 裕 度 : 

参数 (质量 ) 增 益 裕 度 一 lkg 

参数 (质量 ) 位 相 裕 度 = 60° 
lkg 的 增益 裕 度 对 应 车 质量 减少 至 零 仍 保持 稳定 性 ， 参 数 曲 
面 同样 表明 对 稳定 性 说 来 车 的 最 大 质量 为 2.125kg，60° 的 位 
相 裕 度 指 出 闭环 系统 有 充分 的 阻尼 , 


7.4 进一步 的 研究 


可 以 想到 一 个 系统 在 其 正常 工作 时 有 一 个 参数 承受 大 变 
化 的 情况 下 ， 在 为 该 系统 而 用 的 参数 相关 控制 器 的 设计 中 参 
数 曲面 可 能 是 相当 有 用 的 .例如 ,对 飞机 引擎 的 控制 器 说 来 ， 
就 要 求 在 相当 广 的 高 度 范围 内 合适 地 工作 .一 个 可 能 的 设计 
方案 可 以 是 

(i) 在 一 特定 高 度 上 设计 一 实 常数 控制 器 . 

Gi) 用 " 德 阵 插值 "手法 求 得 一 高 度 相 关 控制 器 , 

Gii) 应 用 ”高度 曲面 ”分 析 在 整个 工作 范围 内 系统 的 稳定 
HE. 

如 同 在 本 章 所 指出 的 ， 在 这 一 书 中 所 提出 的 方法 不 仅 可 
应 用 于 增益 与 顷 率 , 而 且 可 用 于 任何 单一 的 系统 参数 与 频率 ， 
要 获得 多 于 一 个 参数 改变 下 的 稳定 性 判 据 是 一 个 复杂 的 问 
题 ,但 却 有 很 大 的 实际 意义 ,可 以 感到 从 多 复 变 函 数论 的 研究 
中 可 能 会 有 有 价值 的 看 法 . 

在 最 近 几 年 里 对 二 维 及 多 维 数值 滤波 器 中 稳定 性 检验 得 
到 了 发 展 [5], 二 维 数值 滤波 器 正 广泛 地 应 用 于 各 领域 ; 例如 

+ 109+ 


求 象 过 程 \ 地 震 过 程 的 地 球 物理 学 ,重力 与 磁性 资料 。 可 以 期 
望 高 维 滤波 器 同样 可 以 找到 应 用 的 领域 :例如 全 息 三 维 滤波 ， 
由 用 Nyquist 方式 表现 出 来 的 稳定 性 判 泥 可 以 想到 这 些 发 展 
经 过 控制 工程 师 的 更 深刻 的 理解 可 以 为 在 控制 领域 内 的 应 用 
得 到 合适 的 改进 ， 例 如 可 能 在 多 参数 变化 下 稳定 性 的 研究 中 
得 到 应 用 。 

贯 串 这 一 工作 的 经 常 性 的 论题 是 多 变量 系统 与 一 个 或 可 
能 多 个 定义 在 合适 的 Rieman 曲面 上 代数 函数 之 间 的 联系 。 
对 一 代数 函数 而 言 Riemann 曲面 是 拓扑 等 价 于 具 环 柄 的 球 ， 
而 环 柄 的 数目 是 当 作 曲面 的 亏 格 数 来 理解 的 [6]。 污 格 数 可 
以 证 明 是 系统 的 一 个 重要 特征 , 它 可 能 与 解 簿 有 关 , 即 变换 多 
变量 系统 到 一 个 系统 ， 实 际 上 它 是 一 单 输入 单 输出 系统 的 集 
A. 一 单 输 入 单 输出 系统 有 一 对 应 零 亏 格 的 显然 (一 页 ) 频 率 
曲面 ， 所 以 解 耦 多 变量 系统 将 等 价 于 化 简 一 个 具 亏 格 数 大 于 
MEER Ym Riemann 曲面 到 六 个 显然 Riemann 曲面 ( 具 
有 亏 格 数 零 ) 的 集合 ， 于 是 在 分 支点 奇异 性 ( 它 的 性 质 和 位 置 
确定 如 同 亏 格 数 一 样 的 拓扑 性 质 ) 与 解 耦 之 间 的 关系 将 可 能 
是 一 成 效 显 著 的 研究 路 线 。 在 附录 7 上 发 展 了 一 个 有 兴趣 的 
关系 ， 它 表明 在 增益 曲面 上 的 一 个 分 支点 对 应 一 个 在 对 应 频 
率 曲面 上 的 分 支点 或 驻 定点 ,而 反 过 来 也 对 . 
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附 录 
1 代数 函数 的 定义 


令 4(9q,y) 是 9 的 形 如 

4(4:2) = fog” H Ao) H o H mC) CALL) 
的 多 项 式 , 其 中 每 个 系数 (f): i= 0,1, +++, m} 本 身 是 v 
的 复 系数 多 项 式 ， 则 一 定义 在 复 "平面 上 的 函数 4(v), 它 满 
足 形 如 

A(q,v) 一 0 (Al.2) 

的 方程 ,9(v) 就 是 一 代数 淫 数 ， 多 项 式 Al, v) 也 可 写成 具 
有 系数 为 4 的 多 项 式 的 的 多 项 式 . 车 按 此 考虑 , 式 (A1.2) 
也 定义 一 代数 函数 "(9)。 

对 于 ”的 一 确定 值 mw， 方程 (A 1.2) 有 mw 个 解 称 为 a) 
的 分 支 , 且 在 w 的 邻 域内 分 支 可 用 宕 级 数 展 开 表 示 [1]. 

ELBE MBI ACG, v) 是 在 (9, v) 上 的 不 可 约 多 
项 式 , 即 A(4,v》 不 是 两 个 或 多 个 (9，v) 的 多 项 式 的 乘积 .， 
若 A(9,v) 是 :个 (49,v) 多 项 式 的 乘积 , 则 方程 (A1.2) 将 定 
X :个 形 如 qlv) 的 代数 函数 ， 或 定义 :个 形 如 v(9) 的 代数 
mm 


2 化 简 至 不 可 约 有 理 正则 型 


在 本 附录 中 给 出 一 个 方法 以 化 简 m% 阶 方 阵 至 其 有 理 正则 
型 ， 给 出 的 方法 与 Ayres 从 方 阵 求 有 理 正 则 型 的 办 法 略 有 区 
别 ， 方 阵 G 的 有 理 正则 型, 除 其 对 角 线 块 对 应 81-G 的 不 变 
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因子 而 不 是 不 可 约 因子 外 将 与 不 可 约 型 相似 。 由 Ayres 给 
出 的 为 寻求 有 理 正 则 型 的 办 法 以 及 某 些 必 要 的 定义 将 在 下 面 
给 出 ， 


定义 AM 
X, GX, @X, +++, GX (A2.1) 
是 线性 无 关 的 ,但 ， 
X, GX, GX, +++, GX, GX (A2.2) 


不 是 , 则 式 (A2.1) 称 为 具有 XX 为 首 ,长 度 为 :的 一 个 链 . 

方法 ”对 给 定 任 一 域 上 的 mw 阶 方 阵 G 

(i) 令 Xm 是 对 下 上 所 有 w 维 向 量 说 来 ， 具 有 最 大 长 度 
REC, 的 首 ; 

Gi) 令 X。-: 是 具有 最 大 长 度 链 C。-: 的 首 〈 该 链 的 元 与 
前 述 Cm 的 元 线性 独立 ), Xm 是 在 与 Cm 的 向 量 线 性 独立 的 
所 有 F 上 ,m 维 向 量 中 选取 的 ; 

Gi) 令 Xm 是 最 大 长 度 链 Cm-; 的 首 (该 链 的 任何 元 均 
线性 独立 于 前 述 Cm 与 Cw- 的 元 ), Xm—2 是 在 与 Cn 和 Cm-1 
的 向 量 线性 独立 的 所 有 上, m 维 向 量 中 选取 的 ; 

如 此 等 等 ， 于 是 对 
E = [Xi, GXi, +++) GI IXi; Xins GXins “°° 
CGOHTIXi 2093 Xms GXms ***» G'm Xm] 
RA EGE 是 G 的 有 理 正则 型. 

在 这 一 方法 中 ， 最 大 长 度 链 是 被 用 来 寻求 不 变 因 子 的 。 
现在 不 变 因 子 已 经 由 不 可 约 有 理 正 则 型 所 要 求 的 不 可 约 特征 
方程 的 科 积 中 产生 。 因 此 如 果 寻 求 最 小 长 度 链 代替 最 大 长 度 
链 , 这 样 形成 的 变换 矩阵 E 就 可 以 给 出 不 可 约 有 理 正则 型 。 

上 述 方法 的 问题 在 于 没 法 指出 给 定 什么 就 能 确 知 什么 时 
修一 个 链 有 最 大 或 最 小 的 长 度 。 此 外 在 相反 情况 下 将 出 现 较 
长 与 较 短 的 链 。 
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指出 最 大 或 最 小 长 度 的 链 刚 好 组 成 G 不 变 子 空间 中 最 大 
或 最 小 子 空间 的 基 这 一 点 是 有 益 的 。 因 此 找 最 小 长 度 链 就 等 
价 于 找 最 小 的 不 变 子 空间 。 同 样 可 以 发 生 的 是 刚好 一 维 不 变 
子 空间 (特征 向 量 ) 对 应 一 特征 值 以 及 : 维 不 变 子 空间 (从 具 
最 小 长 度 中 寻求 ) 对 应 一 确定 + 个 特征 值 的 :次 不 可 约 方程 ， 
例 求 


(十 3) (s+2) 


CGH1IP (+17 

GG) = (s — 3) 一 2 0 
s G+iy GFIY 

-1 -—G+2) 1 


(stl)? 2G+1) s+1 
PARTERNE: & X=[001], 其 中 了 表示 转 置 ，CX = 
ETT 5 一 一 ~ XX, 因而 X 是 一 具有 最 小 长 度 的 链 . 
&Y =[010)7, a 


-| 5 十 2 一 2 . zea 7 
G+1 (6+1 2541 


与 


ey = [+2 $72) —G+2)]' 

: G+ Cty Us + 1) 

i Gx Ty 

由 于 不 能 再 找到 较 小 长 度 的 链 ,因此 Y, GY 构成 一 最 小 长 度 

链 的 集 。 从 而 变换 矩阵 EC) 可 取 为 

Pe +2 

(s+1y 

—2 

GIy 

1 0 一 (人 十 2) 
26s 十 ] 


E(s)={X¥YGY]=j0 1 
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s+ 
0() = ENOCE =| 0 0 GF 
ae 
(s +1) 
由 此 推出 不 可 约 特征 方程 是 
i 
Ales) = 8 — 1 
Sige pa 
Aig, 5) =g a! Gee 


注意 由 在 寻求 最 小 长 度 链 时 得 到 的 相关 关系 。 同样 是 明显 
的 。 


3 判别 式 


在 本 附录 中 ,对 给 定 的 方程 
Cg, s) = bols)g’ 十 六 (8 后 十 二 2) 一 0 
给 出 了 两 个 寻求 其 判别 式 的 方法 . 
方法 1 (Barnett [3]) 
给 定 两 个 多 项 式 : 
a(g) = aog” + aig” * + +++ + an 
elg) = cog” 十 cg +--+ +C,, 


其 中 a;, c; E C, 则 该 丙 个 多 项 式 的 结 式 Rla(g), clg)] 是 下 
述 行列 式 : 
.115 。 


0 a a an-ı Gn 
. Gn; On 
Rlalg), c(g)] = Co ca Cmi Cm 
Co cl C cn 0 
Co Cy G ttt 


多 项 式 o(g) 5 (8) 具有 公 因 子 ( 次 数 大 于 零 )， 当 且 仅 当 
(a + m) 阶 行列 式 为 零 ,其 中 已 设 a, co 不 同时 为 零 . 
设 用 a'《g) 表示 多 项 式 ale) 对 8 的 微 商 , 则 多 项 式 ale 
的 判别 式 是 行列 式 De(ao, ars ***, an)» CH 
De(aoy ars ***, an) = Rla(g), a'(g)] 
KEX. 多项式 a(g) 具有 重 因子 , 当 且 仅 当 判别 式 Delas 
a, "tran 为 零 ， 现 考虑 形式 为 
P(g, 5) = bils)gt + hls)e + 
+b), :>0 
的 多 项 式 ， 其 中 系数 {b)i = 0,1, 2, +++, 2} 均 为 的 
多 项 式 。 这 个 作为 & 的 多 项 式 的 判别 式 可 以 用 前 面 确 定 的 
De (boy bis + 8% 5 bi) L bls), eero bls) 对 应 地 取代 bos 
esb RRE 于 是 存在 一 函数 Ds), 它 也 称 为 判别 式 且 
由 下 式 定义 : 
Dil) = FF RIRE, 8), OC, o) 
其 中 @'(g, s) E Dle, s) 对 8 的 微 商 . 
方法 2 (Sansone 与 Gerretsen [4]) 
考虑 多 项 式 
a(g) = ag? + ag? Heret any n>O 


+1166 


其 中 weC: Male) 的 判别 式 可 按 表达 式 


De(oo。……，an) = aP 


sai Hh PETAR 
ro a Gni 
P=|" a On 


Oni On *** Om J 
而 行列 式 的 元 (oj? = 1,2, +++, 2n 一 2} 是 系数 {qi:i = 0, 
1,2,-+-, n} 的 函数 , 且 
ao 一 力 
TEs ty Oo 满足 下 述 方程 : 
, a; + ao, = 0 
2a: + ayo, + ao, = 0 

(n 一 1)ao + an- 十 十 ac 一 0 
而 ony Ons? ona 可 由 下 述 方程 求 得 : 

Ano 十 an_i0 十 … .十 an 一 0 


OnT, 十 an0 十 … + anny, 一 0 


现 考虑 多 项 式 
P(g,s) = bse + blg + +++ + bs) 
这 一 作为 8 的 多 项 式 的 判别 式 可 以 用 前 面 定义 的 Ds(&,*…， 
b) 通过 用 bls), oes D) 对 应 地 取代 bo, ++, be 来 得 
到 *)， 多 项 式 Og, s) 具有 重 因子 , 当 且 仅 当 
D,(s) = 好 (5)P 


P 按 方法 2， 此 处 boy …， bi 实际 应 为 6s aiy ty ae, 译 考 注 
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HE, Hoh P 为 一 矩阵 的 行列 式 。 其 元 是 系数 (0): = 0, 
1,2，.…， 分 的 函数 


4 一 个 构造 对 应 开 环 增益 矩阵 G(s) 的 代数 
函数 的 Riemann 曲面 定义 域 的 方法 


在 第 3.3-4 节 已 经 阐明 对 于 + 次 的 特征 增益 函数 ,对 应 
的 Riemann 曲面 是 由 * 页 复 * 平面 沿 在 一 些 分 支点 与 无 穷 远 
点 之 间 的 割 口 相互 粘 在 一 起 而 构成 的 ， 虽 然 在 某 种 意义 之 下 
割 口 是 任意 的 ( 即 不 存在 唯一 的 割 口 集合 ), 但 去 选择 一 个 相 
容 的 集合 仍然 是 一 个 问题 ， 这 里 相 容 是 指 能 辨识 以 什么 次 序 
可 以 将 这 些 页 联 在 一 起 。 为 了 解决 这 一 问题 本 附录 给 出 一 个 
系统 的 方法 ， 这 一 方法 的 优点 在 于 其 作出 的 割 口 除了 沿 实 轴 
的 割 口外 其 余 的 割 口 均 以 实 轴 为 对 称 且 经 常平 行 于 虚 轴 , 
同样 ， 由 于 这 一 方法 直接 应 用 于 开 环 增益 矩阵 从 而 不 需要 去 
求 特 征 方程 ， 并 且 最 后 由 相互 连接 的 页 组 成 的 彼此 不 关联 的 
集合 表示 的 是 不 可 约 特 征 方程 的 Riemann 曲面 。 

方法 的 基础 在 于 对 覆盖 * 平面 的 一 个 网 格 上 的 值 寻求 对 
应 矩阵 的 特征 值 ， 然 后 将 求 得 的 特征 值 沿 网 格 的 线 整理 成 连 
续 的 形式 ， 落 传递 函数 是 mx m 阶 的 ， 将 可 以 有 效 地 表示 
Riemana 曲面 页 数 的 mw 对 应 的 ;数组 用 来 存放 特征 值 ， 其 过 
程 与 第 3.3-4 节 中 描述 的 解析 开拓 过 程 类 似 ， 那 里 圆 盘 上 的 
个 体 点 可 作为 唯一 函数 值 的 承受 点 ， 这 里 :平面 各 页 (数组 》 
上 的 个 体 点 也 是 作为 函数 值 的 承受 点 . 

第 一 条 沿 着 它 计 算 特 征 值 并 加 以 整理 的 线 是 实 轴 ， 然 后 
的 计算 与 整理 将 沿 着 由 实 轴 上 点 出 发 平行 于 虚 轴 的 线 进 行 ， 
如 图 42 所 示 。 由 于 下 半 平 面 的 特征 值 仅仅 是 上 半 平 面 特征 
ABA, 因此 计算 工作 仅 需 在 上 半 平 面 进行 ， 继续 这 一 
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过 程 , 则 s 平面 被 计算 点 所 覆盖 , SEAN m BSR, 
它们 在 沿 着 实 轴 与 和 虚 轴 平行 的 线 上 是 连续 的 但 末 必 在 与 实 
轴 相 交 的 那些 地 方 连续 . 


图 42 计算 与 整理 特征 值 所 要 沿 的 线 


通过 观测 每 一 数组 和 页 次 则 必要 的 割 口 就 清楚 了 .如 果 
平行 于 虚 轴 的 一 个 特征 值 的 线 是 与 邻近 的 一 个 线 之 间 不 连 
续 ， 则 它们 之 间 必 然 被 一 个 由 一 个 分 支点 出 发 终结 在 另 一 个 
分 支点 或 无 穷 远 点 的 割 口 所 分 离 ;如 图 43(a), (b) 所 示 。 

x 一 一 x 指出 特征 值 之 


闻 的 连续 性 = 


图 43 ”平行 于 碟 轴 的 割 口 
(a) ARA (b) EJN 


对 应 上 半 平 面 : 值 的 特征 值 是 其 在 下 半 平 面 对 应 值 的 揽 
共 圈 ,因此 当 在 实 轴 上 算出 的 特征 值 为 复数 时 ,通过 实 轴 而 要 
求 特征 值 连续 就 是 不 可 能 的 。 因 此 当 实 轴 上 特征 值 为 复数 
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44 沿 实 铀 的 制品 
IN, 一 个 沿 实 轴 的 割 口 就 是 必要 的 了 。 又 所 有 的 割 口 均 是 由 
分 支点 出 发 而 终结 在 分 支点 或 无 穷 远 点 ,如 图 44 所 示 。 

用 对 比 在 割 口 一 边 一 页 上 的 特征 值 与 盾 口 男 一 边 的 特征 
值 的 办 法 ,把 各 页 互相 进行 粘 合 也 就 比较 显然 了 .一 般 讲 , 当 
这 种 对 比 过 程 完成 时 ， 一 些 由 彼此 联通 的 页 组 成 的 集 就 形成 
了 ,每 一 个 集 均 定义 一 个 Rieman Hi. 

使 辨识 割 口 和 对 比 页 次 的 工作 简化 对 于 在 各 页 上 描画 等 
位 相 与 等 辐 值 回 线 的 工作 是 十 分 有 用 的 ， 这 已 由 正文 中 给 出 
的 例子 所 说 明 。 由 于 实 轴 极点 的 似 然 性 , 沿 着 实 轴 计 算 与 整 
理 特征 值 可 能 有 困难 。 当 第 一 根 计 算 与 整理 的 线 变 至 平行 于 
实 轴 且 保持 一 较 大 距离 时 ,上 述 困难 就 可 望 得 到 克服 。 当 沿 
此 线 的 计算 与 整理 的 工作 终止 后 ,进一步 的 计算 与 整理 ,将 沿 
由 这 根 新 线 出 发 终止 在 实 轴 上 而 与 虐 轴 平行 的 线 进行 ， 如 图 


Im ts 外 


AD 计算 与 整理 符 征 值 记 要 褒 的 线 ( 改 进 方法 ) 
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45 所 示 。 利用 这 一 新 方法 可 能 的 割 口 仅仅 是 或 者 在 实 轴 上 
或 者 在 复 共 孝 分 支点 之 冶 , 如 图 46 所 示 。 


图 47 第 1 页 Riemann 曲面 


对 于 


Gs) = 1 ice ee 
1.25¢s + 1)(s +2) [3 — 18 s— 8 


用 原来 的 办 法 所 作 的 Riemann 曲面 的 构造 如 图 47 与 48 所 
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图 49 第 1 页 Riemann 曲面 《改进 的 构造 方法 ) 


— yj 


图 50 第 2 页 Riemann 曲面 (改进 的 构造 方法 ) 
示 ， 图 49 与 50 是 应 用 改进 了 的 办 法 画 的 并 以 此 说 明 割 口 的 
任意 性 。 


5 扩展 了 的 辐 角 原 理 


辐 角 原理 按 我 们 所 要 求 的 扩展 在 文献 中 看 来 还 不 易 找 
到 ,虽然 对 于 一 般 多 值 解析 函数 说 来 ,原理 的 合适 的 提 法 已 经 
在 有 关 书 中 出 现 了 。 [Evgrafov 5,p. 98]。 因此 ,为 在 第 四 、 
五 章 中 运用 恰当 ,这 里 将 给 出 辐 角 原理 一 个 合适 的 扩展 。 


A5.1 引言 


我 们 所 要 求 的 扩展 是 不 明显 的 这 需要 解决 两 个 问题 。 第 
一 个 是 由 这 样 一 个 事实 提出 来 的 ， 即 一 般 讲 代 数 函 数 的 Rie- 
mann 曲面 是 多 连通 的 。 这 里 举 一 个 表明 困难 缘由 的 例 , 考 虑 
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与 2X2 开 环 增益 矩阵 G(s) 相 联 系 的 特征 增益 函数 g(s), 设 
8(s) 具 有 4 个 分 支点 且 每 个 分 支点 均 联系 到 两 页 的 循环 [6]， 
于 是 联系 到 Riemann 曲面 的 亏 格 数 是 116, 7]， 进 而 设 这 些 
支点 按 图 51 的 办 法 设置 在 频率 平面 上 (s PE) 取 两 层 的 


Im fs) 复数 球 (Riemann HER), Æ 
as= 分 支点 

分 支点 间作 割 口 ， 于 是 用 常 

: 见 的 方法 这 个 Riemann 曲面 


Bets) 将 拓扑 等 价 于 一 个 圆 环 面 

[7] ,如 图 52 与 53 所 示 。 对 

在 环 面 上 边界 为 00 的 具有 

阴影 的 区 域 2， 它 对 应 的 是 

图 51 频率 平面 后 来 在 其 上 构造 出 环 面 的 原 

复数 球 上 一 对 Nyquist D 回 线 的 内 部 (如 图 52 所 示 )。 面 对 

这 样 一 种 情形 ,我 们 必须 保证 辐 角 原理 的 扩展 能 对 于 Riemann 

曲面 上 这 样 的 区 域 成 立 ， 邑 它 不 是 单 连通 的 而 且 其 边界 可 以 

由 几 个 互相 分 离 的 闭 Jordan 曲线 所 构成 。 原理 的 这 样 一 个 

基本 的 进展 是 利用 Cauchy 留 数 定理 的 适当 推广 来 完成 的 ,而 
后 者 在 第 A5.2 节 中 论述 ， 


a DRAR 
SMa 


图 52 两 层 复数 球 


在 导出 合适 的 扩展 了 的 辐 角 原理 的 过 程 中 ， 要 克服 的 第 
二 个 障碍 是 直接 与 代数 函数 的 分 支点 相 联 系 的 情况 。 在 这 一 
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情况 下 将 影响 到 用 绕 边 界 689 计算 回路 积分 以 算出 留 数 来 ， 
02 包围 的 区 域 2 以 分 支点 为 自己 的 内 点 。 这 个 问题 可 以 通 
过 复 变量 的 改变 的 办 法 加 以 克服 ,这 将 在 第 A5.3 节 中 提出 。. 


图 53 所 示 的 Riemann 曲面 拓扑 等 价 于 图 环 面 


最 后 在 第 A5.4 节 将 联合 第 A5.2 与 第 A5.3 两 节 的 结果 
给 出 所 要 求 的 扩展 了 的 辐 角 原理 。 


AS.2 Cauchy 留 数 定理 的 广义 形式 


基本 的 要 求 是 将 Cauchy 留 数 定理 推广 至 定义 在 某 些 已 
知 代数 函数 的 Riemann 曲面 上 的 复 变 函 数 上 ,这 里 积分 区 域 
可 以 是 非 单 连通 的 而 且 其 边界 可 以 是 由 一 个 或 几 个 闭 Jordan 
曲线 所 组 成 。 这 一 推广 可 以 在 Bliss [6] 的 著作 中 找到 。 E 
要 定理 将 给 在 下 面 ， 但 为 了 适应 这 本 书 的 要 求 将 在 符号 与 写 
法 上 作 些 变更 ， 为 了 帮助 不 熟悉 Riemann 曲面 的 读者 了 解 
定理 ,指出 下 述 两 点 是 有 用 的 . l 

(i) 在 代数 函数 f(s) 的 Riemann 曲面 上 的 每 个 非 奇异 
位 置 (或 点 ) 唯一 地 被 一 对 值 (f, s) 所 确定 ,这 一 对 值 满足 一 
已 知 代数 方程 


Af, 5) = 0 
因此 在 Riemann 曲面 上 的 每 一 点 , 复 变量 * 与 了 这 一 对 值 的 
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ERM uv, s) 是 已 定义 了 的 . 

Gi) 对 于 一 代数 函数 说 来 Riemann 曲面 是 一 有 向 曲面 ， 
于 是 可 以 在 其 上 定义 一 个 正 向 以 使 在 此 意义 下 一 边界 89 被 
扫 过 就 和 它 围绕 某 区域 一致， 实现 这 一 手段 最 简单 的 方法 
就 是 想象 有 人 在 9 所 在 的 二 维 流 形 上 (GRERE-ABN 
三 维 空间 中 ) 顺 着 边界 00 行走 ,其 绕 行 的 正 向 系 指 在 此 意义 
下 顺 着 60 的 任 一 部 分 前 进 总 使 0 在 行走 者 的 右 方 .在 图 53 
上 边界 69 就 是 按 正 向 画 出 的 . 
广义 留 数 定理 

令 尽 系 一 代数 函数 f(s) 的 Riemann 曲面 ，0 是 尺 的 一 
部 分 其 边界 由 一 个 或 几 个 不 经 过 Riemann 曲面 的 奇异 位 置 的 
Jordan 闭 曲线 所 组 成 , 但 其 本 身 未 必要 求 单 连通 #— R E 
WEERA wf, s) 在 9 与 00 上 除 在 2 内 可 能 出 现 的 有 限 
个 奇 点 外 均 解 折 , 则 有 


[8s Da = Zot, 9) 在 0 内 的 留 数 《A5.2.1) 
2xj Joa 


其 中 边界 00 是 在 相对 9 为 正 的 意义 下 扫 过 的 . 
为 了 推导 扩展 的 辐 角 原理 ， 首 先 将 广义 留 数 定理 应 用 于 
函数 /5)1Ks) ,其 中 fs) 表 As) 对 :的 导数 ,这 样 就 给 出 了 
1 G "(Cs 
ls ra ds = Siow 的 留 数 
注意 fs) 在 边界 99 上 必须 没有 极点 与 分 支点 (它们 二 者 均 
作为 奇异 位 置 )， 为 了 推导 出 辐 角 原理 我 们 也 认为 iC) 在 80 


上 没有 零点 。 下 一 步 就 来 计算 在 0 wie 的 留 数 . 


A5.3” 留 数 的 计算 


f(s) 的 极点 ,零点 与 分 支点 是 ro 的 奇 点 这 一 事实 在 
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本 节 的 讨论 将 是 清楚 的 。 我 们 将 考虑 所 有 Ms) 的 分 支点 , 但 
结果 仅仅 是 那些 又 是 极点 或 零点 的 才 是 重要 的 . 

在 极点 或 零点 (不 包括 那些 与 分 支点 相 联系 的 ) 的 邻 域内 
is) 可 以 用 一 个 确定 的 单 值 函数 的 级 数 表 示 . 对 于 一 个 tv; 阶 
的 极点 Pi 代数 函数 可 以 表示 为 

Hs)= Dp ens —pi)", a, 0 (A5.3.1) 


ntp, 


即 
K) = een (45.3.2) 
其 中 Cp) = 00 H. ACs) 在 ;的 邻 域内 解析 %， 这 就 给 出 
iG) i LO (A5.3.3) 


K) GP) 6) 


BT SO Hp, mannm want Es=p 有 
一 简单 极点 其 留 数 为 一 ty,。 
对 于 re, 阶 零 点 z, 代数 函数 可 用 级 数 
K= D eas — ais cn, 0 (A5.3.4) 


mt, 
ži 


表示 , 即 
16) = Cs — zi) 06) (A5.3.5) 

其 中 0(%) ~ 0 H OCs) 在 s 的 邻 域内 解析 ， 这 就 给 出 
f(D) i 十 82 (A5.3.6) 


K) Gan) 66) 
于 是 有 ns 在 :一 z; 有 一 简单 极点 其 留 数 为 rs。 
s 
为 了 确定 在 分 支点 的 留 数 , 作 变量 的 替换 [6, p.81] 是 必 


V (1) 还 应 不 为 零 . 一 一 译 者 注 
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要 的 ,其 过 程 如 下 : 
在 有 限 分 支点 b 的 邻 域内 代数 函数 可 以 素 示 成 如 下 形 
式 的 级 数 [8]: 


is) = 5 dal Vs — b (A5.3.7) 


ro-o 


其 中 + 是 在 分 支点 组 成 -一 个 循环 系 的 页 数 ， 如 果 分 支点 是 
砍 阶 的 极点 且 组 成 一 个 r 页 的 循环 , 则 可 以 有 级 数 表 达 式 


1:)= DY) a (Vs — bi)", d_a = 0 (A5.3.8) 


"Ë 
"a-i, 


相仿 若 分 支点 是 h 阶 的 零点 , 则 有 
K) = >) ds (Vs— b)’, da 0 (A5.3.9) 


nate 
ži 


否则 
Ks) = 5 dn (Vs — bi)”, do% 0 (A5.3.10) 
n=0 


这 些 表 达 式 是 多 值 的 因而 不 便于 确定 留 数 ， 但 如 果 作 蔡 
换 
sabe (A5.3.11) 
则 可 找到 级 数 表达 式 (A5.3.8, 9 与 10) 对 应 地 为 


te) = >) dex", d, 50 (A5.3.12) 


= 
f(x) = È dax", dy 0 (A5.3.13) 
5 i l 
f(x) = > dax’, dox0 (A5.3.14) 
其 中 f(x) 理解 为 os 
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f(x) = f(2, + x’) (A5.3.15) 
这 个 替换 已 将 代数 函数 映射 到 x 平面 ， 在 这 里 已 可 用 一 
个 单 值 级 数 表达 式 表示 。 
由 定义 [6] 在 b; HBRA CH b; 自身 可 能 外 ) 解 折 的 六 
在 5 的 留 数 是 


1 f@,, 
(BBO, a f ED a (A5.3.16) 


其 中 C 是 在 RR 上 的 一 闭 正 向 Jordan 曲线 , SRE b 的 邻 域 


NN, 在 内 1 除 可 能 bs 一 点 外 均 解析 。 若 对 。 用 替换 式 
(A5.3.11), WA 
wana | E a 
= +f £@) dz ds 
2xj Jez f(x) ds 
-1| 7) 5 
pr k LO as (A5.3.17) 


其 中 C: 是 在 * 平 面 上 界定 原点 一 邻 域 的 闭 正 向 Jordan HR. 


于 是 人 人 在 分 支点 4 的 留 数 就 是 工科 在 原点 的 留 数 ， 而 
1) 1 


TE 的 留 数 是 容易 由 其 级 数 表达 式 得 到 的 。 MATAN 


方法 我 们 注意 到 在 原点 的 邻 域内 oe 


O 若 分 支点 是 za 阶 的 极点 , 则 Sa 有 一 个 具有 留 数 
—t}, 的 简单 极点 ;或 

Gi) 若 分 支点 为 t MEEA, M 有 一 个 具有 留 数 
15, 的 简单 极点 ; 
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fx) 
au ie) 是 解析 的 。 


AS4 扩展 了 的 辐 角 原理 
将 第 A5.3 节 的 结果 结合 到 方程 (A5.2.2), WA 
+f a= (Dat D4 ,| 


2nj Joa f(s) 
一 [B+ >| 
一 4Z —*p (A5.4.1) 
其 中 #Z 与 #P 对 应 的 是 1G) 在 区 域 9 内 的 总 零点 数 与 总 
MAM, 2 是 用 边界 00 按 正 向 围绕 的 在 代数 函数 C) 的 
Riemann 曲面 上 的 区 域 。 注 意 多 重 极 点 与 零点 必须 按 其 阶 次 
计算 个 数 .。 
方程 (A5.4.1) 的 左 端 等 价 于 曲线 集 T 围绕 了 平面 上 原 
点 顺 时 针 绕 图 数 的 净 和 ， 而 TT 是 OO 在 f(s) 映射 下 的 象 ， 表 
示 这 一 顺 时 针 绕 数 的 净 和 为 N(T，0)， 则 最 终 可 得 所 要 求 的 
扩展 了 的 辐 角 原理 的 形式 : 
N(T, 0) 一 #Z — *P (A5.4.2) 


6 对 于 特征 方程 A(g8,s)=0 的 多 变量 枢 点 


考虑 特征 方程 
A(g,s) =0 (A6.1) 
设 我 们 已 求 得 特征 频率 轨 线 的 一 个 分 支 在 《一 c 时 一 o 
周围 的 一 个 近似 的 形式 是 
s = bk (A6.2) 
若 p 是 对 应 该 浙 近 线 的 枢 点 , 则 对 《二 co 有 
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s=p + bk (46.3) 


或 
k= (=Y 46.4) 
FE 
闭环 极点 依赖 于 不 的 关系 由 式 (6.1.5) 给 出 , 即 
T(k, 5) =0 (A6.5) 


因此 在 方程 (A6.5) PRA Ms, MAS s Se 间 的 关系 或 
4 与 间 的 关系 。 由 于 在 方程 (A6.5) 中 s 的 次 数 通常 高 于 
RRR, RIVARA ER k 将 式 (A6.4) RAY k RAR 
(A6.5), 则 有 方程 
Z(e,s)=0 (A6.6) 
对 * = co K p, 因而 在 方程 (A6.6) PS :二 za*， 从 而 给 出 
Z(ess) = Z(p, zt)= 2 'Q(p, 2)=0 (A6.7) 
其 中 :是 式 (A6.6) 中 :的 最 高 次 数 ， 方 程 


alp, 0)=0 (A6.8) 

可 给 出 多 变量 的 枢 点 值 p。 

例 

考虑 开 环 增 益 矩 阵 

1 
SN) 
—s— 11st — 295+ g2 一 202 十 355 十 70 ] 
412 一 :一 91 332 — 170s 十 118 

它 具 有 -- 特 征 方程 


A(g,s) = (s'— è + 28 — 25s + 29)g? 
—(—8 + 22° — 199s + 210g 
+ (—33s + 594) = 0 
应 用 第 6.2 节 描 述 的 方法 ,我们 发 现 当 4 > co 时 特征 频 
率 轨 线 的 无 穷 分 支 由 下 述 近 似 给 出 : 
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szk 和 sx +/ (334) 
为 寻求 对 应 二 阶 Butterworth 和 模型 j/ (334) 的 枢 点 , 令 


s=p+vV(—33k) 
以 便 给 出 
k= (s = ey 
一 33 


由 特征 方程 A(g, s) = 0 则 得 到 
T(s,k) = (st — $ + 2s? — 25s + 29) 
+ (—s' + 228 — 199s + 210)k 
+ (一 335 + 594)K? = 0 
在 其 中 替换 掉 则 有 
Z(p,s) = 33°%(s* — s + 2° — 25s + 29) 
— 33(—s' + 22° — 199s + 210) 
x (—2ps + ø”) — 335°(22s? — 199s + 210) 
+ (一 33* + 594)(—4s'p + 670’ — 4s + p') 
+ 5945' = 0 
sar", WA 
Oe, 2) = 331 — z + 22? — 252° + 292") 
一 33( 一 1 十 22z —.199z? + 210z?)(—20 + pzz) 
一 33(22 一 1992 + 21027) + (—33 + 594z) 
X( 一 4o + 6zp? — 42’p' + zp) + 594 = 0 
因此 
Q(p,0) = 33? — 66p — 33 X 22 + 33 
xX 40 + 594 =0 


由 此 可 以 求 得 
p 一 一 14.5 
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设 从 开 环 增益 矩阵 G(s) 得 到 一 个 用 复 增益 变量 8 与 复 
频率 变量 :的 形 如 


FSf(g,5)=0 (A7.1) 
的 代数 方程 。 考 虑 下 述 由 五 的 微 商 形成 的 方程 : 
aF\ _ 
its (a7) 
与 
OF 
ai (a73) 


同时 满足 方程 (A7.1 与 2) 的 s Lecce 上 的 分 支点 ; 
而 同时 满足 方程 (A7.1 与 3) 的 &8 的 值 是 增益 曲面 上 的 分 支 


点 . 
aes 的 函数 , 则 F 关 于: 的 全 微分 为 
SF ) 4s 
= = 18 +, | s ds AN 
或 交换 之 , 考 起 :为 8 的 函数 , 则 有 
dF _ (OF OF\ ds. 
dg \ Og i} aa or dg MAG) 


由 方程 (A7.1), 将 恒 为 零 ,因此 其 全 微 商 必 为 零 , 于 是 由 方 
程 (A7.4 与 A7.5) WA 


/ ƏF\ F Xy 
Bite ae A716) 
与 
OF /OF\ d > 
(Ca) aes ), Fa VATI) 


若 在 增益 曲面 上 有 一 分 支点 , 则 方程 (A7.3) Ri 于 是 
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方程 (A7.6) & 


OF d 
一 | =0 LE 一 0 
a is 或 ds 


E 


由 此 推出 在 频率 平面 上 有 一 个 分 支点 或 一 个 驻 点 对 应 增益 曲 
面 上 的 分 支点 。 

反之 , 若 在 频率 曲面 上 有 一 分 支点 , 则 方程 (A7.2) 成 立 ， 
于 是 由 方程 (A7.7) 有 


BF Os \ _ 

(外 ), 一 ?或 (多)= 
这 可 推出 在 增益 曲面 上 有 一 个 分 支点 或 有 一 个 驻 点 对 应 频率 
曲面 上 的 分 支点 . 
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割 口 (在 复 平面 上 ) Cuts (in complex plan) 33,36, 118 
循环 未 Cyclical system 86 
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十 三 画 
解析 开拓 Analytic continuation 33,86 
解析 函数 Analytic function 34,86 
WE Decoupling 110 
频率 曲面 ”Frequency surface 34,36 
辐 角 原理 Principle of the Argument 3,45,50,69,123 
输入 -输出 稳定 性 ”lnput-output stability g 
零点 引 理 Zero lemma 27 
堆 点 多 项 式 法 则 Zero polynomial rule 28 
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HEHE Stability 9,13 


十 五 画 


增益 曲面 ”Gaio surface 38,42 HE 
Butterworth 模型 Burterworth patteros 85 
-加 曲面 g-surface 38 

Jordan 正则 型 Jordan canonical form 10,21 
Kalman-Bucy 滤波 器 Kalman-Bucy filter 1 
Laplace 变换 Laplace transforms 7 

Newton 图 Newton diagram 81,84 

Nyquist D 回 线 Nyquist D-contour 41,51 
Riemann 数 球 Riemann number sphere 124 
Riemann 曲面 ”Riemann surface 3,31,32 
Smith-McMillan 型 Smith-McMillao form 25,62 
s 曲面 s-surface 34 
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